Hamilton-Jacobi Equations and Mean Field Games: From Modeling to Numerics via Analysis

Oct 1, 2026, 8:00 AM → Oct 2, 2026, 6:00 PM Europe/Rome

Aula Dini (Palazzo del Castelletto)

Fabio Camilli (Università "G. d'Annunzio" Chieti-Pescara), Elisabetta Carlini (Sapienza Università di Roma), Annalisa Cesaroni (Università degli Studi di Padova), Claudio Marchi (Università degli Studi di Padova)

The conference focuses on the theory, analysis, and numerical treatment of Hamilton-Jacobi (HJ) equations and Mean Field Games (MFG), exploring the interplay between modeling, rigorous analysis, and computational methods. The event aims to bring together experts from analysis, control theory, applied mathematics, and computational sciences to discuss recent advances in the field. More precisely we propose to investigate, through conferences, lectures and discussions, recent advances in:

Analytical approaches: including turnpike properties, homogenization and singular perturbations problems for MFG systems and control problems in infinite dimension, a-priori estimates and regularity properties for solutions to elliptic and parabolic systems possibly degenerate;

Modelling aspects: including Extended MFG for optimal production of exhaustible resources, urban planning problems, micro-macro derivation of fluidynamic models for pedestrian and vehicular traffic, multi-agents systems and flocking;

Numerical methods: including Finite Differences, Semi-Lagrangian, Policy Iterations, Newton Iterations, Machine Learning based Methods for efficiently approximating and solving the large scale HJ equations and MFG systems.

The conference will permit to emphasize the connections between rigorous analysis, modelling and numerical analysis, highlighting how deep theoretical insights can inform the design of efficient and accurate numerical methods.

LIST OF  SPEAKERS

Martino Bardi (Università degli Studi di Padova)

Alessio Basti (Università G. D’Annunzio di Chieti-Pescara)

Jules Berry (Université Paris-Saclay)

Ariela Briani (Institut Denis Poisson, Université de Tours)

Simone Cacace (Sapienza Università di Roma)

Elisa Continelli (Università degli Studi di Padova)

Valentina Coscetti (Sapienza Università di Roma)

Adriano Festa (Politecnico di Torino)

Genni Fragnelli (Università di Siena)

Daria Ghilli (Università di Pavia)

Alessandro Goffi (Università degli Studi di Firenze)

Diogo A. Gomes (KAUST)

Olivier Ley (INSA de Rennes)

Paola Mannucci (Università degli Studi di Padova)

Cristian Mendico (Université Bourgogne Europe)

Antonio Siconolfi (Sapienza Università di Roma)

Francisco J. Silva (Université de Limoges)

Daniela Tonon (Università degli Studi di Padova)

Hasnaa Zidani (INSA Rouen Normandie)

 

The deadline for registration is the 24th of August 2026.

More information on the application procedure at the REGISTRATION PAGE.

The number of participants attending the Workshop is limited, therefore the registration will be confirmed after its expiration.

 

Funded by:

-INdAM-GNAMPA, Gruppo Nazionale per l'Analisi Matematica, la Probabilità e le loro applicazioni

-Dipartimento di Matematica Dipartimento di Eccellenza 2023-2027, Sapienza Università di Roma 

-Dipartimento di Matematica "Tullio Levi-Civita", Università degli Studi di  Padova. 

 

 

Thu, October 1

8:50 AM → 9:30 AM Registration

9:30 AM → 10:00 AM TBA

Speaker: Antonio Siconolfi (Università degli Studi di Roma "La Sapienza")

10:00 AM → 10:30 AM On Hamilton Jacobi equations with time measurable Hamiltonians posed on a 1-dimensional junction

I will describe my recent contribution to the theory of Hamilton-Jacobi equations and optimal control problems posed on networks (see [1]).
I will study evolutive Hamilton–Jacobi equations with Hamiltonians that are discontinuous in time, posed on a simple network consisting of two edges on the real line connected at a single junction. I will introduce a notion of (flux-limited) viscosity solution for Hamiltonians H_i=H_i(t,x,p), (i=1,2) that are assumed to be only measurable in t. Moreover the flux limiter,  A=A(t), acting at the junction, is not required to be continuous but only in L-infinity.
My first motivation is the optimal control problem studied by P. Cardaliaguet et P. Souganidis in [2], where they investigate how to optimally control traffic flow through a junction by acting only on the junction (i.e. with speed reduction or a traffic light) using a merely measurable flux limiter A(t).
In the case of convex Hamiltonians, I prove a comparison principle and establish an existence result via the construction of an optimal control problem. Generalizations to the nonconvex case and to more general networks are also discussed.

[1] A. Briani, On Hamilton Jacobi equations with time measurable Hamiltonians posed on a 1-dimensional junction , Preprint (2026), ArXiv https://arxiv.org/abs/2603.04183.
[2] P. Cardaliaguet, P. E. Souganidis, An optimal control problem of traffic flow on a junction, preprint, (2023), https://arxiv.org/abs/2312.15418.

Speaker: Ariela Briani (Institut Denis Poisson, Université de Tours)

10:30 AM → 11:00 AM Quantitative vanishing viscosity approximation of fully nonlinear, non-convex, Hamilton-Jacobi equations with Hölder data

I will discuss new quantitative estimates of the vanishing viscosity process for evolutionary Hamilton-Jacobi PDEs that are neither concave nor convex in the gradient and Hessian entries. I will describe a novel approach that exploits the regularizing properties of sup/inf-convolutions for viscosity solutions combined with the comparison principle. This method provides explicit sharp constants without assuming differentiability properties neither on solutions nor on the Hamiltonian. This is a joint work with Alekos Cecchin (Padova).

Speaker: Alessandro Goffi (Università di Firenze)

11:00 AM → 11:30 AM Coffee break

11:30 AM → 12:00 PM TBA

Speaker: Hasnaa Zidani (INSA Rouen Normandie)

12:00 PM → 12:30 PM From Discrete Weak KAM Theory to Computable Aubry–Mather Sets

Aubry and Mather sets describe the invariant structures selected by action minimization in Hamiltonian dynamics. Their computation is particularly relevant in celestial mechanics, where resonances and transitions between regular and chaotic motion create intricate dynamical patterns that are difficult to detect using trajectories alone.
In this talk, I will present joint work with Fabio Camilli on semi-discrete and fully discrete approximations of these sets. Our approach combines discrete Lax–Oleinik operators with a finite-dimensional variational formulation for minimizing measures. The resulting schemes approximate critical values, weak KAM solutions, and the supports of invariant minimizing measures, while preserving the variational structure of the continuous problem.
I will also discuss the numerical challenges created by multiple minimizing components and outline how calibration defects and structure-preserving learning may help identify the variational skeleton of resonant Hamiltonian systems.

Speaker: Cristian Mendico (Université Bourgogne Europe)

12:30 PM → 1:00 PM MFG with large discount and agent-based models

We focus on a class of Mean Field Games with discount. If one assumes that the discount factor is large, these models turn out to behave similarly to agent-based models, where individuals just react to the population distribution according to a given rule. Inspired by a recent work by Bardi and Cardaliaguet, we will elaborate on this connection, showing in particular the uniqueness of equilibria in MFG with large discount, and their converge to certain agent-based models, in a quantitative sense, as the discount factor goes to infinity.
Joint work with Marco Cirant.

Speaker: Elisa Continelli (Università degli Studi di Padova)

1:00 PM → 2:30 PM Lunch

2:30 PM → 3:00 PM Hamilton-Jacobi equations on Wasserstein spaces

In this talk, we present a unified approach to Hamilton–Jacobi equations on the Wasserstein space of probability measures arising from deterministic and stochastic dynamics. We introduce a viscosity-solution framework encompassing both first-order equations and semilinear equations driven by idiosyncratic noise. The framework relies on a suitable notion of subdifferential designed to ensure comparison and stability. We then discuss the vanishing-viscosity limit for semilinear Hamilton–Jacobi equations, proving convergence to the corresponding first-order equation as the noise intensity vanishes, with an optimal convergence rate. These results provide a PDE-level characterization of the zero-noise transition.

Speaker: Daniela Tonon (Università degli Studi di Padova)

3:00 PM → 3:30 PM Building on Lax’s Idea: Quantitative Compactness Estimates for Numerical Solutions of First-Order Nonlinear PDEs

In my talk, I will discuss recent results on Kolmogorov ε-entropy for numerical approximations of scalar conservation laws and Hamilton-Jacobi equations. The analysis provides quantitative compactness estimates and entropy-transfer principles linking exact and discrete solution sets, yielding an information-theoretic perspective on complexity preservation and the resolution of numerical schemes in the sense introduced by P. D. Lax.

Speaker: Alessio Basti (Università G. D'Annunzio di Chieti-Pescara)

3:30 PM → 4:00 PM Tea Break

4:00 PM → 4:30 PM A semi-Lagrangian scheme for First-Order Mean Field Games based on monotone operators

We construct a semi-Lagrangian scheme for first-order, time-dependent, and non-local Mean Field Games. The convergence of the scheme to a weak solution of the system is analyzed by exploiting a key monotonicity property. To solve the resulting discrete problem, we implement a Learning Value Algorithm, prove its convergence, and propose an acceleration strategy based on a Policy iteration method. Finally, we present numerical experiments that validate the effectiveness of the proposed schemes and show that the accelerated version significantly improves performance.

Speaker: Valentina Coscetti (Università degli Studi di Roma "La Sapienza")

4:30 PM → 5:00 PM SLTK: A Semi-Lagrangian ToolKit for PDEs on Unstructured Grids

In this talk, I will present recent advancements in the development of SLTK, a custom research library designed for the numerical solution of PDEs on complex domains. By employing semi-Lagrangian schemes on unstructured grids, SLTK adapts various techniques from computational geometry and computer graphics to the context of partial differential equations.
Specifically, I will introduce a generalized point location method for the efficient tracking of characteristics in advection and advection-diffusion problems, including R-tree spatial indexing for the treatment of boundary conditions.
Furthermore, I will present a mass-conservative method for the tracking and remapping of mesh elements deformed by a transport field in continuity and Fokker-Planck equations, also introducing an exact geometric folding algorithm for managing flux-type boundary conditions.
Finally, I will present some numerical results obtained by applying SLTK to optimal control problems, Mean Field Games, and fluid dynamics on possibly non-convex domains with arbitrary holes.
This is a joint work with R. Ferretti (Roma Tre University) and G. Tatafiore (Sapienza University of Rome).

Speaker: Simone Cacace (Università degli Studi di Roma "La Sapienza")

Fri, October 2

9:00 AM → 9:30 AM Monotonicity Methods for time-dependent Mean Field Games:

Abstract. We prove the existence of solutions to first-order, local, time-dependent mean-field games with periodic boundary conditions on the d-dimensional torus. Our approach extends the Banach space monotone operator framework, previously developed for the stationary case, to the substantially harder time-dependent setting. We reformulate the MFG system as a variational inequality for a monotone operator and introduce a low-order p-Laplacian regularization that restores coercivity without high-order smoothing. For the regularized problems, existence follows from an abstract monotone operator theorem. We then derive uniform a priori estimates, including energy bounds, higher integrability via the nonlinear adjoint method, and, by passing to the limit using Minty's method, obtain variational inequality solutions. We further show that these solutions are strong solutions to the original MFG system, adapted to the BV framework. Compared to earlier Hilbert space approaches that rely on high-order elliptic regularization, our method operates in natural Banach spaces, yields stronger solutions, and provides a framework better suited to numerical algorithms.

Speaker: Diogo A. Gomes (KAUST)

9:30 AM → 10:00 AM Mean Field Games with state constraints and Grushin type dynamics

We consider a class of finite horizon deterministic mean field games with nonlocal coupling where the agents must follow Grushin type dynamics with state constraints. We require some assumptions on the local interplay between the set of state constraints and the dynamics.

As a first step, we consider the associated optimal control problem and we establish some properties as: the existence of an optimal trajectory for any starting point $(x,t)$, the closed graph property for the multivalued map which associates to each point $(x,t)$ the set of optimal trajectories starting from that point, the continuity of the value function.
Moreover, using a penalization method and the Maximum Principle, we obtain an uniform bound for the optimal controls. This allows us to obtain a comparison principle and that the value function is the unique constrained viscosity solution of the associated HJ equation.

Afterwards, we tackle the mean field games; taking advantage of the aforementioned properties, we prove the existence of a relaxed equilibrium (which describes the evolution of the game in terms of a probability on the set of admissible trajectories) and derive the existence of a mild solution (which is a couple formed by the value function for the generic player and a family of time dependent measures on the state).
Research project in collaboration with Alessandra Cutri' (Roma Tor Vergata), Claudio Marchi (Universita' di Padova), Nicoletta Tchou (Universite' de Rennes)

Speaker: Paola Mannucci (Università degli Studi di Padova)

10:00 AM → 10:30 AM Approximation of deterministic mean field games under state constraints

In this talk, we study the approximation of equilibria for continuous-time and continuous-space first-order (deterministic) mean field games (Mean Field Games, MFGs) by equilibria of finite MFGs, i.e., games in which the number of time steps and the state space are finite. In the first part, we recall the framework and results obtained in Hadikhanloo–S (2019). In the second part, based on a join work with J. Gianatti and F. Mierez, we extend the analysis to the case where constraints are imposed on the agents’ states. We also present numerical simulations of approximate MFG equilibria in both the unconstrained and constrained cases.

Speaker: Francisco J. Silva (Université de Limoges)

10:30 AM → 11:00 AM Rate of convergence for singular perturbations of Hamilton-Jacobi equations in unbounded spaces

We present rate of convergence results for singular perturbations of Hamilton-Jacobi equations in unbounded spaces where the fast operator is linear, uniformly elliptic and has an Ornstein-Uhlenbeck-type drift. Our achievements are the first rate of convergence results without any periodicity assumption. More in detail, in our HJ equation, the slow operator is a fully nonlinear elliptic operator while the source term is assumed only locally Hölder continuous in both fast and slow variables. We obtain several rates of convergence according to the regularity of the source term.

Speaker: Daria Ghilli (Università di Pavia)

11:00 AM → 11:30 AM Coffee break

11:30 AM → 12:00 PM Degenerate second-order PDEs on networks with Kirchhoff conditions

The aim of this talk is to present results concerning nonlinear degenerate elliptic partial differential equations (PDEs) on networks, with Kirchhoff conditions at the vertices. While elliptic PDEs on networks have long been known, Hamilton-Jacobi equations on networks (the case of totally degenerate PDEs) have been the subject of extensive study over the past decade. In recent work conducted with Guy Barles (Tours) and Erwin Topp (Rio), we obtained results for the general case of degenerate second-order elliptic PDEs.

Speaker: Olivier Ley (INSA de Rennes)

12:00 PM → 12:30 PM On the null controllability of a degenerate Fokker–Planck equation with a drift term

The Fokker–Planck equation describes the time evolution of the probability density function of the velocity for a particle under the influence of drag forces and random forces. In particular, this equation has multiple applications in information theory, graph theory, data science, finance, economics... In one spatial dimension the Fokker–Planck equation for the probability density $p(t,x)$ can be rewritten as
[
p_t(t,x) - (a(t,x)p(t.x))_{xx} + (\mu(t,x)p(t,x))_x= f(t,x)
]
where $t\in [0, T]$, $T>$ is fixed and $x \in (0,1)$.

In this talk we assume that $a$ is a function degenerating at $x=0$; the purpose is to study the null-controllability of the solution, namely the possibility to drive the solution $p$ to rest completely at time $T$.

Speaker: Genni Fragnelli (Università di Siena)

12:30 PM → 1:00 PM Approximation of stable solutions of second order mean field games

In this talk, we are interested in the stationary MFG system:
\begin{equation}
\begin{cases}
-\Delta u + H(x,Du) + \lambda u = F[m] & \text{in } \mathbb{T}^d, \
-\Delta m - \mathrm{div}(m H_p(x,Du)) + \lambda m = \lambda m_0 & \text{in } \mathbb{T}^d,
\end{cases}
\end{equation}
where $u$ is the value function and $m$ the distribution of the players. Our goal is to establish the existence of finite element approximations of the solutions to the MFG ssytem and to obtain the associated error estimates.

Our approach consists in reformulating system in the form $F(u,m) = 0$, where $F$ is a nonlinear map defined on a suitably chosen Banach space. In the case where this map is of class $C^1$, we show that the stable solutions, in the sense of Briani-Cardaliaguet, correspond to the regular zeros of $F$, i.e., those for which $dF[u,m]$ is invertible. This property makes it possible to apply the Brezzi-Rappaz-Raviart (BRR) approximation theorem in order to obtain existence and error estimates for the finite element approximations of the stable solutions.

However, the $C^1$ regularity of $F$ requires the Hamiltonian $H$ to be of class $C^2$, which is rather restrictive from the point of view of optimal control. In order to generalize our approach to the more natural class of Hamiltonians with $C^{1,1}$ regularity, we generalize the BRR theorem to the case of maps $F$ that are Lipschitz continuous and satisfy a metric regularity assumption. This last assumption, coming from variational analysis, replaces the invertibility of $dF[\bar{x}]$. This generalization allows us to extend our error estimates for the finite element approximations of the stable solutions to the MFG system to the case of Hamiltonians of class $C^{1,1}$.

Speaker: Jules Berry (Université Paris-Saclay)

1:00 PM → 2:30 PM Lunch

2:30 PM → 3:00 PM Semi-Lagrangian approximation of viscous transport and conservative equations with one sided Lipschitz velocity fields

The aim of this work is to investigate semi-Lagrangian approximation schemes on unstructured grids for viscous transport and conservative equations with measurable coefficients that satisfy a one-sided Lipschitz condition. To establish the convergence of the schemes, we exploit the characterization of the solution for these equations expressed in terms of measurable time-dependent viscosity solution and, respectively, duality solution. We supplement our theoretical analysis with various numerical examples to illustrate the features of the schemes.

Speakers: Adriano Festa (Politecnico di Torino), Fabio Camilli (Università "G. d'Annunzio" Chieti-Pescara), Dr Luciano Marzufero (Libera Università di Bolzano)

3:00 PM → 3:30 PM Mean convexity and the truncated Laplacian

We propose a notion of mean convexity of a function in R^n with respect to all k-dimensional subspaces, k integer between 1 and n. We show its connection with the truncated Laplacian of order k of the function, i.e., the sum of the first k eigenvalues of the Hessian matrix.
We consider the Hessian PDE that prescribes the k-truncated Laplacian, in viscosity sense. We study existence and uniqueness of solutions to the Dirichlet problem under suitable conditions on the data, in particular convexity properties of the domain. Here we exploit some results by Birindelli, Galise and Ishii.
Our main result is about the inverse problem of the single-valuedness of the truncated Laplacian in viscosity sense. Its solution allows us to characterise k-mean convexity of continuous functions in terms of inequalities for the Hessian equation.
Joint work with P. Mannucci

Speaker: Martino Bardi (Università degli studi di Padova)

3:30 PM → 4:00 PM Tea Break