Source author record

Yoichiro Mori

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

ResearcherUnclaimed source record

Cell Behavior math.AP Populations and Evolution cond-mat.soft math.DS physics.flu-dyn Tissues and Organs Biological Physics nlin.PS physics.chem-ph physics.soc-ph

Catalog footprint

What is connected

9works

11topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Long-Time Behavior of a PDE Replicator Equation for Multilevel Selection in Group-Structured Populations

In many biological systems, natural selection acts simultaneously on multiple levels of organization. This scenario typically presents an evolutionary conflict between the incentive of individuals to cheat and the collective incentive to establish cooperation within a group. Generalizing previous work on multilevel selection in evolutionary game theory, we consider a hyperbolic PDE model of a group-structured population, in which members within a single group compete with each other for individual-level replication; while the group also competes against other groups for group-level replication. We derive a threshold level of the relative strength of between-group competition such that defectors take over the population below the threshold while cooperation persists in the long-time population above the threshold. Under stronger assumptions on the initial distribution of group compositions, we further prove that the population converges to a steady state density supporting cooperation for between-group selection strength above the threshold. We further establish long-time bounds on the time-average of the collective payoff of the population, showing that the long-run population cannot outperform the payoff of a full-cooperator group even in the limit of infinitely-strong between-group competition. When the group replication rate is maximized by an intermediate level of within-group cooperation, individual-level selection casts a long shadow on the dynamics of multilevel selection: no level of between-group competition can erase the effects of the individual incentive to defect. We further extend our model to study the case of multiple types of groups, showing how the games that groups play can coevolve with the level of cooperation.

preprint2022arXiv

Selfish optimization and collective learning in populations

A selfish learner seeks to maximize their own success, disregarding others. When success is measured as payoff in a game played against another learner, mutual selfishness typically fails to produce the optimal outcome for a pair of individuals. However, learners often operate in populations, and each learner may have a limited duration of interaction with any other individual. Here, we compare selfish learning in stable pairs to selfish learning with stochastic encounters in a population. We study gradient-based optimization in repeated games like the prisoner's dilemma, which feature multiple Nash equilibria, many of which are suboptimal. We find that myopic, selfish learning, when distributed in a population via ephemeral encounters, can reverse the dynamics that occur in stable pairs. In particular, when there is flexibility in partner choice, selfish learning in large populations can produce optimal payoffs in repeated social dilemmas. This result holds for the entire population, not just for a small subset of individuals. Furthermore, as the population size grows, the timescale to reach the optimal population payoff remains finite in the number of learning steps per individual. While it is not universally true that interacting with many partners in a population improves outcomes, this form of collective learning achieves optimality for several important classes of social dilemmas. We conclude that naïve learning can be surprisingly effective in populations of individuals navigating conflicts of interest.

preprint2021arXiv

Evolution of prosocial behavior in multilayer populations

Human societies include diverse social relationships. Friends, family, business colleagues, and online contacts can all contribute to one's social life. Individuals may behave differently in different domains, but success in one domain may engender success in another. Here, we study this problem using multilayer networks to model multiple domains of social interactions, in which individuals experience different environments and may express different behaviors. We provide a mathematical analysis and find that coupling between layers tends to promote prosocial behavior. Even if prosociality is disfavored in each layer alone, multilayer coupling can promote its proliferation in all layers simultaneously. We apply this analysis to six real-world multilayer networks, ranging from the socio-emotional and professional relationships in a Zambian community, to the online and offline relationships within an academic University. We discuss the implications of our results, which suggest that small modifications to interactions in one domain may catalyze prosociality in a different domain.

preprint2020arXiv

Accuracy of slender body theory in approximating force exerted by thin fiber on viscous fluid

We consider the accuracy of slender body theory in approximating the force exerted by a thin fiber on the surrounding viscous fluid when the fiber velocity is prescribed. We term this the slender body inverse problem, as it is known that slender body theory converges to a well-posed PDE solution when the force is prescribed and the fiber velocity is unknown. From a PDE perspective, the slender body inverse problem is simply the Dirichlet problem for the Stokes equations, but from an approximation perspective, nonlocal slender body theory exhibits high wavenumber instabilities which complicate analysis. Here we consider two methods for regularizing the slender body approximation: spectral truncation and the $δ$-regularization of Tornberg and Shelley (2004). For a straight, periodic fiber with constant radius $ε>0$, we explicitly calculate the spectrum of the operator mapping fiber velocity to force for both the PDE and the approximations. We show that the spectrum of the original slender body approximation agrees closely with the PDE solution at low wavenumbers but differs at high frequencies, allowing us to define a truncated approximation with a wavenumber cutoff $\sim1/ε$. For both the truncated and $δ$-regularized approximations, we obtain similar convergence results to the PDE solution as $ε\to0$: a fiber velocity with $H^1$ regularity gives $O(ε)$ convergence, while a fiber velocity with at least $H^2$ regularity yields $O(ε^2)$ convergence. Moreover, we determine the dependence of the $δ$-regularized error estimate on the regularization parameter $δ$.

preprint2020arXiv

An error bound for the slender body approximation of a thin, rigid fiber sedimenting in Stokes flow

We investigate the motion of a thin rigid body in Stokes flow and the corresponding slender body approximation used to model sedimenting fibers. In particular, we derive a rigorous error bound comparing the rigid slender body approximation to the classical PDE for rigid motion in the case of a closed loop with constant radius. Our main tool is the slender body PDE framework established by the authors and D. Spirn in [18,19], which we adapt to the rigid setting.

preprint2016arXiv

Rhythmomimetic drug delivery: modeling, analysis and numerical simulation

We develop, analyse and numerically simulate a model of a prototype, glucose-driven, rhythmic drug delivery device, aimed at hormone therapies, and based on chemomechanical interaction in a polyelectrolyte gel membrane. The pH-driven interactions trigger volume phase transitions between the swollen and collapsed states of the gel.For a robust set of material parameters, we find a class of solutions of the governing system that oscillate between such states, and cause the membrane to rhythmically swell, allowing for transport of the drug, fuel and reaction products across it, and collapse, hampering all transport across it. The frequency of the oscillations can be adjusted so that it matches the natural frequency of the hormone to be released. The work is linked to extensive laboratory experimental studies of the device built by Siegel's team. The thinness of the membrane and its clamped boundary, together with the homogeneously held conditions in the experimental apparatus, justify neglecting spatial dependence on the fields of the problem. Upon identifying the forces and energy relevant to the system, and taking into account its dissipative properties, we apply Rayleigh's variational principle to obtain the governing equations. The material assumptions guarantee the monotonicity of the system and lead to the existence of a three dimensional limit cycle. By scaling and asymptotic analysis, this limit cycle is found to be related to a two-dimensional one that encodes the volume phase transitions of the model. The identification of the relevant parameter set of the model is aided by a Hopf bifurcation study of steady state solutions.

preprint2014arXiv

A Multidomain Model for Ionic Electrodiffusion and Osmosis with an Application to Cortical Spreading Depression

Ionic electrodiffusion and osmotic water flow are central processes in many physiological systems. We formulate a system of partial differential equations that governs ion movement and water flow in biological tissue. A salient feature of this model is that it satisfies a free energy identity, ensuring the thermodynamic consistency of the model. A numerical scheme is developed for the model in one spatial dimension and is applied to a model of cortical spreading depression, a propagating breakdown of ionic and cell volume homeostasis in the brain.

preprint2011arXiv

Mathematical Properties of Pump-Leak Models of Cell Volume Control and Electrolyte Balance

Homeostatic control of cell volume and intracellular electrolyte content is a fundamental problem in physiology and is central to the functioning of epithelial systems. These physiological processes are modeled using pump-leak models, a system of differential algebraic equations that describes the balance of ions and water flowing across the cell membrane. Despite their widespread use, very little is known about their mathematical properties. Here, we establish analytical results on the existence and stability of steady states for a general class of pump-leak models. We treat two cases. When the ion channel currents have a linear current-voltage relationship, we show that there is at most one steady state, and that the steady state is globally asymptotically stable. If there are no steady states, the cell volume tends to infinity with time. When minimal assumptions are placed on the properties of ion channel currents, we show that there is an asymptotically stable steady state so long as the pump current is not too large. The key analytical tool is a free energy relation satisfied by a general class of pump-leak models, which can be used as a Lyapunov function to study stability.

preprint2010arXiv

Asymptotic and bifurcation analysis of wave-pinning in a reaction-diffusion model for cell polarization

We describe and analyze a bistable reaction-diffusion (RD) model for two interconverting chemical species that exhibits a phenomenon of wave-pinning: a wave of activation of one of the species is initiated at one end of the domain, moves into the domain, decelerates, and eventually stops inside the domain, forming a stationary front. The second ("inactive") species is depleted in this process. This behavior arises in a model for chemical polarization of a cell by Rho GTPases in response to stimulation. The initially spatially homogeneous concentration profile (representative of a resting cell) develops into an asymmetric stationary front profile (typical of a polarized cell). Wave-pinning here is based on three properties: (1) mass conservation in a finite domain, (2) nonlinear reaction kinetics allowing for multiple stable steady states, and (3) a sufficiently large difference in diffusion of the two species. Using matched asymptotic analysis, we explain the mathematical basis of wave-pinning, and predict the speed and pinned position of the wave. An analysis of the bifurcation of the pinned front solution reveals how the wave-pinning regime depends on parameters such as rates of diffusion and total mass of the species. We describe two ways in which the pinned solution can be lost depending on the details of the reaction kinetics: a saddle-node or a pitchfork bifurcation.