Source author record

Marcus A. M. de Aguiar

Marcus A. M. de Aguiar appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

25works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

25 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.

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.

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.

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.

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.

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.

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.

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.

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' location, breaking the radial symmetry. The two-cluster pattern is stable and active, with the clusters slowly rotating around the source while exchanging particles.

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.

preprint2016arXiv

Multiconfigurational quantum propagation with trajectory-guided generalized coherent states

A generalized version of the coupled coherent states method for coherent states of arbitrary Lie groups is developed. In contrast to the original formulation, which is restricted to frozen-Gaussian basis sets, the extended method is suitable for propagating quantum states of systems featuring diversified physical properties, such as spin degrees of freedom or particle indistinguishability. The approach is illustrated with simple models for interacting bosons trapped in double- and triple-well potentials, most adequately described in terms of SU(2) and SU(3) bosonic coherent states, respectively.

preprint2016arXiv

Voting Contagion

Social influence plays an important role in human behavior and decisions. The sources of influence can be generally divided into external, which are independent of social context, or as originating from peers, such as family and friends. An important question is how to disentangle the social contagion by peers from external influences. While a variety of experimental and observational studies provided insight into this problem, identifying the extent of social contagion based on large-scale observational data with an unknown network structure remains largely unexplored. By bridging the gap between the large-scale complex systems perspective of collective human dynamics and the detailed approach of the social sciences, we present a parsimonious model of social influence, and apply it to a central topic in political science -- elections and voting behavior. We provide an analytical expression of the county vote-share distribution in a two party system, which is in excellent agreement with 92 years of observed U.S. presidential election data. Analyzing the social influence topography over this period reveals an abrupt transition in the patterns of social contagion -- from low to high levels of social contagion. The results from our analysis reveal robust differences among regions of the United States in terms of their social influence index. In particular, we identify two regions of 'hot' and 'cold spots of social influence, each comprising states that are geographically close. These results suggest that social contagion effects are becoming more instrumental in shaping large scale collective political behavior, which is at the core of democratic societies.

preprint2015arXiv

Binary dynamics on star networks under external perturbations

We study a binary dynamical process that is a representation of the voter model with opinion makers. The process models an election with two candidates but can also describe the frequencies of a biallelic gene in a population or atoms with two spin orientations in a magnetic material. The system is represented by a network whose nodes have internal states labeled 0 or 1, and nodes that are connected can influence each other. The network is perturbed by a set of external nodes whose states are fixed in 0 or 1 and that can influence all nodes of the network. The fixed nodes play the role of opinion makers in the voter model, mutation rates in population genetics or temperature in a magnetic material. The quantity of interest is the probability $ρ(m,t)$ that $m$ nodes are in state 1 at time $t$. Here we study this process on star networks and compare the results with those obtained for networks that are fully connected. In both cases a transition from disordered to ordered equilibrium states is observed as the number of external fixed nodes becomes small, but it differs significantly between the two network topologies. For fully connected networks the probability distribution becomes uniform at the critical point, which is independent of the network size. For star networks, on the other hand, the equilibrium distribution $ρ(m)$ splits in two peaks, reflecting the two possible states of the central node. We obtain approximate analytical solutions that hold near the transition and that clarify the role of the central node in the process. We show that the different character of the transition also manifest itself in the magnetization of system, obtained in the limit of large $N$. Finally, extending the analysis to two star networks we compare our results with simulations in scalefree networks, detecting the presence of hubs.

preprint2015arXiv

Semiclassical Propagator in the Generalized Coherent-State Representation

A detailed derivation of the semiclassical propagator in the generalized coherent-state representation is performed by applying the saddle-point method to a path integral over the classical phase space. With the purpose of providing greater accessibility and applicability to the developed formalism, a brief review of the generalized concept of coherent states is presented, in which three examples of coherent-state sets are examined, namely, the canonical, spin, and SU(n) bosonic coherent states.

preprint2015arXiv

Topological structures in the Husimi flow

We study the topological properties of the quantum phase space current in the Husimi representation, focusing on the dynamical differences, induced by these properties, between the quantum and the classical flows. We show that the zeros of the Husimi function are stagnation points of the current and have a nonzero topological charge. Due to overall charge conservation, new stagnation points with opposite charge appear in pairs in the Husimi current and they have important roles in dynamical processes. As an example we show the topological effect of the zeros in the transmission rate of particle through a potential barrier.

preprint2014arXiv

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

Why do a market'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' 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.

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.

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.

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.

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.

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.

preprint2011arXiv

The Moran model as a dynamical process on networks and its implications for neutral speciation

In genetics the Moran model describes the neutral evolution of a bi-allelic gene in a population of haploid individuals subjected to mutations. We show in this paper that this model can be mapped into an influence dynamical process on networks subjected to external influences. The panmictic case considered by Moran corresponds to fully connected networks and can be completely solved in terms of hypergeometric functions. Other types of networks correspond to structured populations, for which approximate solutions are also available. This new approach to the classic Moran model leads to a relation between regular networks based on spatial grids and the mechanism of isolation by distance. We discuss the consequences of this connection for topopatric speciation and the theory of neutral speciation and biodiversity. We show that the effect of mutations in structured populations, where individuals can mate only with neighbors, is greatly enhanced with respect to the panmictic case. If mating is further constrained by genetic proximity between individuals, a balance of opposing tendencies take place: increasing diversity promoted by enhanced effective mutations versus decreasing diversity promoted by similarity between mates. Stabilization occurs with speciation via pattern formation. We derive an explicit relation involving the parameters characterizing the population that indicates when speciation is possible.

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.

preprint2002arXiv

Solving the boundary value problem for finite Kirchhoff rods

The Kirchhoff model describes the statics and dynamics of thin rods within the approximations of the linear elasticity theory. In this paper we develop a method, based on a shooting technique, to find equilibrium configurations of finite rods subjected to boundary conditions and given load parameters. The method consists in making a series of small changes on a {\it trial} solution satisfying the Kirchhoff equations but not necessarily the boundary conditions. By linearizing the differential equations around the trial solution we are able to push its end point to the desired position, step by step. The method is also useful to obtain configurations of rods with fixed end points but different mechanical parameters, such as tension, components of the moment or inhomogeneities.