Researcher profile

Jose A. Cuesta

Jose A. Cuesta contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

10 published item(s)

preprint2010arXiv

Catastrophic regime shifts in model ecological communities are true phase transitions

Ecosystems often undergo abrupt regime shifts in response to gradual external changes. These shifts are theoretically understood as a regime switch between alternative stable states of the ecosystem dynamical response to smooth changes in external conditions. Usual models introduce nonlinearities in the macroscopic dynamics of the ecosystem that lead to different stable attractors among which the shift takes place. Here we propose an alternative explanation of catastrophic regime shifts based on a recent model that pictures ecological communities as systems in continuous fluctuation, according to certain transition probabilities, between different micro-states in the phase space of viable communities. We introduce a spontaneous extinction rate that accounts for gradual changes in external conditions, and upon variations on this control parameter the system undergoes a regime shift with similar features to those previously reported. Under our microscopic viewpoint we recover the main results obtained in previous theoretical and empirical work (anomalous variance, hysteresis cycles, trophic cascades). The model predicts a gradual loss of species in trophic levels from bottom to top near the transition. But more importantly, the spectral analysis of the transition probability matrix allows us to rigorously establish that we are observing the fingerprints, in a finite size system, of a true phase transition driven by background extinctions.

preprint2010arXiv

Huge progeny production during the transient of a quasi-species model of viral infection, reproduction and mutation

Eigen's quasi-species model describes viruses as ensembles of different mutants of a high fitness "master" genotype. Mutants are assumed to have lower fitness than the master type, yet they coexist with it forming the quasi-species. When the mutation rate is sufficiently high, the master type no longer survives and gets replaced by a wide range of mutant types, thus destroying the quasi-species. It is the so-called "error catastrophe". But natural selection acts on phenotypes, not genotypes, and huge amounts of genotypes yield the same phenotype. An important consequence of this is the appearance of beneficial mutations which increase the fitness of mutants. A model has been recently proposed to describe quasi-species in the presence of beneficial mutations. This model lacks the error catastrophe of Eigen's model and predicts a steady state in which the viral population grows exponentially. Extinction can only occur if the infectivity of the quasi-species is so low that this exponential is negative. In this work I investigate the transient of this model when infection is started from a small amount of low fitness virions. I prove that, beyond an initial regime where viral population decreases (and can go extinct), the growth of the population is super-exponential. Hence this population quickly becomes so huge that selection due to lack of host cells to be infected begins to act before the steady state is reached. This result suggests that viral infection may widespread before the virus has developed its optimal form.

preprint2010arXiv

Neutral networks of genotypes: Evolution behind the curtain

Our understanding of the evolutionary process has gone a long way since the publication, 150 years ago, of "On the origin of species" by Charles R. Darwin. The XXth Century witnessed great efforts to embrace replication, mutation, and selection within the framework of a formal theory, able eventually to predict the dynamics and fate of evolving populations. However, a large body of empirical evidence collected over the last decades strongly suggests that some of the assumptions of those classical models necessitate a deep revision. The viability of organisms is not dependent on a unique and optimal genotype. The discovery of huge sets of genotypes (or neutral networks) yielding the same phenotype --in the last term the same organism--, reveals that, most likely, very different functional solutions can be found, accessed and fixed in a population through a low-cost exploration of the space of genomes. The 'evolution behind the curtain' may be the answer to some of the current puzzles that evolutionary theory faces, like the fast speciation process that is observed in the fossil record after very long stasis periods.

preprint2010arXiv

Species assembly in model ecosystems, I: Analysis of the population model and the invasion dynamics

Recently we have introduced a simplified model of ecosystem assembly (Capitan et al., 2009) for which we are able to map out all assembly pathways generated by external invasions in an exact manner. In this paper we provide a deeper analysis of the model, obtaining analytical results and introducing some approximations which allow us to reconstruct the results of our previous work. In particular, we show that the population dynamics equations of a very general class of trophic-level structured food-web have an unique interior equilibrium point which is globally stable. We show analytically that communities found as end states of the assembly process are pyramidal and we find that the equilibrium abundance of any species at any trophic level is approximately inversely proportional to the number of species in that level. We also find that the per capita growth rate of a top predator invading a resident community is key to understand the appearance of complex end states reported in our previous work. The sign of these rates allows us to separate regions in the space of parameters where the end state is either a single community or a complex set containing more than one community. We have also built up analytical approximations to the time evolution of species abundances that allow us to determine, with high accuracy, the sequence of extinctions that an invasion may cause. Finally we apply this analysis to obtain the communities in the end states. To test the accuracy of the transition probability matrix generated by this analytical procedure for the end states, we have compared averages over those sets with those obtained from the graph derived by numerical integration of the Lotka-Volterra equations. The agreement is excellent.

preprint2010arXiv

Species assembly in model ecosystems, II: Results of the assembly process

In the companion paper of this set (Capitan and Cuesta, 2010) we have developed a full analytical treatment of the model of species assembly introduced in Capitan et al. (2009). This model is based on the construction of an assembly graph containing all viable configurations of the community, and the definition of a Markov chain whose transitions are the transformations of communities by new species invasions. In the present paper we provide an exhaustive numerical analysis of the model, describing the average time to the recurrent state, the statistics of avalanches, and the dependence of the results on the amount of available resource. Our results are based on the fact that the Markov chain provides an asymptotic probability distribution for the recurrent states, which can be used to obtain averages of observables as well as the time variation of these magnitudes during succession, in an exact manner. Since the absorption times into the recurrent set are found to be comparable to the size of the system, the end state is quickly reached (in units of the invasion time). Thus, the final ecosystem can be regarded as a fluctuating complex system where species are continually replaced by newcomers without ever leaving the set of recurrent patterns. The assembly graph is dominated by pathways in which most invasions are accepted, triggering small extinction avalanches. Through the assembly process, communities become less resilient (e.g., have a higher return time to equilibrium) but become more robust in terms of resistance against new invasions.

preprint2010arXiv

The struggle for space: Viral extinction through competition for cells

The design of protocols to suppress the propagation of viral infections is an enduring enterprise, especially hindered by limited knowledge of the mechanisms through which extinction of infection propagation comes about. We here report on a mechanism causing extinction of a propagating infection due to intraspecific competition to infect susceptible hosts. Beneficial mutations allow the pathogen to increase the production of progeny, while the host cell is allowed to develop defenses against infection. When the number of susceptible cells is unlimited, a feedback runaway co-evolution between host resistance and progeny production occurs. However, physical space limits the advantage that the virus can obtain from increasing offspring numbers, thus infection clearance may result from an increase in host defenses beyond a finite threshold. Our results might be relevant to better understand propagation of viral infections in tissues with mobility constraints, and the implications that environments with different geometrical properties might have in devising control strategies.

preprint2009arXiv

Phase diagram of a two-dimensional lattice gas model of a ramp system

Using Monte Carlo Simulation and fundamental measure theory we study the phase diagram of a two-dimensional lattice gas model with a nearest neighbor hard core exclusion and a next-to-nearest neighbors finite repulsive interaction. The model presents two competing ranges of interaction and, in common with many experimental systems, exhibits a low density solid phase, which melts back to the fluid phase upon compression. The theoretical approach is found to provide a qualitatively correct picture of the phase diagram of our model system.

preprint2009arXiv

Promotion of cooperation on networks? The myopic best response case

We address the issue of the effects of considering a network of contacts on the emergence of cooperation on social dilemmas under myopic best response dynamics. We begin by summarizing the main features observed under less intellectually demanding dynamics, pointing out their most relevant general characteristics. Subsequently we focus on the new framework of best response. By means of an extensive numerical simulation program we show that, contrary to the rest of dynamics considered so far, best response is largely unaffected by the underlying network, which implies that, in most cases, no promotion of cooperation is found with this dynamics. We do find, however, nontrivial results differing from the well-mixed population in the case of coordination games on lattices, which we explain in terms of the formation of spatial clusters and the conditions for their advancement, subsequently discussing their relevance to other networks.

preprint2009arXiv

Smectic and columnar ordering in length-polydisperse fluids of parallel hard cylinders

We apply a recently proposed density functional for mixtures of parallel hard cylinders, based on Rosenfeld's fundamental measure theory, to study the effect of length-polydispersity on the relative stability between the smectic and columnar liquid crystal phases.To this purpose we derive from this functional an expression for the direct correlation function and use it to perform a bifurcation analysis. We compare the results with those obtained with a second and a third virial approximation of this function. All three approximations lead to the same conclusion: there is a terminal polydispersity beyond which the smectic phase is less stable than the columnar phase. This result is in agreement with previous Monte Carlo simulations conducted on a freely rotating length-polydisperse hard spherocylinder fluid, although the theories always overestimate the terminal polydispersity because the nematic-columnar phase transition is first order and exhibits a wide coexistence gap. Both, the fundamental-measure functional and the third virial approximation, predict a metastable nematic-nematic demixing. Conversely, according to second virial approximation this demixing might be stable at high values of the polydispersity, something that is observed neither in simulations nor in experiments. The results of the fundamental-measure functional are quantitatively superior to those obtained from the other two approximations. Thus this functional provides a promising route to map out the full phase diagram of this system.

preprint1997arXiv

Dimensional crossover of the fundamental-measure functional for parallel hard cubes

We present a regularization of the recently proposed fundamental-measure functional for a mixture of parallel hard cubes. The regularized functional is shown to have right dimensional crossovers to any smaller dimension, thus allowing to use it to study highly inhomogeneous phases (such as the solid phase). Furthermore, it is shown how the functional of the slightly more general model of parallel hard parallelepipeds can be obtained using the zero-dimensional functional as a generating functional. The multicomponent version of the latter system is also given, and it is suggested how to reformulate it as a restricted-orientation model for liquid crystals. Finally, the method is further extended to build a functional for a mixture of parallel hard cylinders.