BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

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

Network

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

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

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

Researcher profile

José A. Carrillo

José A. Carrillo contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.APmath.NANumerical Analysiscond-mat.stat-mechNeurons and Cognitionphysics.comp-phcond-mat.softmath.DSmath.OCmath.PRnlin.AOphysics.flu-dyn

Trust snapshot

Quick read

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

17 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&#39;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