Researcher profile

Sagar Chakraborty

Sagar Chakraborty contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

23 published item(s)

preprint2026arXiv

Oscillatory evolutionarily stable state and limit cycle in replicator dynamics

The idea of evolutionarily stable state (ESS) of a population is a cornerstone of evolutionary game theory; moreover, it coincides with the game-theoretic concept of Nash equilibrium. Such a state corresponds to a strategy adopted by the population such that a rare mutant strategy cannot invade the population. In parallel, the dynamical formulation of evolutionary game theory -- particularly through replicator dynamics embodying the tenet of survival of the fittest -- provides a framework for modelling frequency-dependent selection over time. While it is well known that an ESS corresponds to stable fixed point in replicator dynamics, the evolutionary game-theoretic characterization of limit cycles is unknown. Here we fill this lacuna by defining oscillatory ESS (OESS) which we prove to be a stable limit cycle. We also show when an OESS is unique and if there are multiple OESSes, then what their locations are in the phase space.

preprint2022arXiv

Dynamic Scaling in Rotating Turbulence: A Shell Model Study

We investigate the scaling form of appropriate time-scales extracted from time-dependent correlation functions in rotating, turbulent flows. In particular, we obtain precise estimates of the dynamic exponents $z_p$, associated with the time-scales, and their relation with the more commonly measured equal-time exponents $ζ_p$. These theoretical predictions, obtained by using the multifractal formalism, are validated through extensive numerical simulations of a shell model for such rotating flows.

preprint2022arXiv

Eco-evolutionary games for harvesting self-renewing common resource: Effect of growing harvester population

The tragedy of the commons (TOC) is a ubiquitous social dilemma witnessed in interactions between a population of living entities and shared resources available to them: The individuals in the population tend to selfishly overexploit a common resource as it is arguably the rational choice, or in case of non-human beings, it may be an evolutionarily uninvadable action. How to avert the TOC is a significant problem related to the conservation of resources. It is not hard to envisage situations where the resource could be self-renewing and the size of the population may be dependent on the state of the resource through the fractions of the population employing different exploitation rates. If the self-renewal rate of the resource lies between the maximum and the minimum exploitation rates, it is not a priori obvious under what conditions the TOC can be averted. In this paper, we address this question analytically and numerically using the setup of an evolutionary game theoretical replicator equation that models the Darwinian tenet of natural selection. Through the replicator equation, while we investigate how a population of replicators exploit the shared resource, the latter's dynamical feedback on the former is also not ignored. We also present a transparent bottom-up derivation of the game-resource feedback model to facilitate future studies on the stochastic effects on the findings presented herein.

preprint2021arXiv

Cooperators overcome migration dilemma through synchronization

Synchronization, cooperation, and chaos are ubiquitous phenomena in nature. In a population composed of many distinct groups of individuals playing the prisoner's dilemma game, there exists a migration dilemma: No cooperator would migrate to a group playing the prisoner's dilemma game lest it should be exploited by a defector; but unless the migration takes place, there is no chance of the entire population's cooperator-fraction to increase. Employing a randomly rewired coupled map lattice of chaotic replicator maps, modelling replication-selection evolutionary game dynamics, we demonstrate that the cooperators -- evolving in synchrony -- overcome the migration dilemma to proliferate across the population when altruism is mildly incentivized making few of the demes play the leader game.

preprint2021arXiv

Deciphering chaos in evolutionary games

Discrete-time replicator map is a prototype of evolutionary selection game dynamical models that have been very successful across disciplines in rendering insights into the attainment of the equilibrium outcomes, like the Nash equilibrium and the evolutionarily stable strategy. By construction, only the fixed point solutions of the dynamics can possibly be interpreted as the aforementioned game-theoretic solution concepts. Although more complex outcomes like chaos are omnipresent in the nature, it is not known to which game-theoretic solutions they correspond. Here we construct a game-theoretic solution that is realized as the chaotic outcomes in the selection monotone game dynamic. To this end, we invoke the idea that in a population game having two-player--two-strategy one-shot interactions, it is the product of the fitness and the heterogeneity (the probability of finding two individuals playing different strategies in the infinitely large population) that is optimized over the generations of the evolutionary process.

preprint2021arXiv

Effect of chaotic agent dynamics on coevolution of cooperation and synchronization

The effect of the chaotic dynamical states of the agents on the coevolution of cooperation and synchronization in a structured population of the agents remains unexplored. With a view to gaining insights into this problem, we construct a coupled map lattice of the paradigmatic chaotic logistic map by adopting the Watts--Strogatz network algorithm. The map models the agent's chaotic state dynamics. In the model, an agent benefits by synchronizing with its neighbours and in the process of doing so, it pays a cost. The agents update their strategies (cooperation or defection) by using either a stochastic or a deterministic rule in an attempt to fetch themselves higher payoffs than what they already have. Among some other interesting results, we find that beyond a critical coupling strength, that increases with the rewiring probability parameter of the Watts--Strogatz model, the coupled map lattice is spatiotemporally synchronized regardless of the rewiring probability. Moreover, we observe that the population does not desynchronize completely -- and hence finite level of cooperation is sustained -- even when the average degree of the coupled map lattice is very high. These results are at odds with how a population of the non-chaotic Kuramoto oscillators as agents would behave. Our model also brings forth the possibility of the emergence of cooperation through synchronization onto a dynamical state that is a periodic orbit attractor.

preprint2021arXiv

Evolutionary dynamics of the delayed replicator-mutator equation: Limit cycle and cooperation

Game theory deals with strategic interactions among players and evolutionary game dynamics tracks the fate of the players' populations under selection. In this paper, we consider the replicator equation for two-player-two-strategy games involving cooperators and defectors. We modify the equation to include the effect of mutation and also delay that corresponds either to the delayed information about the population state or in realizing the effect of interaction among players. By focusing on the four exhaustive classes of symmetrical games -- the Stag Hunt game, the Snowdrift game, the Prisoners' Dilemma game, and the Harmony game -- we analytically and numerically analyze delayed replicator-mutator equation to find the explicit condition for the Hopf bifurcation bringing forth stable limit cycle. The existence of the asymptotically stable limit cycle imply the coexistence of the cooperators and the defectors; and in the games, where defection is a stable Nash strategy, a stable limit cycle does provide a mechanism for evolution of cooperation. We find that while mutation alone can never lead to oscillatory cooperation state in two-player-two-strategy games, the delay can change the scenario. On the other hand, there are situations when delay alone cannot lead to the Hopf bifurcation in the absence of mutation in the selection dynamics.

preprint2021arXiv

Periodic Orbit can be Evolutionarily Stable: Case Study of Discrete Replicator Dynamics

In evolutionary game theory, it is customary to be partial to the dynamical models possessing fixed points so that they may be understood as the attainment of evolutionary stability, and hence, Nash equilibrium. Any show of periodic or chaotic solution is many a time perceived as a shortcoming of the corresponding game dynamic because (Nash) equilibrium play is supposed to be robust and persistent behaviour, and any other behaviour in nature is deemed transient. Consequently, there is a lack of attempt to connect the non-fixed point solutions with the game theoretic concepts. Here we provide a way to render game theoretic meaning to periodic solutions. To this end, we consider a replicator map that models Darwinian selection mechanism in unstructured infinite-sized population whose individuals reproduce asexually forming non-overlapping generations. This is one of the simplest evolutionary game dynamic that exhibits periodic solutions giving way to chaotic solutions (as parameters related to reproductive fitness change) and also obeys the folk theorems connecting fixed point solutions with Nash equilibrium. Interestingly, we find that a modified Darwinian fitness -- termed heterogeneity payoff -- in the corresponding population game must be put forward as (conventional) fitness times the probability that two arbitrarily chosen individuals of the population adopt two different strategies. The evolutionary dynamics proceeds as if the individuals optimize the heterogeneity payoff to reach an evolutionarily stable orbit, should it exist. We rigorously prove that a locally asymptotically stable period orbit must be heterogeneity stable orbit -- a generalization of evolutionarily stable state.

preprint2021arXiv

Replicator equations induced by microscopic processes in nonoverlapping population playing bimatrix games

This paper is concerned with exploring the microscopic basis for the discrete versions of the standard replicator equation and the adjusted replicator equation. To this end, we introduce frequency-dependent selection -- as a result of competition fashioned by game-theoretic consideration -- into the Wright--Fisher process, a stochastic birth-death process. The process is further considered to be active in a generation-wise nonoverlapping finite population where individuals play a two-strategy bimatrix population game. Subsequently, connections among the corresponding master equation, the Fokker--Planck equation, and the Langevin equation are exploited to arrive at the deterministic discrete replicator maps in the limit of infinite population size.

preprint2021arXiv

Weight of fitness deviation governs strict physical chaos in replicator dynamics

Replicator equation -- a paradigm equation in evolutionary game dynamics -- mathematizes the frequency dependent selection of competing strategies vying to enhance their fitness (quantified by the average payoffs) with respect to the average fitnesses of the evolving population under consideration. In this paper, we deal with two discrete versions of the replicator equation employed to study evolution in a population where any two players, interaction is modeled by a two-strategy symmetric normal-form game. There are twelve distinct classes of such games, each typified by a particular ordinal relationship among the elements of the corresponding payoff matrix. Here we find the sufficient conditions for the existence of asymptotic solutions of the replicator equations such that the solutions -- fixed points, periodic orbits, and chaotic trajectories -- are all strictly physical, meaning that the frequency of any strategy lies inside the closed interval zero to one at all times. Thus, we elaborate which of the twelve types of games are capable of showing meaningful physical solutions and for which of the two types of replicator equation. Subsequently, we introduce the concept of the weight of fitness deviation that is the scaling factor in a positive affine transformation connecting two payoff matrices such that the corresponding one-shot games have exactly same Nash equilibria and evolutionary stable states. The weight also quantifies how much the excess of fitness of a strategy over the average fitness of the population affects the per capita change in the frequency of the strategy. Intriguingly, the weight's variation is capable of making the Nash equilibria and the evolutionary stable states useless by introducing strict physical chaos in the replicator dynamics based on the normal-form game.

preprint2020arXiv

Bridging Inertial and Dissipation Range Statistics in Rotating Turbulence

We investigate the connection between the inertial range and the dissipation range statistics of rotating turbulence through detailed simulations of a helical shell model and a multifractal analysis. In particular, by using the latter, we find an explicit relation between the (anomalous) scaling exponents of equal-time structure functions in the inertial range in terms of the generalised dimensions associated with the energy dissipation rate. This theoretical prediction is validated by detailed simulations of a helical shell model for various strengths of rotation from where the statistics of dissipation rate, and thus the generalised dimensions, as well as the inertial range, in particular the anomalous scaling exponents, are extracted. Our work also underlines a surprisingly good agreement---such as in the spatial structure of the energy dissipation rates and the decrease in inertial range intermittency with increasing strengths of rotation---between solutions of the Navier--Stokes equation in a rotating frame with those obtained from low-dimensional, dynamical systems such as the shell model which are not explicitly anisotropic. Finally, we perform direct numerical simulations of the Navier--Stokes equation, with the Coriolis force incorporated, to confirm the robustness of the conclusions drawn from our multifractal and shell model studies.

preprint2020arXiv

Comprehending deterministic and stochastic occasional uncoupling induced synchronizations through each other

In this paper, we numerically study the stochastic and the deterministic occasional uncoupling methods of effecting identical synchronized states in low dimensional, dissipative, diffusively coupled, chaotic flows that are otherwise not synchronized when continuously coupled at the same coupling strength parameter. In the process of our attempt to understand the mechanisms behind the success of the occasional uncoupling schemes, we devise a hybrid between the transient uncoupling and the stochastic on-off coupling, and aptly name it the transient stochastic uncoupling---yet another stochastic occasional uncoupling method. Our subsequent investigation on the transient stochastic uncoupling allows us to surpass the effectiveness of the stochastic on-off coupling with very fast on-off switching rate. Additionally, through the transient stochastic uncoupling, we establish that the indicators quantifying the local contracting dynamics in the corresponding transverse manifold are generally not useful in finding the optimal coupling region of the phase space in the case of the deterministic transient uncoupling. In fact, we highlight that the autocorrelation function---a non-local indicator of the dynamics---of the corresponding response system's chaotic time-series dictates when the deterministic uncoupling could be successful. We illustrate all our heuristic results using a few well-known examples of diffusively coupled chaotic oscillators.

preprint2020arXiv

Giant Tunable Mechanical Nonlinearity in Graphene-Silicon Nitride Hybrid Resonator

High quality factor mechanical resonators have shown great promise in developing classical or quantum technologies. Simultaneously, progress has been made in developing controlled mechanical nonlinearity. Here we combine these two directions of progress in a single platform consisting of coupled Silicon Nitride (SiNx) and graphene mechanical resonators. We show that nonlinear response can be induced on a large area SiNx resonator mode and can be efficiently controlled by coupling it to a gate-tunable, freely suspended graphene mode. The induced nonlinear response of the hybrid modes, as measured on the SiNx resonator surface is giant, with one of the highest measured Duffing constants. We observe a novel phononic frequency comb which we use as an alternate validation of the measured values, along with numerical simulations which are in overall agreement with measurements.

preprint2020arXiv

Suppressing birhythmicity by parametrically modulating nonlinearity in limit cycle oscillators

Multirhythmicity, a form of multistability, in an oscillator is an intriguing phenomenon found across many branches of science. From an application point of view, while the multirhythmicity is sometimes desirable as it presents us with many possible coexisting stable oscillatory states to tap into, it can also be a nuisance because a random perturbation may make the system settle onto an unwanted stable state. Consequently, it is not surprising that there are many natural and artificial mechanisms available that can control the multirhythmicity. We propose in this paper the idea of incorporating parametric (periodic) modulation of the nonlinear damping in the limit cycle oscillators with a view to exciting resonance and antiresonance responses at particular angular driving frequencies, and controlling the resulting birhythmicity by changing the amplitude of the modulation. To this end, we employ analytical (perturbative) and numerical techniques on the van der Pol oscillator---a paradigmatic limit cycle system---having additional position dependent time delay term and its modified autonomous birhythmic version. We also bring the fact to the fore that introduction of delay---a commonly adopted method of controlling multirhythmicity---in such a system can sometimes bring forth unwanted birhythmicity; and interestingly, our method of controlling birhythmicity through periodic modulation can suppress such a delay induced birhythmic response.

preprint2013arXiv

The Stability of Weakly Collisional Plasmas with Thermal and Composition Gradients

Over the last decade, substantial efforts have been devoted to understanding the stability properties, transport phenomena, and long-term evolution of weakly-collisional, magnetized plasmas which are stratified in temperature. These studies have improved our understanding of the physics governing the intra-cluster medium (ICM), but assumed that ICM is a homogeneous. This, however, might not be a good approximation if heavy elements sediment in the inner region of the galaxy cluster. In this paper, we analyze the stability of a weakly-collisional, magnetized plane-parallel atmosphere which is stratified in both temperature and composition. This allows us to discuss for the first time the dynamics of weakly-collisional environments where heat conduction, momentum transport, and ion-diffusion are anisotropic with respect to the direction of the magnetic field. We show that, depending on the relative signs and magnitudes of the gradients in the temperature and the mean molecular weight, the plasma can be subject to a wide variety of unstable modes which include modifications to the magnetothermal instability (MTI), the heat-flux-driven buoyancy instability (HBI), and overstable gravity modes previously studied in homogeneous media. We discuss the astrophysical implications of our findings for a representative galaxy cluster where helium has sedimented. Our findings suggest that the core insulation that results from the magnetic field configurations that arise as a natural consequence of the HBI, which would be MTI stable in a homogeneous medium, could be alleviated if the mean molecular weight gradient is steep enough, i.e., $(\nabla μ)/μ> (\nabla T)/T$. This study constitutes a first step toward understanding the interaction between magnetic turbulence and the diffusion of heavy elements in the ICM. (abridged)

preprint2012arXiv

Nelkin scaling for the Burgers equation and the role of high-precision calculations

Nelkin scaling, the scaling of moments of velocity gradients in terms of the Reynolds number, is an alternative way of obtaining inertial-range information. It is shown numerically and theoretically for the Burgers equation that this procedure works already for Reynolds numbers of the order of 100 (or even lower when combined with a suitable extended self-similarity technique). At moderate Reynolds numbers, for the accurate determination of scaling exponents, it is crucial to use higher than double precision. Similar issues are likely to arise for three-dimensional Navier--Stokes simulations.

preprint2012arXiv

Oscillators with Parametrically Excited Nonlinearity: Resonance, Anti-resonance and Switch

We discover presence of a hitherto unexplored type of resonance in a parametrically excited Van der Pol oscillator. The oscillator also possesses a state of anti-resonance. In the weak non-linear limit, we explain how to practically get a complete picture of different states of limiting oscillations present in the oscillator when the non-linear term therein is excited by an arbitrary $2π$ periodic function of time. We also illustrate how two such oscillators can be coupled to behave like a two-state switch allowing an sharp change of value of amplitude for stable oscillations from one constant to another.

preprint2011arXiv

The Stability of Dilute Plasmas with Thermal and Composition Gradients. I. The Slow Conduction Limit: Overstable Gravity Modes

We analyze the stability of a dilute plasma with thermal and composition gradients in the limit where conduction is slow compared to the dynamical timescale. We find necessary and sufficient conditions for stability when the background magnetic field is either parallel or perpendicular to the thermal and composition gradients that are parallel to the gravitational field. We provide approximate solutions for all the relevant modes involved, which are driven by gravity, conduction, and diffusion. We discuss the astrophysical implications of our findings for a representative galaxy cluster where helium has sedimented.

preprint2010arXiv

Three-dimensional turbulent relative dispersion by the Gledzer-Ohkitani-Yamada shell model

We study pair dispersion in a three-dimensional incompressible high Reynolds number turbulent flow generated by Fourier transforming the dynamics of the Gledzer-Ohkitani-Yamada (GOY) shell model into real space. We show that GOY shell model can successfully reproduce both the Batchelor and the Richardson-Obukhov regimes of turbulent relative dispersion. We also study how the crossover time scales with the initial separations of a particle pair and compare it to the prediction by Batchelor.

preprint2009arXiv

Extended Self Similarity works for the Burgers equation and why

Extended Self-Similarity (ESS), a procedure that remarkably extends the range of scaling for structure functions in Navier--Stokes turbulence and thus allows improved determination of intermittency exponents, has never been fully explained. We show that ESS applies to Burgers turbulence at high Reynolds numbers and we give the theoretical explanation of the numerically observed improved scaling at both the infrared and ultraviolet end, in total a gain of about three quarters of a decade: there is a reduction of subdominant contributions to scaling when going from the standard structure function representation to the ESS representation. We conjecture that a similar situation holds for three-dimensional incompressible turbulence and suggest ways of capturing subdominant contributions to scaling.

preprint2009arXiv

Why does the Sun's torsional oscillation begin before the sunspot cycle?

Although the Sun's torsional oscillation is believed to be driven by the Lorentz force associated with the sunspot cycle, this oscillation begins 2--3 yr before the sunspot cycle. We provide a theoretical explanation of this with the help of a solar dynamo model having a meridional circulation penetrating slightly below the bottom of the convection zone, because only in such dynamo models the strong toroidal field forms a few years before the sunspot cycle and at a higher latitude.