BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Marcus A. M. de Aguiar

Marcus A. M. de Aguiar contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
Populations and Evolutionnlin.AOquant-phmath-phmath.MPnlin.PSphysics.class-phphysics.data-anphysics.soc-phq-fin.GN

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
18works
0followers
10topics
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

18 published item(s)

preprint2026arXiv

Interplay of sync and swarm: Theory and application of swarmalators

Swarmalators, entities that combine the properties of swarming particles with synchronized oscillations, represent a novel and growing area of research in the study of collective behavior. This review provides a comprehensive overview of the current state of swarmalator research, focusing on the interplay between spatial organization and temporal synchronization. After a brief introduction to synchronization and swarming as separate phenomena, we discuss the various mathematical models that have been developed to describe swarmalator systems, highlighting the key parameters that govern their dynamics. The review also discusses the emergence of complex patterns, such as clustering, phase waves, and synchronized states, and how these patterns are influenced by factors such as interaction range, coupling strength, and frequency distribution. Recently, some minimal models were proposed that are solvable and mimic real-world phenomena. The effect of predators in the swarmalator dynamics is also discussed. Finally, we explore potential applications in fields ranging from robotics to biological systems, where understanding the dual nature of swarming and synchronization could lead to innovative solutions. By synthesizing recent advances and identifying open challenges, this review aims to provide a foundation for future research in this interdisciplinary field.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

The multi-allelic Moran process as a multi-zealot voter model: exact results and consequences for diversity thresholds

The Moran process is a foundational model of genetic drift and mutation in finite populations. In its standard two-allele form with population size $n$, allele counts, and hence allele frequencies, change through stochastic replacement and mutation, yet converge to a stationary distribution. This distribution undergoes a qualitative transition at the \emph{critical mutation rate} $μ_c=1/(2n)$: at $μ=μ_c$ it is exactly uniform, so that the probability of observing $k$ copies of allele~1 (and $n-k$ of allele~2) is $π(k)=1/(n+1)$ for $k=0,\dots,n$. For $μ<μ_c$ diversity is low: the stationary distribution places most of its mass near $k=0$ and $k=n$, and the population is therefore typically dominated by one allele. For $μ>μ_c$, on the other hand, diversity is high: the distribution concentrates around intermediate values, so that both alleles are commonly present at comparable frequencies. Recently, the two-allele Moran process was shown to be exactly equivalent to the voter model with two candidates and $α_1$ and $α_2$ committed voters (\emph{zealots}) in a population of $n+α_1+α_2$, where mutation is played by zealot influence. Here we extend this equivalence to multiple alleles and multiple candidates. Using the mapping, we derive the exact stationary distribution of allele counts for well-mixed populations with an arbitrary number $m$ of alleles, and obtain the critical mutation rate $μ_c = 1/(m+2n-2)$, which depends explicitly on $m$. We then analyze the Moran process on randomly connected populations and show that both the stationary distribution and $μ_c$ are invariant to network structure and coincide with the well-mixed results. Finally, simulations on general network topologies show that structural heterogeneity can substantially reshape the stationary allele distribution and, consequently, the level of genetic diversity.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Laplacian Matrices for Extremely Balanced and Unbalanced Phylogenetic Trees

Phylogenetic trees are important tools in the study of evolutionary relationships between species. Measures such as the index of Sackin, Colless, and Total Cophenetic have been extensively used to quantify tree balance, one key property of phylogenies. Recently a new proposal has been introduced, based on the spectrum of the Laplacian matrix associated with the tree. In this work, we calculate the Laplacian matrix analytically for two extreme cases, corresponding to fully balanced and fully unbalanced trees. For maximally balanced trees no closed form for the Laplacian matrix was derived, but we present an algorithm to construct it. We show that Laplacian matrices of fully balanced trees display self-similar patterns that result in highly degenerated eigenvalues. Degeneracy is the main signature of this topology, since it is totally absent in fully unbalanced trees. We also establish some analytical and numerical results about the largest eigenvalue of Laplacian matrices for these topologies.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Measuring Tree Balance with Normalized Tree Area

The study of species organization and their clustering by genetic or phenotypic similarity is carried out with the tools of phylogenetic trees. An important structural property of phylogenetic trees is the balance, which measures how taxa are distributed among clades. Tree balance can be measured using indices such as the Sackin ($S$) and the Total Cophenetic ($Φ$), which are based on the distance between nodes of the tree and its root. Here, we propose a new metric for tree balance, $\bar{d}$, the Area per Pair (APP) of the tree, which is a re-scaled version of the so called tree area. We compute $\bar{d}$ for the rooted caterpillar and maximally balanced trees and we also obtain exact formulas for its expected value and variance under the Yule model. The variance of APP for Yule trees has the remarkable property of converging to an asymptotic constant value for large trees. We compare the Sackin, Total Cophenetic and APP indices for hundreds of empirical phylogenies and show that APP represents the observed distribution of tree balances better than the two other metrics.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantifying the effects of quarantine using an IBM SEIR model on scalefree networks

The COVID-19 pandemic led several countries to resort to social distancing, the only known way to slow down the spread of the virus and keep the health system under control. Here we use an individual based model (IBM) to study how the duration, start date and intensity of quarantine affect the height and position of the peak of the infection curve. We show that stochastic effects, inherent to the model dynamics, lead to variable outcomes for the same set of parameters, making it crucial to compute the probability of each result. To simplify the analysis we divide the outcomes in only two categories, that we call {best and worst scenarios. Although long and intense quarantine is the best way to end the epidemic, it is very hard to implement in practice. Here we show that relatively short and intense quarantine periods can also be very effective in flattening the infection curve and even killing the virus, but the likelihood of such outcomes are low. Long quarantines of relatively low intensity, on the other hand, can delay the infection peak and reduce its size considerably with more than 50% probability, being a more effective policy than complete lockdown for short periods.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Sympatric speciation based on pure assortative mating

Although geographic isolation has been shown to play a key role in promoting reproductive isolation, it is now believed that speciation can also happen in sympatry and with considerable gene flow. Here we present a model of sympatric speciation based on assortative mating that does not require a genetic threshold for reproduction, i.e., that does not directly associate genetic differences between individuals with reproductive incompatibilities. In the model individuals mate with the most similar partner in their pool of potential mates, irrespective of how dissimilar it might be. We show that assortativity alone can lead to the formation of clusters of genetically similar individuals. The absence of a minimal genetic similarity for mating implies the constant generation of hybrids and brings up the old problem of species definition. Here, we define species based on clustering of genetically similar individuals but allowing genetic flow among different species. We show that the results obtained with the present model are in good agreement with empirical data, in which different species can still reproduce and generate hybrids.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Diverse interactions and ecosystem engineering stabilize community assembly

The complexity of an ecological community can be distilled into a network, where diverse interactions connect species in a web of dependencies. Species interact not only with each other but indirectly through environmental effects, however the role of these ecosystem engineers has not yet been considered in models of ecological networks. Here we explore the dynamics of ecosystem assembly, where the colonization and extinction of species within a community depends on the constraints imposed by trophic, service, and engineering dependencies. We show that our assembly model reproduces many key features of ecological systems, such as the role of generalists during assembly, realistic maximum trophic levels, and increased nestedness with higher frequencies of mutualisms. We find that ecosystem engineering has large and nonlinear effects on extinction rates, facilitating robustness by creating niche space, but at the same time increasing the magnitude of extinction cascades. We emphasize the importance of redundancies in engineered effects and show that such redundancy lowers the barriers to colonization, promoting community diversity. Together, our results suggest that ecological engineers may enhance community diversity while increasing persistence by facilitating colonization and limiting competitive exclusion.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Modular structure in C. elegans neural network and its response to external localized stimuli

Synchronization plays a key role in information processing in neuronal networks. Response of specific groups of neurons are triggered by external stimuli, such as visual, tactile or olfactory inputs. Neurons, however, can be divided into several categories, such as by physical location, functional role or topological clustering properties. Here we study the response of the electric junction C. elegans network to external stimuli using the partially forced Kuramoto model and applying the force to specific groups of neurons. Stimuli were applied to topological modules, obtained by the ModuLand procedure, to a ganglion, specified by its anatomical localization, and to the functional group composed of all sensory neurons. We found that topological modules do not contain purely anatomical groups or functional classes, corroborating previous results, and that stimulating different classes of neurons lead to very different responses, measured in terms of synchronization and phase velocity correlations. In all cases, however, the modular structure hindered full synchronization, protecting the system from seizures. More importantly, the responses to stimuli applied to topological and functional modules showed pronounced patterns of correlation or anti-correlation with other modules that were not observed when the stimulus was applied to ganglia.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Synchronization and spatial patterns in forced swarmalators

Swarlamators are particles capable of synchronize and swarm. Here we study the effects produced by an external periodic stimulus over a system of swarmalators that move in two dimensions. When the particles are fixed and interact with equal strength (Kuramoto oscillators) their phases tend to synchronize and lock to the external stimulus if its intensity is sufficiently large. Here we show that in a system of swarmalators the force also shifts the phases and angular velocities leading to synchronization with the external frequency. However, the correlation between phase and spatial location decreases with the intensity of the force, going to zero in what appears to be a second order phase transition. In the regime of zero correlation the particles form a static symmetric circular distribution, following a simple model of aggregation. Interestingly, for intermediate values of the force intensity, a different pattern emerges, with the particles splitting in two clusters centered at opposite sides of the stimulus&#39; location, breaking the radial symmetry. The two-cluster pattern is stable and active, with the clusters slowly rotating around the source while exchanging particles.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Revealing biases in the sampling of ecological interaction networks

The structure of ecological interactions is commonly understood through analyses of interaction networks. However, these analyses may be sensitive to sampling biases in both the interactors (the nodes of the network) and interactions (the links between nodes). These issues may affect the accuracy of empirically constructed ecological networks. We explore the properties of sampled ecological networks by simulating large-scale ecological networks with predetermined topologies, and sampling them with different mathematical procedures. Several types of modular networks were generated, intended to represent a wide variety of communities that vary in size and types of ecological interactions. We sampled these networks with different sampling designs that may be encountered in field experiments. The observed networks generated by each sampling process were analyzed with respect to number and size of components. We show that the sampling effort needed to estimate underlying network properties depends both on the sampling design and on network topology. Networks with random or scale-free modules require more complete sampling compared to networks whose modules are nested or bipartite. Overall, the structure of nested modules was the easiest to detect, regardless of sampling design. Sampling according to species degree was consistently found to be the most accurate strategy to estimate network structure. Conversely, sampling according to module results in an accurate view of certain modules, but fails to provide a global picture of the underlying network. We recommend that these findings are incorporated into the design of projects aiming to characterize large networks of species interactions in the field, to reduce sampling biases. The software scripts developed to construct and sample networks are provided for further explorations of network structure and comparisons to real interaction data.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The $500.00 AAPL close: Manipulation or hedging? A quantitative analysis

Why do a market&#39;s prices move up or down? Claims about causes are made without actual information, and accepted or dismissed based upon poor or non-existent evidence. Here we investigate the price movements that ended with Apple stock closing at \$500.00 on January 18, 2013. There is a ready explanation for this price movement: market manipulation by those who sold stock options, who stood to directly benefit from this closing price. Indeed, one web commentator predicted this otherwise unlikely event publicly. This explanation was subsequently dismissed by press articles that claim that stock prices end near such round numbers based upon legitimate hedging activity. But how can we know? We show that the accepted model that points to hedging as the driving cause of prices is not quantitatively consistent with the price movement on that day. The price moved upward too quickly over a period in which the hedgers&#39; position would require selling rather than buying. Under these conditions hedgers would have driven the price away from the strike price rather than toward it. We also show that a long published theory of the role of hedging is incomplete mathematically, and that the correct theory results in much weaker price movements. This evidence substantially weakens the case of those who claim hedging as cause of anomalous market price movements. The explanation that market manipulation is responsible for the final close cannot be dismissed based upon unsubstantiated, even invalid, hedging claims. Such proffered explanations shield potential illegal activity from further inquiry even though the claims behind those explanations have not been demonstrated.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Evolutionary consequences of assortativeness in haploid genotypes

We study the evolution of allele frequencies in a large population where random mating is violated in a particular way that is related to recent works on speciation. Specifically, we consider non-random encounters in haploid organisms described by biallelic genes at two loci and assume that individuals whose alleles differ at both loci are incompatible. We show that evolution under these conditions leads to the disappearance of one of the alleles and substantially reduces the diversity of the population. The allele that disappears, and the other allele frequencies at equilibrium, depend only on their initial values, and so does the time to equilibration. However, certain combinations of allele frequencies remain constant during the process, revealing the emergence of strong correlation between the two loci promoted by the epistatic mechanism of incompatibility. We determine the geometrical structure of the haplotype frequency space and solve the dynamical equations, obtaining a simple rule to determine equilibrium solution from the initial conditions. We show that our results are equivalent to selection against double heterozigotes for a population of diploid individuals and discuss the relevance of our findings to speciation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Robustness Against Extinction by Stochastic Sex Determination in Small Populations

Sexually reproducing populations with small number of individuals may go extinct by stochastic fluctuations in sex determination, causing all their members to become male or female in a generation. In this work we calculate the time to extinction of isolated populations with fixed number $N$ of individuals that are updated according to the Moran birth and death process. At each time step, one individual is randomly selected and replaced by its offspring resulting from mating with another individual of opposite sex; the offspring can be male or female with equal probability. A set of $N$ time steps is called a generation, the average time it takes for the entire population to be replaced. The number k of females fluctuates in time, similarly to a random walk, and extinction, which is the only asymptotic possibility, occurs when k=0 or k=N. We show that it takes only one generation for an arbitrary initial distribution of males and females to approach the binomial distribution. This distribution, however, is unstable and the population eventually goes extinct in 2^N/N generations. We also discuss the robustness of these results against bias in the determination of the sex of the offspring, a characteristic promoted by infection by the bacteria Wolbachia in some arthropod species or by temperature in reptiles.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Energy Dissipation Via Coupling With a Finite Chaotic Environment

We study the flow of energy between a harmonic oscillator (HO) and an external environment consisting of N two-degrees of freedom non-linear oscillators, ranging from integrable to chaotic according to a control parameter. The coupling between the HO and the environment is bilinear in the coordinates and scales with system size with the inverse square root of N. We study the conditions for energy dissipation and thermalization as a function of N and of the dynamical regime of the non-linear oscillators. The study is classical and based on single realization of the dynamics, as opposed to ensemble averages over many realizations. We find that dissipation occurs in the chaotic regime for a fairly small N, leading to the thermalization of the HO and environment a Boltzmann distribution of energies for a well defined temperature. We develop a simple analytical treatment, based on the linear response theory, that justifies the coupling scaling and reproduces the numerical simulations when the environment is in the chaotic regime.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Initial value representation for the SU(n) semiclassical propagator

The semiclassical propagator in the representation of SU(n) coherent states is characterized by isolated classical trajectories subjected to boundary conditions in a doubled phase space. In this paper we recast this expression in terms of an integral over a set of initial-valued trajectories. These trajectories are monitored by a filter that collects only the appropriate contributions to the semiclassical approximation. This framework is suitable for the study of bosonic dynamics in n modes with fixed total number of particles. We exemplify the method for a Bose-Einstein condensate trapped in a triple-well potential, providing a detailed discussion on the accuracy and efficiency of the procedure.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Semiclassical propagator for SU(n) coherent states

We present a detailed derivation of the semiclassical propagator in the SU(n) coherent state representation. In order to provide support for immediate physical applications, we restrict this work to the fully symmetric irreducible representations, which are suitable for the treatment of bosonic dynamics in n modes, considering systems with conservation of total particle number. The derivation described here can be easily extended to other classes of coherent states, thus representing an alternative approach to previously published methods.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The role of sex separation in neutral speciation

Neutral speciation mechanisms based on isolation by distance and sexual selection, termed topopatric, have recently been shown to describe the observed patterns of abundance distributions and species-area relationships. Previous works have considered this type of process only in the context of hermaphrodic populations. In this work we extend a hermaphroditic model of topopatric speciation to populations where individuals are explicitly separated into males and females. We show that for a particular carrying capacity speciation occurs under similar conditions, but the number of species generated decreases as compared to the hermaphroditic case. Evolution results in fewer species having more abundant populations.

0 citations0 reviewsNo code linkedOpen access