BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

José A. Carrillo

José A. Carrillo appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.APmath.NANumerical AnalysisNeurons and Cognitioncond-mat.softcond-mat.stat-mechmath.PRphysics.comp-phphysics.flu-dynmath.DSmath.FAmath.OCnlin.AOnlin.PSPopulations and Evolution

Catalog footprint

What is connected

42works
15topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

42 published item(s)

preprint2026arXiv

Asymptotic stability of steady states for the compressible Navier-Stokes-Riesz system in the presence of vacuum

We consider a one-dimensional physical vacuum free boundary problem on the compressible Navier-Stokes-Riesz system for an attractive Riesz potential $|x|^{2s-1}/(2s-1)$ with $0<s<1/2$. It is proved that for the adiabatic constant $γ$ satisfying $2(1-s)<γ<1+2s/3$ under the additional condition that $3/8<s<1/2$, there exists a unique global-in-time strong solution. Specifically, we establish the Lyapunov-type stability of the compactly supported steady states in the Lagrangian coordinates and we also obtain the time rate of convergence for the strong solution to steady states with the same mass in weighted Sobolev spaces where the weights indicate the behavior of solutions near the vacuum free boundary. The difficulties and challenges in the proof are caused not only by the degeneracy due to the vacuum free boundary but also by the non-local feature of the Riesz potential.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Can the clocks tick together despite the noise? Stochastic simulations and analysis

The suprachiasmatic nucleus (SCN), also known as the circadian master clock, consists of a large population of oscillator neurons. Together, these neurons produce a coherent signal that drives the body's circadian rhythms. What properties of the cell-to-cell communication allow the synchronization of these neurons, despite a wide range of environmental challenges such as fluctuations in photoperiods? To answer that question, we present a mean-field description of globally coupled neurons modeled as Goodwin oscillators with standard Gaussian noise. Provided that the initial conditions of all neurons are independent and identically distributed, any finite number of neurons becomes independent and has the same probability distribution in the mean-field limit, a phenomenon called propagation of chaos. This probability distribution is a solution to a Vlasov-Fokker-Planck type equation, which can be obtained from the stochastic particle model. We study, using the macroscopic description, how the interaction between external noise and intercellular coupling affects the dynamics of the collective rhythm, and we provide a numerical description of the bifurcations resulting from the noise-induced transitions. Our numerical simulations show a noise-induced rhythm generation at low noise intensities, while the SCN clock is arrhythmic in the high noise setting. Notably, coupling induces resonance-like behavior at low noise intensities, and varying coupling strength can cause period locking and variance dissipation even in the presence of noise.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A simplified voltage-conductance kinetic model for interacting neurons and its asymptotic limit

The voltage-conductance kinetic model for the collective behavior of neurons has been studied by scientists and mathematicians for two decades, but the rigorous analysis of its solution structure has been only partially obtained in spite of plenty of numerical evidence in various scenarios. In this work, we consider a simplified voltage-conductance model in which the velocity field in the voltage variable is in a separable form. The long time behavior of the simplified model is fully investigated leading to the following dichotomy: either the density function converges to the global equilibrium, or the firing rate diverges as time goes to infinity. Besides, the fast conductance asymptotic limit is justified and analyzed, where the solution to the limit model either blows up in finite time, or globally exists leading to time periodic solutions. An important implication of these results is that the non-separable velocity field, or physically the leaky mechanism, is a key element for the emergence of periodic solutions in the original model based on the available numerical evidence.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

From radial symmetry to fractal behavior of aggregation equilibria for repulsive-attractive potentials

For the interaction energy with repulsive-attractive potentials, we give generic conditions which guarantee the radial symmetry of the local minimizers in the infinite Wasserstein distance. As a consequence, we obtain the uniqueness of local minimizers in this topology for a class of interaction potentials. We introduce a novel notion of concavity of the interaction potential allowing us to show certain fractal-like behavior of the local minimizers. We provide a family of interaction potentials such that the support of the associated local minimizers has no isolated points and any superlevel set has no interior points.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Minimizers of 3D anisotropic interaction energies

We study a large family of axisymmetric Riesz-type singular interaction potentials with anisotropy in three dimensions. We generalize some of the results of our recent work in two dimensions to the present setting. For potentials with linear interpolation convexity, their associated global energy minimizers are given by explicit formulas whose supports are ellipsoids. We show that for less singular anisotropic Riesz potentials, the global minimizer may collapse into one or two dimensional concentrated measures which minimize restricted isotropic Riesz interaction energies. Some partial aspects of these questions are also tackled in the intermediate range of singularities in which one dimensional vertical collapse is not allowed. Collapse to lower dimensional structures is proved at the critical value of the convexity but not necessarily to vertically or horizontally concentrated measures, leading to interesting open problems.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A Finite-Volume Method for Fluctuating Dynamical Density Functional Theory

We introduce a finite-volume numerical scheme for solving stochastic gradient-flow equations. Such equations are of crucial importance within the framework of fluctuating hydrodynamics and dynamic density functional theory. Our proposed scheme deals with general free-energy functionals, including, for instance, external fields or interaction potentials. This allows us to simulate a range of physical phenomena where thermal fluctuations play a crucial role, such as nucleation and other energy-barrier crossing transitions. A positivity-preserving algorithm for the density is derived based on a hybrid space discretization of the deterministic and the stochastic terms and different implicit and explicit time integrators. We show through numerous applications that not only our scheme is able to accurately reproduce the statistical properties (structure factor and correlations) of the physical system, but, because of the multiplicative noise, it allows us to simulate energy barrier crossing dynamics, which cannot be captured by mean-field approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A second-order numerical method for the aggregation equations

Inspired by so-called TVD limiter-based second-order schemes for hyperbolic conservation laws, we develop a second-order accurate numerical method for multi-dimensional aggregation equations. The method allows for simulations to be continued after the first blow-up time of the solution. In the case of symmetric, lambda-convex potentials with a possible Lipschitz singularity at the origin we prove that the method converges in the Monge--Kantorovich distance towards the unique gradient flow solution. Several numerical experiments are presented to validate the second-order convergence rate and to explore the performance of the scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A consensus-based global optimization method for high dimensional machine learning problems

We improve recently introduced consensus-based optimization method, proposed in [R. Pinnau, C. Totzeck, O. Tse and S. Martin, Math. Models Methods Appl. Sci., 27(01):183--204, 2017], which is a gradient-free optimization method for general non-convex functions. We first replace the isotropic geometric Brownian motion by the component-wise one, thus removing the dimensionality dependence of the drift rate, making the method more competitive for high dimensional optimization problems. Secondly, we utilize the random mini-batch ideas to reduce the computational cost of calculating the weighted average which the individual particles tend to relax toward. For its mean-field limit--a nonlinear Fokker-Planck equation--we prove, in both time continuous and semi-discrete settings, that the convergence of the method, which is exponential in time, is guaranteed with parameter constraints {\it independent} of the dimensionality. We also conduct numerical tests to high dimensional problems to check the success rate of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An Optimal Mass Transport Method for Random Genetic Drift

We propose and analyze an optimal mass transport method for a random genetic drift problem driven by a Moran process under weak-selection. The continuum limit, formulated as a reaction-advection-diffusion equation known as the Kimura equation, inherits degenerate diffusion from the discrete stochastic process that conveys to the blow-up into Dirac-delta singularities hence brings great challenges to both the analytical and numerical studies. The proposed numerical method can quantitatively capture to the fullest possible extent the development of Dirac-delta singularities for genetic segregation on one hand, and preserves several sets of biologically relevant and computationally favored properties of the random genetic drift on the other. Moreover, the numerical scheme exponentially converges to the unique numerical stationary state in time at a rate independent of the mesh size up to a mesh error. Numerical evidence is given to illustrate and support these properties, and to demonstrate the spatio-temporal dynamics of random generic drift.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Convergence of a Finite Volume Scheme for a System of Interacting Species with Cross-Diffusion

In this work we present the convergence of a positivity preserving semi-discrete finite volume scheme for a coupled system of two non-local partial differential equations with cross-diffusion. The key to proving the convergence result is to establish positivity in order to obtain a discrete energy estimate to obtain compactness. We numerically observe the convergence to reference solutions with a first order accuracy in space. Moreover we recover segregated stationary states in spite of the regularising effect of the self-diffusion. However, if the self-diffusion or the cross-diffusion is strong enough, mixing occurs while both densities remain continuous.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Numerical study of Bose-Einstein condensation in the Kaniadakis-Quarati model for bosons

Kaniadakis and Quarati (1994) proposed a Fokker--Planck equation with quadratic drift as a PDE model for the dynamics of bosons in the spatially homogeneous setting. It is an open question whether this equation has solutions exhibiting condensates in finite time. The main analytical challenge lies in the continuation of exploding solutions beyond their first blow-up time while having a linear diffusion term. We present a thoroughly validated time-implicit numerical scheme capable of simulating solutions for arbitrarily large time, and thus enabling a numerical study of the condensation process in the Kaniadakis--Quarati model. We show strong numerical evidence that above the critical mass rotationally symmetric solutions of the Kaniadakis--Quarati model in 3D form a condensate in finite time and converge in entropy to the unique minimiser of the natural entropy functional at an exponential rate. Our simulations further indicate that the spatial blow-up profile near the origin follows a universal power law and that transient condensates can occur for sufficiently concentrated initial data.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Phase transitions for nonlinear nonlocal aggregation-diffusion equations

We are interested in studying the stationary solutions and phase transitions of aggregation equations with degenerate diffusion of porous medium-type, with exponent $1 < m < \infty$. We first prove the existence of possibly infinitely many bifurcations from the spatially homogeneous steady state. We then focus our attention on the associated free energy proving existence of minimisers and even uniqueness for sufficiently weak interactions. In the absence of uniqueness, we show that the system exhibits phase transitions: we classify values of $m$ and interaction potentials $W$ for which these phase transitions are continuous or discontinuous. Finally, we comment on the limit $m \to \infty$ and the influence that the presence of a phase transition has on this limit.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantifying the hydrodynamic limit of Vlasov-type equations with alignment and nonlocal forces

In this paper, we quantify the asymptotic limit of collective behavior kinetic equations arising in mathematical biology modeled by Vlasov-type equations with nonlocal interaction forces and alignment. More precisely, we investigate the hydrodynamic limit of a kinetic Cucker--Smale flocking model with confinement, nonlocal interaction, and local alignment forces, linear damping and diffusion in velocity. We first discuss the hydrodynamic limit of our main equation under strong local alignment and diffusion regime, and we rigorously derive the isothermal Euler equations with nonlocal forces. We also analyze the hydrodynamic limit corresponding to strong local alignment without diffusion. In this case, the limiting system is pressureless Euler-type equations. Our analysis includes the Coulombian interaction potential for both cases and explicit estimates on the distance towards the limiting hydrodynamic equations. The relative entropy method is the crucial technology in our main results, however, for the case without diffusion, we combine a modulated macroscopic kinetic energy with the bounded Lipschitz distance to deal with the nonlocality in the interaction forces. For the sake of completeness, the existence of weak and strong solutions to the kinetic and fluid equations are also established.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantitative error estimates for the large friction limit of Vlasov equation with nonlocal forces

We study an asymptotic limit of Vlasov type equation with nonlocal interaction forces where the friction terms are dominant. We provide a quantitative estimate of this large friction limit from the kinetic equation to a continuity type equation with a nonlocal velocity field, the so-called aggregation equation, by employing $2$-Wasserstein distance. By introducing an intermediate system, given by the pressureless Euler equations with nonlocal forces, we can quantify the error between the spatial densities of the kinetic equation and the pressureless Euler system by means of relative entropy type arguments combined with the $2$-Wasserstein distance. This together with the quantitative error estimate between the pressureless Euler system and the aggregation equation in $2$-Wasserstein distance in [Commun. Math. Phys, 365, (2019), 329--361] establishes the quantitative bounds on the error between the kinetic equation and the aggregation equation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Traveling Bands, Clouds, and Vortices of Chiral Active Matter

We consider stochastic dynamics of self-propelled particles with nonlocal normalized alignment interactions subject to phase lag. The role of the lag is to indirectly generate chirality into particle motion. To understand large scale behavior, we derive a continuum description of an active Brownian particle (ABP) flow with macroscopic scaling in the form of a partial differential equation (PDE) for a one-particle probability density function (DF). Due to indirect chirality, we find a new spatially homogeneous nonstationary analytic solution for this class of equations. Our development of kinetic and hydrodynamic theories towards such a solution reveals the existence of a wide variety of spatially nonhomogeneous patterns reminiscent of the traveling bands, clouds, and vortical structures of linear active matter. Our model may thereby serve as the basis for understanding the nature of chiral active media and designing multiagent swarms with designated behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

On the singularity formation and relaxation to equilibrium in 1D Fokker-Planck model with superlinear drift

We consider a class of Fokker--Planck equations with linear diffusion and superlinear drift enjoying a formal Wasserstein-like gradient flow structure with convex mobility function. In the drift-dominant regime, the equations have a finite critical mass above which the measure minimising the associated entropy functional displays a singular component. Our approach, which addresses the one-dimensional case, is based on a reformulation of the problem in terms of the pseudo-inverse distribution function. Motivated by the structure of the equation in the new variables, we establish a general framework for global-in-time existence, uniqueness and regularity of monotonic viscosity solutions to a class of nonlinear degenerate (resp. singular) parabolic equations, using as a key tool comparison principles and maximum arguments. We then focus on a specific equation and study in more detail the regularity and dynamics of solutions. In particular, blow-up behaviour, formation of condensates (i.e. Dirac measures at zero) and long-time asymptotics are investigated. As a consequence, in the mass-supercritical case, solutions will blow up in $L^\infty$ in finite time and---understood in a generalised, measure sense---they will eventually have condensate. We further show that the singular part of the measure solution does in general interact with the density and that condensates can be transient. The equations considered are motivated by a model for bosons introduced by Kaniadakis and Quarati (1994), which has a similar entropy structure and a critical mass if $d\ge3$.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Convergence to Equilibrium in Wasserstein distance for damped Euler equations with interaction forces

We develop tools to construct Lyapunov functionals on the space of probability measures in order to investigate the convergence to global equilibrium of a damped Euler system under the influence of external and interaction potential forces with respect to the 2-Wasserstein distance. We also discuss the overdamped limit to a nonlocal equation used in the modelling of granular media with respect to the 2-Wasserstein distance, and provide rigorous proofs for particular examples in one spatial dimension.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Compactly supported stationary states of the degenerate Keller-Segel system in the diffusion-dominated regime

We first show the existence of unique global minimizer of the free energy for all masses associated to a nonlinear diffusion version of the classical Keller-Segel model when the diffusion dominates over the attractive force of the chemoattractant. The strategy uses an approximation of the variational problem in the whole space by the minimization problem posed on bounded balls with large radii. We show that all stationary states in a wide class coincide up to translations with the unique compactly supported radially decreasing and smooth inside its support global minimizer of the free energy. Our results complement and show alternative proofs with respect to \cite{Strohmer,LY,Kim-Yao}.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Displacement convexity for the entropy in semidiscrete nonlinear Fokker-Planck equations

The displacement $λ$-convexity of a nonstandard entropy with respect to a nonlocal transportation metric in finite state spaces is shown using a gradient flow approach. The constant $λ$ is computed explicitly in terms of a priori estimates of the solution to a finite-difference approximation of a nonlinear Fokker-Planck equation. The key idea is to employ a new mean function, which defines the Onsager operator in the gradient flow formulation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Explicit Equilibrium Solutions For the Aggregation Equation with Power-Law Potentials

Despite their wide presence in various models in the study of collective behaviors, explicit swarming patterns are difficult to obtain. In this paper, special stationary solutions of the aggregation equation with power-law kernels are constructed by inverting Fredholm integral operators or by employing certain integral identities. These solutions are expected to be the global energy stable equilibria and to characterize the generic behaviors of stationary solutions for more general interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Mean-field limit for collective behavior models with sharp sensitivity regions

We rigorously show the mean-field limit for a large class of swarming individual based models with local sharp sensitivity regions. For instance, these models include nonlocal repulsive-attractive forces locally averaged over sharp vision cones and Cucker-Smale interactions with discontinuous communication weights. We construct global-in-time defined notion of solutions through a differential inclusion system corresponding to the particle descriptions. We estimate the error between the solutions to the differential inclusion system and weak solutions to the expected limiting kinetic equation by employing tools from optimal transport theory. Quantitative bounds on the expansion of the 1-Wasserstein distance along flows based on a weak-strong stability estimate are obtained. We also provide different examples of realistic sensitivity sets satisfying the assumptions of our main results.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Numerical simulation of nonlinear continuity equations by evolving diffeomorphisms

In this paper we present a numerical scheme for nonlinear continuity equations, which is based on the gradient flow formulation of an energy functional with respect to the quadratic transportation distance. It can be applied to a large class of nonlinear continuity equations, whose dynamics are driven by internal energies, given external potentials and/or interaction energies. The solver is based on its variational formulation as a gradient flow with respect to the Wasserstein distance. Positivity of solutions as well as energy decrease of the semi-discrete scheme are guaranteed by its construction. We illustrate this properties with various examples in spatial dimension one and two.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the pressureless damped Euler-Poisson equations with non-local forces: Critical thresholds and large-time behavior

We analyse the one-dimensional pressureless Euler-Poisson equations with a linear damping and non-local interaction forces. These equations are relevant for modelling collective behavior in mathematical biology. We provide a sharp threshold between the supercritical region with finite-time breakdown and the subcritical region with global-in-time existence of the classical solution. We derive an explicit form of solution in Lagrangian coordinates which enables us to study the time-asymptotic behavior of classical solutions with the initial data in the subcritical region.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A local version of the Hughes model for pedestrian flow

Roger Hughes proposed a macroscopic model for pedestrian dynamics, in which individuals seek to minimize their travel time but try to avoid regions of high density. One of the basic assumptions is that the overall density of the crowd is known to every agent. In this paper we present a modification of the Hughes model to include local effects, namely limited vision, and a conviction towards decision making. The modified velocity field enables smooth turning and temporary waiting behavior. We discuss the modeling in the micro- and macroscopic setting as well as the efficient numerical simulation of either description. Finally we illustrate the model with various numerical experiments and evaluate the behavior with respect to the evacuation time and the overall performance.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Convergence of a linearly transformed particle method for aggregation equations

We study a linearly transformed particle method for the aggregation equation with smooth or singular interaction forces. For the smooth interaction forces, we provide convergence estimates in $L^1$ and $L^\infty$ norms depending on the regularity of the initial data. Moreover, we give convergence estimates in bounded Lipschitz distance for measure valued solutions. For singular interaction forces, we establish the convergence of the error between the approximated and exact flows up to the existence time of the solutions in $L^1 \cap L^p$ norm.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Critical thresholds in 1D Euler equations with nonlocal forces

We study the critical thresholds for the compressible pressureless Euler equations with pairwise attractive or repulsive interaction forces and non-local alignment forces in velocity in one dimension. We provide a complete description for the critical threshold to the system without interaction forces leading to a sharp dichotomy condition between global in time existence or finite-time blow-up of strong solutions. When the interaction forces are considered, we also give a classification of the critical thresholds according to the different type of interaction forces. We also analyze conditions for global in time existence when the repulsion is modeled by the isothermal pressure law.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Measure solutions for some models in population dynamics

We give a direct proof of well-posedness of solutions to general selection-mutation and structured population models with measures as initial data. This is motivated by the fact that some stationary states of these models are measures and not $L^1$ functions, so the measures are a more natural space to study their dynamics. Our techniques are based on distances between measures appearing in optimal transport and common arguments involving Picard iterations. These tools provide a simplification of previous approaches and are applicable or adaptable to a wide variety of models in population dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Weak solutions for Euler systems with non-local interactions

We consider several modifications of the Euler system of fluid dynamics including its pressureless variant driven by non-local interaction repulsive-attractive and alignment forces in the space dimension $N=2,3$. These models arise in the study of self-organisation in collective behavior modeling of animals and crowds. We adapt the method of convex integration to show the existence of infinitely many global-in-time weak solutions for any bounded initial data. Then we consider the class of \emph{dissipative} solutions satisfying, in addition, the associated global energy balance (inequality). We identify a large set of initial data for which the problem admits infinitely many dissipative weak solutions. Finally, we establish a weak-strong uniqueness principle for the pressure driven Euler system with non-local interaction terms as well as for the pressureless system with Newtonian interaction.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Finite-Volume Method for Nonlinear Nonlocal Equations with a Gradient Flow Structure

We propose a positivity preserving entropy decreasing finite volume scheme for nonlinear nonlocal equations with a gradient flow structure. These properties allow for accurate computations of stationary states and long-time asymptotics demonstrated by suitably chosen test cases in which these features of the scheme are essential. The proposed scheme is able to cope with non-smooth stationary states, different time scales including metastability, as well as concentrations and self-similar behavior induced by singular nonlocal kernels. We use the scheme to explore properties of these equations beyond their present theoretical knowledge.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A hybrid variational principle for the Keller-Segel system in $\mathbb R^2$

We construct weak global in time solutions to the classical Keller-Segel system cell movement by chemotaxis in two dimensions when the total mass is below the well-known critical value. Our construction takes advantage of the fact that the Keller-Segel system can be realized as a gradient flow in a suitable functional product space. This allows us to employ a hybrid variational principle which is a generalisation of the minimising implicit scheme for Wasserstein distances introduced by Jordan, Kinderlehrer and Otto (1998).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Non-local kinetic and macroscopic models for self-organised animal aggregations

The last two decades have seen a surge in kinetic and macroscopic models derived to investigate the multi-scale aspects of self-organised biological aggregations. Because the individual-level details incorporated into the kinetic models (e.g., individual speeds and turning rates) make them somewhat difficult to investigate, one is interested in transforming these models into simpler macroscopic models, by using various scaling techniques that are imposed by the biological assumptions of the models. Here, we consider three scaling approaches (parabolic, hydrodynamic and grazing collision limits) that can be used to reduce a class of non-local 1D and 2D models for biological aggregations to simpler models existent in the literature. Next, we investigate how some of the spatio-temporal patterns exhibited by the original kinetic models are preserved via these scalings. To this end, we focus on the parabolic scaling for non-local 1D models and apply asymptotic preserving numerical methods, which allow us to analyse changes in the patterns as the scaling coefficient $ε$ is varied from $ε=1$ (for 1D transport models) to $ε=0$ (for 1D parabolic models). We show that some patterns (describing stationary aggregations) are preserved in the limit $ε\to 0$, while other patterns (describing moving aggregations) are lost in this limit. To understand the loss of these patterns, we construct bifurcation diagrams.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nonlinear diffusion: Geodesic Convexity is equivalent to Wasserstein Contraction

It is well known that nonlinear diffusion equations can be interpreted as a gradient flow in the space of probability measures equipped with the Euclidean Wasserstein distance. Under suitable convexity conditions on the nonlinearity, due to R. J. McCann, the associated entropy is geodesically convex, which implies a contraction type property between all solutions with respect to this distance. In this note, we give a simple straightforward proof of the equivalence between this contraction type property and this convexity condition, without even resorting to the entropy and the gradient flow structure.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nonlocal-interaction equations on uniformly prox-regular sets

We study the well-posedness of a class of nonlocal-interaction equations on general domains $Ω\subset \mathbb{R}^d$, including nonconvex ones. We show that under mild assumptions on the regularity of domains (uniform prox-regularity), for $λ$-geodesically convex interaction and external potentials, the nonlocal-interaction equations have unique weak measure solutions. Moreover, we show quantitative estimates on the stability of solutions which quantify the interplay of the geometry of the domain and the convexity of the energy. We use these results to investigate on which domains and for which potentials the solutions aggregate to a single point as time goes to infinity. Our approach is based on the theory of gradient flows in spaces of probability measures.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A numerical study of the Navier-Stokes transport coefficients for 2D granular hydrodynamics

A numerical study is presented to analyze the thermal mechanisms of unsteady, supersonic granular flow, by means of hydrodynamic simulations of the Navier-Stokes granular equations. For this purpose a paradigmatic problem in granular dynamics such as the Faraday instability is selected. Two different approaches for the Navier-Stokes transport coefficients for granular materials are considered, namely the traditional Jenkins-Richman theory for moderately dense quasi-elastic grains, and the improved Garzó-Dufty-Lutsko theory for arbitrary inelasticity, which we also present here. Both solutions are compared with event-driven simulations of the same system under the same conditions, by analyzing the density, the temperature and the velocity field. Important differences are found between the two approaches leading to interesting implications. In particular, the heat transfer mechanism coupled to the density gradient which is a distinctive feature of inelastic granular gases, is responsible for a major discrepancy in the temperature field and hence in the diffusion mechanisms.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Complexity Reduction of Rate-Equations Models for Two-Choice Decision-Making

We are concerned with the complexity reduction of a stochastic system of differential equations governing the dynamics of a neuronal circuit describing a decision-making task. This reduction is based on the slow-fast behavior of the problem and holds on the whole phase space and not only locally around the spontaneous state. Macroscopic quantities, such as performance and reaction times, computed applying this reduction are in agreement with previous works in which the complexity reduction is locally performed at the spontaneous point by means of a Taylor expansion.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Equivalence of gradient flows and entropy solutions for singular nonlocal interaction equations in 1D

We prove the equivalence between the notion of Wasserstein gradient flow for a one-dimensional nonlocal transport PDE with attractive/repulsive Newtonian potential on one side, and the notion of entropy solution of a Burgers-type scalar conservation law on the other. The solution of the former is obtained by spatially differentiating the solution of the latter. The proof uses an intermediate step, namely the $L^2$ gradient flow of the pseudo-inverse distribution function of the gradient flow solution. We use this equivalence to provide a rigorous particle-system approximation to the Wasserstein gradient flow, avoiding the regularization effect due to the singularity in the repulsive kernel. The abstract particle method relies on the so-called wave-front-tracking algorithm for scalar conservation laws. Finally, we provide a characterization of the sub-differential of the functional involved in the Wasserstein gradient flow.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Classical Solutions for a nonlinear Fokker-Planck equation arising in Computational Neuroscience

In this paper we analyze the global existence of classical solutions to the initial boundary-value problem for a nonlinear parabolic equation describing the collective behavior of an ensemble of neurons. These equations were obtained as a diffusive approximation of the mean-field limit of a stochastic differential equation system. The resulting Fokker-Planck equation presents a nonlinearity in the coefficients depending on the probability flux through the boundary. We show by an appropriate change of variables that this parabolic equation with nonlinear boundary conditions can be transformed into a non standard Stefan-like free boundary problem with a source term given by a delta function. We prove that there are global classical solutions for inhibitory neural networks, while for excitatory networks we give local well-posedness of classical solutions together with a blow up criterium. Finally, we will also study the spectrum for the linear problem corresponding to uncoupled networks and its relation to Poincaré inequalities for studying their asymptotic behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Decay rates for a class of diffusive-dominated interaction equations

We analyse qualitative properties of the solutions to a mean-field equation for particles interacting through a pairwise potential while diffusing by Brownian motion. Interaction and diffusion compete with each other depending on the character of the potential. We provide sufficient conditions on the relation between the interaction potential and the initial data for diffusion to be the dominant term. We give decay rates of Sobolev norms showing that asymptotically for large times the behavior is then given by the heat equation. Moreover, we show an optimal rate of convergence in the $L^1$-norm towards the fundamental solution of the heat equation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Mean-field limit for the stochastic Vicsek model

We consider the continuous version of the Vicsek model with noise, proposed as a model for collective behavior of individuals with a fixed speed. We rigorously derive the kinetic mean-field partial differential equation satisfied when the number N of particles tends to infinity, quantifying the convergence of the law of one particle to the solution of the PDE. For this we adapt a classical coupling argument to the present case in which both the particle system and the PDE are defined on a surface rather than on the whole space. As part of the study we give existence and uniqueness results for both the particle system and the PDE.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A well-posedness theory in measures for some kinetic models of collective motion

We present existence, uniqueness and continuous dependence results for some kinetic equations motivated by models for the collective behavior of large groups of individuals. Models of this kind have been recently proposed to study the behavior of large groups of animals, such as flocks of birds, swarms, or schools of fish. Our aim is to give a well-posedness theory for general models which possibly include a variety of effects: an interaction through a potential, such as a short-range repulsion and long-range attraction; a velocity-averaging effect where individuals try to adapt their own velocity to that of other individuals in their surroundings; and self-propulsion effects, which take into account effects on one individual that are independent of the others. We develop our theory in a space of measures, using mass transportation distances. As consequences of our theory we show also the convergence of particle systems to their corresponding kinetic equations, and the local-in-time convergence to the hydrodynamic limit for one of the models.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Analysis of Nonlinear Noisy Integrate\&Fire Neuron Models: blow-up and steady states

Nonlinear Noisy Leaky Integrate and Fire (NNLIF) models for neurons networks can be written as Fokker-Planck-Kolmogorov equations on the probability density of neurons, the main parameters in the model being the connectivity of the network and the noise. We analyse several aspects of the NNLIF model: the number of steady states, a priori estimates, blow-up issues and convergence toward equilibrium in the linear case. In particular, for excitatory networks, blow-up always occurs for initial data concentrated close to the firing potential. These results show how critical is the balance between noise and excitatory/inhibitory interactions to the connectivity parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Global Solutions for the One-Dimensional Vlasov-Maxwell System for Laser-Plasma Interaction

We analyse a reduced 1D Vlasov--Maxwell system introduced recently in the physical literature for studying laser-plasma interaction. This system can be seen as a standard Vlasov equation in which the field is split in two terms: an electrostatic field obtained from Poisson's equation and a vector potential term satisfying a nonlinear wave equation. Both nonlinearities in the Poisson and wave equations are due to the coupling with the Vlasov equation through the charge density. We show global existence of weak solutions in the non-relativistic case, and global existence of characteristic solutions in the quasi-relativistic case. Moreover, these solutions are uniquely characterised as fixed points of a certain operator. We also find a global energy functional for the system allowing us to obtain $L^p$-nonlinear stability of some particular equilibria in the periodic setting.

0 citations0 reviewsNo code linkedOpen access