Keywords
Computer Science and Digital Science
 A1.2.4. QoS, performance evaluation
 A6.1.4. Multiscale modeling
 A6.2.3. Probabilistic methods
 A8.1. Discrete mathematics, combinatorics
 A8.2. Optimization
 A8.3. Geometry, Topology
 A8.6. Information theory
 A8.7. Graph theory
 A8.8. Network science
 A8.9. Performance evaluation
 A9.2. Machine learning
 A9.7. AI algorithmics
Other Research Topics and Application Domains
 B4.3. Renewable energy production
 B6.2.2. Radio technology
 B6.3.4. Social Networks
1 Team members, visitors, external collaborators
Research Scientists
 Bartlomiej Blaszczyszyn [Team leader, Inria, Senior Researcher, HDR]
 François Baccelli [Inria, Senior Researcher, HDR]
 Antoine Baker [Inria, Advanced Research Position, from Nov 2021]
 Ana Busic [Inria, Researcher]
 Christine Fricker [Inria, Researcher, HDR]
 Mir Omid Haji Mirsadeghi [Inria, Advanced Research Position, from Nov 2021]
 Ali Khezeli [Inria, Starting Research Position]
 Marc Lelarge [Inria, Senior Researcher, HDR]
 Laurent Massoulié [Inria, Senior Researcher, HDR]
 Sean Meyn [Inria, Chair, from Jun 2021 until Jul 2021]
 Kevin Scaman [Inria, Starting Faculty Position]
PostDoctoral Fellows
 Simon Coste [Inria, until Sep 2021]
 Ke Feng [Inria, from Nov 2021]
 Sanket Sanjay Kalamkar [Inria, until Feb 2021]
PhD Students
 Claire Bizon Monroc [Inria, from Nov 2021]
 Matthieu Blanke [Inria, from Sep 2021]
 Antoine Brochard [Inria, from Sep 2021]
 Arnaud Cadas [Université Paris Sciences et Lettres, until Nov 2021]
 Eric Daoud [Institut Curie, from Sep 2021]
 Michel Davydov [Ecole normale supérieure ParisSaclay]
 Bastien Dubail [École Normale Supérieure de Paris]
 Mathieu Even [Inria, from Sep 2021]
 Roman Gambelin [Inria]
 Luca Ganassali [Inria]
 Cedric Gerbelot [École Normale Supérieure de Paris, from Sep 2021]
 Hadrien Hendrikx [Inria, until Sep 2021]
 Sayeh Khaniha [Inria]
 Thomas Le Corre [École Normale Supérieure de Paris, from Nov 2021]
 Maxime Leiber [Groupe SAFRAN]
 Pierre Popineau [Inria]
 David Robin [Inria, from Oct 2021]
 Bharath Roy Choudhury [Inria]
 Sébastien Samain [École Normale Supérieure de Paris, until Aug 2021]
 Ilia Shilov [Inria]
 Ludovic Stephan [Sorbonne Université, until Oct 2021]
 Guodong Sun [Nokia, 2021, CIFRE]
 Amaury Triboulin [Inria, from Oct 2021]
 Lucas Weber [DGA, from Oct 2021]
Interns and Apprentices
 Julien Ancel [Inria, from Mar 2021 until Aug 2021]
 Waiss Azizian [École Normale Supérieure de Paris, from Oct 2021]
 Matthieu Blanke [Inria, from Apr 2021 until Aug 2021]
 Eva Bouba [IFPEN, from Apr 2021 until Sep 2021]
 Jean Guillaume Brasier [Inria, from May 2021 until Aug 2021]
 Odilon Duranthon [Inria, from Feb 2021 until Jul 2021]
 Mathieu Even [École Normale Supérieure de Paris, from Mar 2021 until Aug 2021]
 Thomas Le Corre [EDF, from Apr 2021 until Nov 2021]
 Felix Lefebvre [Inria, from Apr 2021 until Aug 2021]
 Bianca MarinMoreno [Inria, from Apr 2021 until Aug 2021]
 Can Pouliquen [Inria, from Jun 2021 until Aug 2021]
 Amaury Triboulin [Inria, from Mar 2021 until Jul 2021]
 Victor Vermès [ENPC, from Apr 2021 until Jul 2021]
Administrative Assistants
 Helene Bessin Rousseau [Inria, from Feb 2021]
 Julien Guieu [Inria, ERC NEMO]
 Helene Milome [Inria]
 Scheherazade Rouag [Inria, until Apr 2021]
Visiting Scientist
 Mathieu Even [École Normale Supérieure de Paris, until Feb 2021]
External Collaborators
 Pierre Bremaud [EPFL and ENS (Honorary Professor)]
 Marc Olivier Buob [Bell Labs (Alcatel)]
 Fabien Mathieu [Nokia, HDR]
2 Overall objectives
The general scientific focus of DYOGENE is on the development of network mathematics. The following theories lie within our research interest: dynamical systems, queuing theory, optimization and control, information theory, stochastic processes, random graphs, stochastic geometry.
Our theoretical developments are motivated by and applied in the context of communication networks (Internet, wireless, mobile, cellular, peertopeer), social and economic networks, power grids, and, recently, infectious diseases.
We collaborate with many industrial partners. Our current industrial relations involve EDF, Huawei, Microsoft, Nokia, Orange, Safran.
More specifically, the scientific focus of DYOGENE defined in 2013 was on geometric network dynamics arising in communications. By geometric networks we understand networks with a nontrivial, discrete or continuous, geometric definition of the existence of links between the nodes. In stochastic geometric networks, this definition leads to random graphs or stochastic geometric models.
A first type of geometric network dynamics is the one where the nodes or the links change over time according to an exogeneous dynamics (e.g. node motion and geometric definition of the links). We will refer to this as dynamics of geometric networks below. A second type is that where links and/or nodes are fixed but harbor local dynamical systems (in our case, stemming from e.g. information theory, queuing theory, social and economic sciences). This will be called dynamics on geometric networks. A third type is that where the dynamics of the network geometry and the local dynamics interplay. Our motivations for studying these systems stem from many fields of communications where they play a central role, and in particular: message passing algorithms; epidemic algorithms; wireless networks and information theory; device to device networking; distributed content delivery; social and economic networks, neural networks, and power grids.
3 Research program
3.1 Initial research axes
The following research axes have been defined in 2013 when the projectteam was created.
 Algorithms for network performance analysis, led by A. Bouillard and A. Busic.
 Stochastic geometry and information theory for wireless network, led by F. Baccelli and B. Blaszczyszyn.
 The cavity method for network algorithms, led by M. Lelarge.
Our scientific interests keep evolving. Research areas which received the most of our attention in 2021 are summarized in the following sections.
3.2 Models of infectious diseases
Over the past years, with several researchers and collaborations with we have looked at mathematical models of the evolution of the epidemics, in particular Covid 19. This year we studied the role of geographic mobility on the propagation of epidemics using point process techniques, see 4.
3.3 Distributed network control and smartgrids
Theory and algorithms for distributed control of networks with applications to the stabilization of power grids subject to high volatility of renewable energy production are being developed by A. Busic in collaboration with Sean Meyn [Prof. at University of Florida and Inria International Chair] and Prabir Barooah [University of Florida]. We extended the KullbackLeiblerquadratic optimal control approach to the case of stochastic disturbance. Also new advances have been obtained for thermostatically controlled loads. Within our collaboration with Vito (Belgium) that started in 2019 the main focus in 2021 has been on the Generalized Nash Equilibrium models for P2P markets.
3.4 Reinforcement learning
This year, in a book chapter 29, we reviewed the foundations of reinforcement learning algorithm design based on recent and more classical results from stochastic approximation. A new collaboration with IFPEN started on reinforcement learning with constraints and applications to wind farm production optimization, with PhD thesis of Claire Bizon Monroc.
3.5 Mathematics of wireless cellular networks
This year we revisited our previous line of thought on probabilistic modeling of geographic caching in wireless networks, by contributing a chapter 30 to the book with editors V. Poor [Princeton University] and W. Chen [Tsinghua University], www.bibliovault.org/BV.book.epl?ISBN=9781839531231. In a more fundamental aspect, in 7 we addressed the problem of spacial scaling of general probabilistic policies (e.g. of caching), going beyond Gibbs framework. Also we continued to work on cellular network dimensioning toolbox in a longterm collaboration between TREC/DYOGENE represented by B. Błaszczyszyn, and Orange Labs, represented by M. K. Karray. Furthemore, a collaboration with the Standardization and Research Lab at Nokia Bell Labs and ERC NEMO led by F. Baccelli, started in 2019 and led to several joint publications including 11 New directions started on vehicular networks and particularly on V2X in collaborations with Chang Sik Choi (Hongik University, Korea) and Nithin Ramesan (UT Austin).
3.6 Highdimensional statistical inference and distributed learning
We computed information theoretic bounds for unsupervised and semisupervised learning and proved complexity bounds for distributed optimization of convex functions using a network of computing units.
3.7 Distributed optimization for machine learning
We obtained new variance reduction results on stochastic Bregman gradient methods 21. We introduced the "continuized" Nesterov acceleration, a close variant of Nesterov acceleration whose variables are indexed by a continuous time parameter 22, and obtained the first rigorous acceleration of asynchronous gossip algorithms. We also obtained concentration results of nonisotropic random tensors and studied their pplications to learning and empirical risk minimization 23.
3.8 Random Geometry
Point processes and stochastic geometry offer a mathematical framework for the analysis of various random structures and dynamics embedded in Euclidean spaces. In particular, it has enabled efficient analysis of truly deployed wireless networks. Behind this are mathematical asymptotic results (involving growing number of nodes subject to random processes) often leading to ”Poissonisation” of network architectures, and enabling probabilistic analysis. More generally, the analysis of large, homogeneous networks can rely on the key concept of a typical network element (user, base station, link) and its local relationships representing the overall behavior of the network. Indeed, the framework of the stochastic geometry, equipped with Palm calculus for stationary models, and its analog framework for random unimodular graphs, essentially relate to the properties of a typical network element through the socalled ”mass transport” principle. This paradigm was extended beyond the Euclidean case, based on the theory of unimodular networks, see 2. We developed the probabilistic machinery and also studied applications in a variery of domains: statistical physics, combinatorial optimization, communications, particle gradient descent model for point process generation, processes on Delaunay neighbors in the PoissonVoronoi tessellation, Dirichlet measures, stochastic games, etc. We collaborated on the matter with V. Anantharam (EECS at UC Berkeley), Ch. Hirsch (University of Groningen), S. Mallat (ENS/Flatiron Institute) and S. Zhang (IRITSC).
3.9 Mathematics of stochastic networks via meanfield analysis
In order to understand the fundamental processes taking place in very large networks, it is clear that we sometimes have to abandon the representation of local details and focus on the more macroscopic properties of higher scales. Indeed, the main idea of Mean Field Theory is to replace all detailed bodytobody interactions with a typical element of the network guided by an average or effective interaction. This year, we have several results in this area related to neural networks, opinion dynamics and carsharing systems.
4 Application domains
4.1 Physical communication networks
Internet, wireless, mobile, cellular networks, transportation networks, distributed systems (cloud, call centers). In collaboration with Nokia Bell Labs and Orange Labs.
4.2 Abstract networks
Social interactions, human communities, economic networks.
4.3 Power grids
Energy networks. In collaboration with EDF and Vito (Belgium).
5 Social and environmental responsibility
6 Highlights of the year
6.1 Awards
 The paper 22 on "A Continuized View on Nesterov Acceleration for Stochastic Gradient Descent and Randomized Gossip" by Mathieu Even, Raphaël Berthier, Francis Bach, Nicolas Flammarion, Pierre Gaillard, Hadrien Hendrikx, Laurent Massoulié and Adrien Taylor received a outstanding paper award at Neurips 2021; cf neurips.cc/virtual/2021/awards_detail.
 Bianca Marin Moreno received a research internship award (section Applied mathematics) from Ecole Polytechnique for her internship entitled "A Stochastic Model and Data Analysis of the Incentive Policy of a FreeFloating Car Sharing System".
 Julien Ancel received congratulations (section Applied mathematics) from Ecole Polytechnique for his research internship entitled "Stochastic modelisation of 2classes bikesharing networks".
7 New results
Participants: All Dyogene.
7.1 Models of infectious diseases
1. The role of mobility on the propagation of epidemics on point processes 4 This paper is focused on SIS (SusceptibleInfectedSusceptible) epidemic dynamics (also known as the contact process) on populations modeled by homogeneous Poisson point processes of the Euclidean plane, where the infection rate of a susceptible individual is proportional to the number of infected individuals in a disc around it. The main focus of the paper is a model where points are also subject to some random motion. Conservation equations for moment measures are leveraged to analyze the stationary regime of the point processes of infected and susceptible individuals. A heuristic factorization of the third moment measure is then proposed to obtain simple polynomial equations allowing one to derive closed form approximations for the fraction of infected individuals in the steady state. These polynomial equations also lead to a phase diagram which tentatively delineates the regions of the space of parameters (population density, infection radius, infection and recovery rate, and motion rate) where the epidemic survives and those where there is extinction. A key takeaway from this phase diagram is that the extinction of the epidemic is not always aided by a decrease in the motion rate. These results are substantiated by simulations on large two dimensional tori. These simulations show that the polynomial equations accurately predict the fraction of infected individuals when the epidemic survives. The simulations also show that the proposed phase diagram accurately predicts the parameter regions where the mean survival time of the epidemic increases (resp. decreases) with motion rate.
7.2 Distributed network control and smartgrids
2. KullbackLeiblerQuadratic Optimal Control in a Stochastic Environment 19 This work presents advances in KullbackLeiblerQuadratic (KLQ) optimal control: a stochastic control framework for Markovian models. The motivation is distributed control of large networked systems. The objective function is composed of a control cost in the form of KullbackLeibler divergence plus a quadratic cost on the sequence of marginal distributions. With this choice of objective function, the optimal probability distribution of a population of agents over a finite time horizon is shown to be an exponential tilting of the nominal probability distribution. The same is true for the controlled transition matrices that induce the optimal probability distribution. However, one limitation of the previous work is that randomness can only be introduced via the control policy; all uncontrolled processes must be modeled as deterministic to render them immutable under an exponential tilting. In this work, only the controlled dynamics are subject to tilting, allowing for more general probabilistic models. Numerical experiments are conducted in the context of power networks. The distributed control techniques described in this paper can transform a large collection of flexible loads into a 'virtual battery' capable of delivering the same grid services as traditional batteries. Additionally, quality of service to the load owner is guaranteed, privacy is preserved, and computation and communication requirements are reduced, relative to alternative centralized control techniques.
3. Control oriented modeling of TCLs 20 Thermostatically controlled loads (TCLs) have the potential to be a valuable resource for the Balancing Authority (BA) of the future. Examples of TCLs include household appliances such as air conditioners, water heaters, and refrigerators. Since the rated power of each TCL is on the order of kilowatts, to provide meaningful service for the BA, it is necessary to control large collections of TCLs. To perform design of a distributed coordination/control algorithm, the BA requires a control oriented model that describes the relevant dynamics of an ensemble. Works focusing on solely modeling the ensemble date back to the 1980's, while works focusing on control oriented modeling are more recent. In this work, we contribute to the control oriented modeling literature. We leverage techniques from computational fluid dynamics (CFD) to discretize a pair of FokkerPlanck equations derived in earlier work. The discretized equations are shown to admit a certain factorization, which makes the developed model useful for control design. In particular, the effects of weather and control are shown to independently effect the system dynamics.
4. A unified framework for coordination of thermostatically controlled loads 35 A collection of thermostatically controlled loads (TCLs) – such as air conditioners and water heaters – can vary their power consumption within limits to help the balancing authority of a power grid maintain demand supply balance. Doing so requires loads to coordinate their on/off decisions so that the aggregate power consumption profile tracks a gridsupplied reference. At the same time, each consumer's quality of service (QoS) must be maintained. While there is a large body of work on TCL coordination, there are several limitations. One is that they do not provide guarantees on the reference tracking performance and QoS maintenance. A second limitation of past work is that they do not provide a means to compute a suitable reference signal for power demand of a collection of TCLs. In this work we provide a framework that addresses these weaknesses. The framework enables coordination of an arbitrary number of TCLs that: (i) is computationally efficient, (ii) is implementable at the TCLs with local feedback and low communication, and (iii) enables reference tracking by the collection while ensuring that temperature and cycling constraints are satisfied at every TCL at all times. The framework is based on a Markov model obtained by discretizing a pair of FokkerPlanck equations derived in earlier work by Malhame and Chong. We then use this model to design randomized policies for TCLs. The balancing authority broadcasts the same policy to all TCLs, and each TCL implements this policy which requires only local measurement to make on/off decisions. Simulation results are provided to support these claims.
5. A Generalized Nash Equilibrium analysis of the interaction between a peertopeer financial market and the distribution grid 28 We consider the interaction between the distribution grid (physical level) managed by the distribution system operator (DSO), and a financial market in which prosumers optimize their demand, generation, and bilateral trades in order to minimize their costs subject to local constraints and bilateral trading reciprocity coupling constraints. We model the interaction problem between the physical and financial levels as a noncooperative generalized Nash equilibrium problem. We compare two designs of the financial level prosumer market: a centralized design and a peertopeer fully distributed design. We prove the Pareto efficiency of the equilibria under homogeneity of the trading cost preferences. In addition, we prove that the pricing structure of our noncooperative game does not permit freelunch behavior. Finally, in the numerical section we provide additional insights on the efficiency loss with respect to the different levels of agents' flexibility and amount of renewables in the network. We also quantify the impact of the prosumers' pricing on the noncooperative game social cost.
6. Privacy Impact on Generalized Nash Equilibrium in PeertoPeer Electricity Market 15 We consider a peertopeer electricity market, where agents hold private information that they might not want to share. The problem is modeled as a noncooperative communication game, which takes the form of a Generalized Nash Equilibrium Problem, where the agents determine their randomized reports to share with the other market players, while anticipating the form of the peertopeer market equilibrium. In the noncooperative game, each agent decides on the deterministic and random parts of the report, such that (a) the distance between the deterministic part of the report and the truthful private information is bounded and (b) the expectation of the privacy loss random variable is bounded. This allows each agent to change her privacy level. We characterize the equilibrium of the game, prove the uniqueness of the Variational Equilibria and provide a closed form expression of the privacy price. In addition, we provide a closed form expression to measure the impact of the privacy preservation caused by inclusion of random noise and deterministic deviation from agents' true values. Numerical illustrations are presented on the 14bus IEEE network.
7. A product form for the general stochastic matching model 12 We consider a stochastic matching model with a general compatibility graph, as introduced by Mairesse and Moyal (2016). We show that the natural necessary condition of stability of the system is also sufficient for the natural ‘firstcome, firstmatched’ matching policy. To do so, we derive the stationary distribution under a remarkable product form, by using an original dynamic reversibility property related to that of Adan, Bušić, Mairesse, and Weiss (2018) for the bipartite matching model.
8. Stochastic dynamic matching: A mixed graphtheory and linearalgebra approach 36 The stochastic dynamic matching problem has recently drawn attention in the stochasticmodeling community due to its numerous applications, ranging from supplychain management to kidney exchange programs. In this paper, we consider a matching problem in which items of different classes arrive according to independent Poisson processes. Unmatched items are stored in a queue, and compatibility constraints are described by a simple graph on the classes, so that two items can be matched if their classes are neighbors in the graph. We analyze the efficiency of matching policies, not only in terms of system stability, but also in terms of matching rates between different classes. Our results rely on the observation that, under any stable policy, the matching rates satisfy a conservation equation that equates the arrival and departure rates of each item class. Our main contributions are threefold. We first introduce a mapping between the dimension of the solution set of this conservation equation, the structure of the compatibility graph, and the existence of a stable policy. In particular, this allows us to derive a necessary and sufficient stability condition that is verifiable in polynomial time. Secondly, we describe the convex polytope of nonnegative solutions of the conservation equation. When this polytope is reduced to a single point, we give a closedform expression of the solution; in general, we characterize the vertices of this polytope using again the graph structure. Lastly, we show that greedy policies cannot, in general, achieve every point in the polytope. In contrast, nongreedy policies can reach any point of the interior of this polytope, and we give a condition for these policies to also reach the boundary of the polytope.
9. A Multiclass Energy Packet Networks with finite capacity energy queues 14 Energy packet Network (EPN) consists of a queueing network formed by blocks, where each of them is formed by one data queue, that handles the workload, and one energy queue, that handles packets of energy. We study an EPN model where the energy packets start the transfer. In this model, energy packets are sent to the data queue of the same block. An energy packet routes one workload packet to the next block if the data queue is not empty, and it is lost otherwise. We assume that the energy queues have a finite buffer size and if an energy packet arrives to the system when the buffer is full, jumpover blocking (JOB) is performed, and therefore with some probability it is sent to the data queue and it is lost otherwise. We first provide a value of the jumpover blocking probability such that the steadystate probability distribution of packets in the queues admits a product form solution. The product form is established for multiserver and multiclass data packet queues under FCFS, preemptive LCFS and PS discipline. Moreover, in the case of a directed tree queueing network, we show that the number of data packets in each subtree decreases as the JOB probability increases for each block.
7.3 Reinforcement learning
10. Fundamental Design Principles for Reinforcement Learning Algorithms 29 Along with the sharp increase in visibility of the field, the rate at which new reinforcement learning algorithms are being proposed is at a new peak. While the surge in activity is creating excitement and opportunities, there is a gap in understanding of two basic principles that these algorithms need to satisfy for any successful application. One has to do with guarantees for convergence, and the other concerns the convergence rate. The vast majority of reinforcement learning algorithms belong to a class of learning algorithms known as stochastic approximation (SA). The objective here is to review the foundations of reinforcement learning algorithm design based on recent and ancient results from SA. In particular, it was established in (Borkar and Meyn, 2000) that both stability and convergence of these algorithms are guaranteed by analyzing the stability of two associated ODEs. Moreover, if the linearized ODE passes a simple eigenvalue test, then an optimal rate of convergence is guaranteed. This chapter contains a survey of these concepts, along with a survey of the new class of Zap reinforcement learning algorithms introduced by the authors. These algorithms can achieve convergence almost universally, while also guaranteeing optimal rate of convergence.
7.4 Mathematics of wireless cellular networks
11. Maximum UtilityAware Capacity Partitioning in Cooperative Computing 16 In many networks, a user has to allocate the link capacity between upload and download. When such networks are used for cooperative computing, the user needs to maintain the division of upload and download capacities at an optimal value to receive the maximum utility. To determine this optimal value, we model uploaddownload partitioning as a resource maximization game. We show that a Nash equilibrium (NE) obtained for this game is socially optimal. Thus this NE acts as an upper bound on capacity partitioning and serves as a benchmark to analyze the efficiency and performance of various capacity partitioning algorithms. Specifically, using this upper bound and simulations, we examine the performance of different partitioning algorithms while considering the dynamics of resource requests.
12. Where to Deploy Reconfigurable Intelligent Surfaces in the Presence of Blockages? 27 Wireless communications aided by reconfigurable intelligent surfaces (RISs) is a promising way to improve the coverage for cellular users. The controlled reflection of the signal from RISs is especially useful in mmwave networks when the direct link between a cellular user and its serving base station (BS) is weak or unavailable due to blockages. But the joint blockage of the userRIS and the userBS links may significantly degrade the performance of RISaided transmissions. This paper aims to study the effect of joint blockages on downlink performance. When the RIS locations are coupled with BS locations, using tools from stochastic geometry, we obtain an optimal placement of RISs either to minimize the joint blockage probability of the userRIS and the userBS links or to maximize the downlink coverage probability. The results show that installing RISs on the street intersections improves the coverage probability. For users associated with BSs that are deployed sufficiently close to intersections, the intersectionmounted RISs offer a better coverage performance compared to BScoupled RISs.
13. Beam Management in Cellular Networks 11 Beam management is central in the operation of beamformed wireless cellular systems such as 5G New Radio (NR) networks. Focusing the energy radiated to mobile terminals (MTs) by increasing the number of beams per cell increases signal power and decreases interference, and has hence the potential to bring major improvements on area spectral efficiency (ASE). We proposed a first systemlevel stochastic geometry model encompassing major aspects of the beam management problem: frequencies, antenna configurations, and propagation; physical layer, wireless links, and coding; network geometry, interference, and resource sharing; sensing, signaling, and mobility management. This model leads to a simple analytical expression for the effective rate that the typical user gets in this context. This in turn allows one to find the number of beams per cell and per MT that maximizes the effective ASE by offering the best tradeoff between beamforming gains and beam management operational overheads and costs, for a wide variety of 5G network scenarios including millimeter wave (mmWave) and sub6 GHz. As part of the systemlevel analysis, we define and analyze several underlying new and fundamental performance metrics that are of independent interest. The numerical results discuss the effects of different systemic tradeoffs and performance optimizations of mmWave and sub6 GHz 5G deployments.
14. A FineGrained Analysis of Radar Detection in Vehicular Networks 26 Automotive radar is a critical feature in advanced driverassistance systems. It is important in enhancing vehicle safety by detecting the presence of other vehicles in the vicinity. The performance of radar detection is, however, affected by the interference from radars of other vehicles as well as the variation in the target radar crosssection (RCS) due to varying physical features of the target vehicle. Considering such interference and random RCS, this work provides a finegrained performance analysis of radar detection. Specifically, using stochastic geometry, we calculate the meta distribution of the signaltointerferenceandnoise ratio that permits the reliability analysis of radar detection at individual vehicles. We also evaluate the delay aspect of radar detection, namely, the mean local delay which is the average number of transmission attempts needed until the first successful target detection. For a given target distance, we obtain the optimal transmit probability that maximizes the density of successful radar detection while keeping the mean local delay below a threshold. We also provide several system design insights in terms of the fraction of reliable radar links, transmission delay, the density of vehicles, and congestion control.
15. Performance of vehicletoeverything (V2X) communications 8 We analyzed the broadcast of safetyrelated V2X communications in cellular networks where base stations and vehicles are assumed to share the same spectrum and vehicles broadcast their safety messages to neighboring users. We model the locations of vehicles as a Poisson line Cox point process and the locations of users as a planar Poisson point process. We assume that users are associated with their closest base stations when there is no vehicle within a certain distance $\rho $. On the other hand, users located within a distance $\rho $ from vehicles are associated with the vehicles to receive their safety messages. We quantify the properties of this vehicleprioritized association using the stochastic geometry framework. We derive the fractions of users that receive safety messages from vehicles. Then, we obtain the expression for the signaltointerference ratio of the typical user evaluated on each association type. To address the impact of vehicular broadcast on the cellular network, we also derived the effective rate offered to the typical user in this setting
16. Randomised Geographic Caching and its Applications in Wireless Networks 30 The randomised (or probabilistic) geographic caching is a proactive content placement strategy that has attracted a lot of attention, because it can simplify a great deal cachemanagement problems at the wireless edge. It diversifies content placement over caches and applies to scenarios where a request can be possibly served by multiple cache memories. Its simplicity and strength is due to randomisation. It allows one to formulate continuous optimisation problems for content placement over large homogeneous geographic areas. These can be solved to optimality by standard convex methods, and can even provide closedform solutions for specific cases. This way the algorithmic obstacles from NPhardness are avoided and optimal solutions can be derived with low computational cost. Randomised caching has a large spectrum of applications in realworld wireless problems, including femtocaching, multitier networks, devicetodevice communications, mobility, mmwave, security, UAVs, and more. In this chapter we will formally present the main policy with its applications in various wireless scenarios. We will further introduce some very useful extensions related to unequal filesizes and content placement with neighbourhood dependence.
17. How Wireless Queues Benefit from Motion 13 This paper considers the time evolution of a queue that is embedded in a Poisson point process of moving wireless interferers. The queue is driven by an external arrival process and is subject to a timevarying service process that is a function of the SINR that it sees. Static configurations of interferers result in an infinite queue workload with positive probability. In contrast, a generic stability condition is established for the queue in the case where interferers possess any nonzero mobility that results in displacements that are both independent across interferers and oblivious to interferer positions. The proof leverages the mixing property of the Poisson point process. The effect of an increase in mobility on queuing metrics is also studied. Convex ordering tools are used to establish that faster moving interferers result in a queue workload that is smaller for the increasingconvex stochastic order. As a corollary, mean workload and mean delay decrease as network mobility increases. This stochastic ordering as a function of mobility is explained by establishing positive correlations between SINR levelcrossing events at different time points, and by determining the autocorrelation function for interference and observing that it decreases with increasing mobility. System behaviour is empirically analyzed using discreteevent simulation and the performance of various mobility models is evaluated using heavytraffic approximations.
18. On Velocitybased Association Policies for Multitier 5G Wireless Networks 43 Mobility is a key challenge for beam management in 5G cellular networks due to the overhead incurred at beam switching and base station (BS) handover events. This paper focuses on a network that has a multitier structure with two types of BSs operating in the same frequency bands, namely macro BSs that are sparser but with higher transmit power, and micro BSs that are denser and with lower transmit power. We propose a downlink user association policy which is a function of the user mobility. Typically, high mobility users should associate with macro BSs so as to incur less beam switching overhead, whereas low mobility ones should be associated with micro BSs. The main contribution of the paper is a formalization of the optimal threshold association policy, when the optimality is understood with respect to the effective Shannon rate. The analysis is based on stochastic geometry and on an exact representation of the effective Shannon rate of the typical user in this beamforming multitier context. Two models are discussed. The simplest one focuses on a singleuser optimization problem. We also discuss a more realistic model with bandwidth sharing between all users in the cell. Finally, we identify the mobility and userdensity patterns where the velocitybased threshold association policy outperforms the classical best mean power association policy.
19. Energy and Delay TradeOffs of EndtoEnd Vehicular Communications using a Hyperfractal Urban Modelling 33 We characterize tradeoffs between the endtoend communication delay and the energy in urban vehicular communications with infrastructure assistance. Our study exploits the selfsimilarity of the location of communication entities in cities by modeling them with the hyperfractal model which charaterize the distribution of mobile nodes and relay nodes by a fractal dimension ${d}_{F}$ and ${d}_{r}$ , both larger than the dimension of the embedded map. We compute theoretical bounds for the endtoend communication hop count considering two different energyminimizing goals: either total accumulated energy or maximum energy per node. Let $\delta >1$ the attenuation factor in the street, we prove that when we aim to a total energy cost of order ${n}^{(1\delta )(1\alpha )}$ the hop count for an endtoend transmission is of order ${n}^{1\alpha /({d}_{F}1)}$ , with $\alpha <1$ is a tunable parameter. This proves that for both goals the energy decreases as we allow choosing routing paths of higher length. The asymptotic limit of the energy becomes significantly small when the number of nodes becomes asymptotically large. A lower bound on the network throughput capacity with constraints on path energy is also given. We show that our model fits real deployments where open data sets are available. The results are confirmed through simulations using different fractal dimensions in a Matlab simulator.
7.5 Highdimensional statistical inference
20. Emergence of extended states at zero in the spectrum of sparse random graphs 9 We confirm the longstanding prediction that $c=e\approx 2.718$ is the threshold for the emergence of a nonvanishing absolutely continuous part (extended states) at zero in the limiting spectrum of the ErdősRenyi random graph with average degree $c$. This is achieved by a detailed secondorder analysis of the resolvent ${(Az)}^{1}$ near the singular point $z=0$, where $A$ is the adjacency operator of the PoissonGaltonWatson tree with mean offspring $c$. More generally, our method applies to arbitrary unimodular GaltonWatson trees, yielding explicit criteria for the presence or absence of extended states at zero in the limiting spectral measure of a variety of random graph models, in terms of the underlying degree distribution.
21. Sparse matrices: convergence of the characteristic polynomial seen from infinity 37 We prove that the reverse characteristic polynomial $det({I}_{n}z{A}_{n})$ of a random $n\times n$ matrix ${A}_{n}$ with iid $\mathrm{Bernoulli}(d/n)$ entries converges in distribution towards the random infinite product ${\prod}_{\ell =1}^{\infty}{(1{z}^{\ell})}^{{Y}_{\ell}}$ where ${Y}_{\ell}$ are independent $\mathrm{Poisson}({d}^{\ell}/\ell )$ random variables. We show that this random function is a Poisson analog of more classical Gaussian objects such as the Gaussian holomorphic chaos. As a byproduct, we obtain new simple proofs of previous results on the asymptotic behaviour of extremal eigenvalues of sparse ErdősRényi digraphs: for every $d>1$, the greatest eigenvalue of ${A}_{n}$ is close to $d$ and the second greatest is smaller than $\sqrt{d}$, a Ramanujanlike property for irregular digraphs. For $d<1$, the only nonzero eigenvalues of ${A}_{n}$ converge to a Poisson multipoint process on the unit circle. Our results also extend to the semisparse regime where $d$ is allowed to grow to $\infty $ with $n$, slower than ${n}^{o\left(1\right)}$. We show that the reverse characteristic polynomial converges towards a more classical object written in terms of the exponential of a logcorrelated real Gaussian field. In the semisparse regime, the empirical spectral distribution of ${A}_{n}/\sqrt{{d}_{n}}$ converges to the circle distribution; as a consequence of our results, the second eigenvalue sticks to the edge of the circle.
22. Expressive Power of Invariant and Equivariant Graph Neural Networks 18 Various classes of Graph Neural Networks (GNN) have been proposed and shown to be successful in a wide range of applications with graph structured data. In this paper, we propose a theoretical framework able to compare the expressive power of these GNN architectures. The current universality theorems only apply to intractable classes of GNNs. Here, we prove the first approximation guarantees for practical GNNs, paving the way for a better understanding of their generalization. Our theoretical results are proved for invariant GNNs computing a graph embedding (permutation of the nodes of the input graph does not affect the output) and equivariant GNNs computing an embedding of the nodes (permutation of the input permutes the output). We show that Folklore Graph Neural Networks (FGNN), which are tensor based GNNs augmented with matrix multiplication are the most expressive architectures proposed so far for a given tensor order. We illustrate our results on the Quadratic Assignment Problem (a NPHard combinatorial problem) by showing that FGNNs are able to learn how to solve the problem, leading to much better average performances than existing algorithms (based on spectral, SDP or other GNNs architectures). On a practical side, we also implement masked tensors to handle batches of graphs of varying sizes.
23. Impossibility of Partial Recovery in the Graph Alignment Problem 24 Random graph alignment refers to recovering the underlying vertex correspondence between two random graphs with correlated edges. This can be viewed as an averagecase and noisy version of the wellknown graph isomorphism problem. For the correlated ErdősRényi model, we prove an impossibility result for partial recovery in the sparse regime, with constant average degree and correlation, as well as a general bound on the maximal reachable overlap. Our bound is tight in the noiseless case (the graph isomorphism problem) and we conjecture that it is still tight with noise. Our proof technique relies on a careful application of the probabilistic method to build automorphisms between tree components of a subcritical ErdősRényi graph.
24. Sharp threshold for alignment of graph databases with Gaussian weights 25 We study the fundamental limits for reconstruction in weighted graph (or matrix) database alignment. We consider a model of two graphs where ${\pi}^{*}$ is a planted uniform permutation and all pairs of edge weights ${({A}_{i,j},{B}_{{\pi}^{*}\left(i\right),{\pi}^{*}\left(j\right)})}_{1\le i<j\le n}$ are i.i.d. pairs of Gaussian variables with zero mean, unit variance and correlation parameter $\rho \in [0,1]$. We prove that there is a sharp threshold for exact recovery of ${\pi}^{*}$: if $n{\rho}^{2}\ge (4+\u03f5)logn+\omega \left(1\right)$ for some $\u03f5>0$, there is an estimator $\widehat{\pi}$ – namely the MAP estimator – based on the observation of databases $A,B$ that achieves exact reconstruction with high probability. Conversely, if $n{\rho}^{2}\le 4lognloglogn\omega \left(1\right)$, then any estimator $\widehat{\pi}$ verifies $\widehat{\pi}=\pi $ with probability $o\left(1\right)$. This result shows that the informationtheoretic threshold for exact recovery is the same as the one obtained for detection in a recent work by Wu et al. (2020): in other words, for Gaussian weighted graph alignment, the problem of reconstruction is not more difficult than that of detection. Though the reconstruction task was already well understood for vectorshaped database alignment (that is taking signal of the form ${({u}_{i},{v}_{{\pi}^{*}\left(i\right)})}_{1\le i\le n}$ where $({u}_{i},{v}_{{\pi}^{*}\left(i\right)})$ are i.i.d. pairs in ${\mathbb{R}}^{{d}_{u}}\times {\mathbb{R}}^{{d}_{v}}$), its formulation for graph (or matrix) databases brings a drastically different problem for which the hard phase is conjectured to be wide. The proofs build upon the analysis of the MAP estimator and the second moment method, together with the study of the correlation structure of energies of permutations.
25. A simpler spectral approach for clustering in directed network 38 We study the task of clustering in directed networks. We show that using the eigenvalue/eigenvector decomposition of the adjacency matrix is simpler than all common methods which are based on a combination of data regularization and SVD truncation, and works very well down to the very sparse regime where the edge density has constant order. This simple approach was largely unnoticed in the mathematics and network science communities. Our analysis is based on a Master Theorem describing sharp asymptotics for isolated eigenvalues/eigenvectors of sparse, nonsymmetric matrices with independent entries. We also describe the limiting distribution of the entries of these eigenvectors; in the task of digraph clustering with spectral embeddings, we provide numerical evidence for the superiority of Gaussian Mixture clustering over the widely used kmeans algorithm.
26. Accelerating Abelian Random Walks with Hyperbolic Dynamics 39 Given integers $d\ge 2,n\ge 1$, we consider affine random walks on torii ${(\mathbb{Z}/n\mathbb{Z})}^{d}$ defined as ${X}_{t+1}=A{X}_{t}+{B}_{t}\phantom{\rule{0.166667em}{0ex}}\mathrm{mod}\phantom{\rule{0.166667em}{0ex}}n$, where $A\in {\mathrm{GL}}_{d}\left(\mathbb{Z}\right)$ is an invertible matrix with integer entries and ${\left({B}_{t}\right)}_{t\ge 0}$ is a sequence of iid random increments on ${\mathbb{Z}}^{d}$. We show that when $A$ has no eigenvalues of modulus 1, this random walk mixes in $O(lognloglogn)$ steps as $n\to \infty $, and mixes actually in $O(logn)$ steps only for almost all $n$. These results generalize those on the socalled ChungDiaconisGraham process, which corresponds to the case $d=1$. Our proof is based on the initial arguments of Chung, Diaconis and Graham, and relies extensively on the properties of the dynamical system $x\mapsto {A}^{\top}x$ on the continuous torus ${\mathbb{R}}^{d}/{\mathbb{Z}}^{d}$. Having no eigenvalue of modulus one makes this dynamical system a hyperbolic toral automorphism, a typical example of a chaotic system known to have a rich behaviour. As such our proof sheds new light on the speedup gained by applying a deterministic map to a Markov chain.
27. Asynchrony and Acceleration in Gossip Algorithms 41 This paper considers the minimization of a sum of smooth and strongly convex functions dispatched over the nodes of a communication network. Previous works on the subject either focus on synchronous algorithms, which can be heavily slowed down by a few slow nodes (the straggler problem), or consider a historical asynchronous setting (Boyd et al., 2006), which relies on a communication model that cannot be readily implemented in practice, as it does not capture important aspects of asynchronous communications such as noninstantaneous computations and communications. We have two main contributions. 1) We introduce a new communication scheme, based on LossNetworks, that is programmable in a fully asynchronous and decentralized fashion. We establish empirically and theoretically that it improves over existing synchronous algorithms by depending on local communication delays in the analysis instead of global worstones. 2) We provide an acceleration of the standard gossip algorithm in the historical asynchronous model without requiring any additional synchronization.
28. Nonbacktracking spectra of weighted inhomogeneous random graphs 44 We study a model of random graphs where each edge is drawn independently (but not necessarily identically distributed) from the others, and then assigned a random weight. When the mean degree of such a graph is low, it is known that the spectrum of the adjacency matrix $A$ deviates significantly from that of its expected value $\mathbb{E}A$. In contrast, we show that over a wide range of parameters the top eigenvalues of the nonbacktracking matrix $B$ – a matrix whose powers count the nonbacktracking walks between two edges – are close to those of $\mathbb{E}A$, and all other eigenvalues are confined in a bulk with known radius. We also obtain a precise characterization of the scalar product between the eigenvectors of $B$ and their deterministic counterparts derived from the model parameters. This result has many applications, in domains ranging from (noisy) matrix completion to community detection, as well as matrix perturbation theory. In particular, we establish as a corollary that a result known as the BaikBen ArousPéché phase transition, previously established only for rotationally invariant random matrices, holds more generally for matrices $A$ as above under a mild concentration hypothesis.
7.6 Distributed optimization for machine learning
29. Fast Stochastic Bregman Gradient Methods: Sharp Analysis and Variance Reduction 21 We study the problem of minimizing a relativelysmooth convex function using stochastic Bregman gradient methods. We first prove the convergence of Bregman Stochastic Gradient Descent (BSGD) to a region that depends on the noise (magnitude of the gradients) at the optimum. In particular, BSGD with a constant stepsize converges to the exact minimizer when this noise is zero (interpolation setting, in which the data is fit perfectly). Otherwise, when the objective has a finite sum structure, we show that variance reduction can be used to counter the effect of noise. In particular, fast convergence to the exact minimizer can be obtained under additional regularity assumptions on the Bregman reference function. We illustrate the effectiveness of our approach on two key applications of relative smoothness: tomographic reconstruction with Poisson noise and statistical preconditioning for distributed optimization.
30. A Continuized View on Nesterov Acceleration for Stochastic Gradient Descent and Randomized Gossip 22 We introduce the "continuized" Nesterov acceleration, a close variant of Nesterov acceleration whose variables are indexed by a continuous time parameter. The two variables continuously mix following a linear ordinary differential equation and take gradient steps at random times. This continuized variant benefits from the best of the continuous and the discrete frameworks: as a continuous process, one can use differential calculus to analyze convergence and obtain analytical expressions for the parameters; and a discretization of the continuized process can be computed exactly with convergence rates similar to those of Nesterov original acceleration. We show that the discretization has the same structure as Nesterov acceleration, but with random parameters. We provide continuized Nesterov acceleration under deterministic as well as stochastic gradients, with either additive or multiplicative noise. Finally, using our continuized framework and expressing the gossip averaging problem as the stochastic minimization of a certain energy function, we provide the first rigorous acceleration of asynchronous gossip algorithms.
31. Concentration of NonIsotropic Random Tensors with Applications to Learning and Empirical Risk Minimization 23 Dimension is an inherent bottleneck to some modern learning tasks, where optimization methods suffer from the size of the data. In this paper, we study nonisotropic distributions of data and develop tools that aim at reducing these dimensional costs by a dependency on an effective dimension rather than the ambient one. Based on nonasymptotic estimates of the metric entropy of ellipsoidsthat prove to generalize to infinite dimensionsand on a chaining argument, our uniform concentration bounds involve an effective dimension instead of the global dimension, improving over existing results. We show the importance of taking advantage of nonisotropic properties in learning problems with the following applications: i) we improve stateoftheart results in statistical preconditioning for communicationefficient distributed optimization, ii) we introduce a nonisotropic randomized smoothing for nonsmooth optimization. Both applications cover a class of functions that encompasses empirical risk minization (ERM) for linear models.
7.7 Stochastic Geometry
32. Optimal stationary markings 7 Many specific problems ranging from theoretical probability to applications in statistical physics, combinatorial optimization and communications can be formulated as an optimal tuning of local parameters in large systems of interacting particles. Using the framework of stationary point processes in the Euclidean space, we pose it as a problem of an optimal stationary marking of a given stationary point process. The quality of a given marking is evaluated in terms of scores calculated in a covariant manner for all points in function of the proposed marked configuration. In the absence of total order of the configurations of scores, we identify intensityoptimality and local optimality as two natural ways for defining optimal stationary marking. We derive tightness and integrability conditions under which intensityoptimal markings exist and further stabilization conditions making them equivalent to locally optimal ones. We present examples motivating the proposed, general framework. Finally, we discuss various possible approaches leading to uniqueness results.
33. Continuum LineofSight Percolation on PoissonVoronoi Tessellations 10 In this work, we study a new model for continuum lineofsight percolation in a random environment given by a PoissonVoronoi tessellation. The edges of this tessellation are the support of a Cox point process, while the vertices are the support of a Bernoulli point process. Taking the superposition $\mathbb{Z}$ of these two processes, two points of $\mathbb{Z}$ are linked by an edge if and only if they are sufficiently close and located on the same edge of the supporting tessellation. We study the percolation of the random graph arising from this construction and prove that a subcritical phase as well as a supercritical phase exist under general assumptions. Our proofs are based on a renormalization argument with some notion of stabilization and asymptotic essential connectedness to investigate continuum percolation for Cox point processes. We also give numerical estimates of the critical parameters of the model. Our model can be seen as a good candidate for modelling telecommunications networks in a random environment with obstructive conditions for signal propagation.
34. On Point Processes Defined by Angular Conditions on Delaunay Neighbors in the PoissonVoronoi Tessellation 3 Consider a homogeneous Poisson point process of the Euclidean plane and its Voronoi tessellation. The present note discusses the properties of two stationary point processes associated with the latter and depending on a parameter $\theta $. The first one is the set of points that belong to some onedimensional facet of the Voronoi tessellation and are such that the angle with which they see the two nuclei defining the facet is $\theta $. The main question of interest on this first point process is its intensity. The second point process is that of the intersections of the said tessellation with a straight line having a random orientation. Its intensity is well known. The intersection points almost surely belong to onedimensional facets. The main question here is about the Palm distribution of the angle with which the points of this second point process see the two nuclei associated with the facet. The note gives answers to these two questions and briefly discusses their practical motivations. It also discusses natural extensions to dimension three.
35. Unimodular Hausdorff and Minkowski dimensions 2 We introduced two new notions of dimension, namely the unimodular Minkowski and Hausdorff dimensions, which are inspired from the classical analogous notions. These dimensions are defined for unimodular discrete spaces, introduced in this work, which provide a common generalization to stationary point processes under their Palm version and unimodular random rooted graphs. The use of unimodularity in the definitions of dimension is novel. Also, a toolbox of results was presented for the analysis of these dimensions. In particular, analogues of Billingsley's lemma and Frostman's lemma are presented. These last lemmas are instrumental in deriving upper bounds on dimensions, whereas lower bounds are obtained from specific coverings. The notions of unimodular Hausdorff size, which is a discrete analogue of the Hausdorff measure, and unimodular dimension function are also introduced. This toolbox allows one to connect the unimodular dimensions to other notions such as volume growth rate, discrete dimension and scaling limits. It is also used to analyze the dimensions of a set of examples pertaining to point processes, branching processes, random graphs, random walks, and selfsimilar discrete random spaces. Further results of independent interest are also presented, like a version of the maxflow mincut theorem for unimodular oneended trees and a weak form of pointwise ergodic theorems for all unimodular discrete spaces.
7.8 Mathematics of stochastic networks via meanfield analysis
36. Pair Replica MeanField Neural Networks 5 Replicameanfield models have been proposed to decipher the activity of otherwise analytically intractable neural networks via a multiplyandconquer approach. In this approach, one considers limit networks made of infinitely many replicas with the same basic neural structure as that of the network of interest, but exchanging spikes in a randomized manner. The key point is that these replicameanfield networks are tractable versions that retain important features of the finite structure of interest. To date, the replica framework has been discussed for firstorder models, whereby elementary replica constituents are single neurons with independent Poisson inputs. In 5, we extend this replica framework to allow elementary replica constituents to be composite objects, namely, pairs of neurons. As they include pairwise interactions, these pairreplica models exhibit nontrivial dependencies in their stationary dynamics, which cannot be captured by firstorder replica models. Our contributions are twofold: (i) We analytically characterize the stationary dynamics of a pair of intensitybased neurons with independent Poisson input. This analysis involves the reduction of a boundaryvalue problem related to a twodimensional transport equation to a system of Fredholm integral equations—a result of independent interest. (ii) We analyze the set of consistency equations determining the full network dynamics of certain replica limits. These limits are those for which replica constituents, be they single neurons or pairs of neurons, form a partition of the network of interest. Both analyses are numerically validated by computing input/output transfer functions for neuronal pairs and by computing the correlation structure of certain pairdominated network dynamics.
37. Replicameanfield limits of fragmentationinteractionaggregation processes 1 Network dynamics with pointprocessbased interactions are of paramount modeling interest. Unfortunately, most relevant dynamics involve complex graphs of interactions for which an exact computational treatment is impossible. To circumvent this difficulty, the replicameanfield approach focuses on randomly interacting replicas of the networks of interest. In the limit of an infinite number of replicas, these networks become analytically tractable under the socalled ‘Poisson hypothesis’. However, in most applications this hypothesis is only conjectured. In this paper we establish the Poisson hypothesis for a general class of discretetime, pointprocessbased dynamics that we propose to call fragmentationinteractionaggregation processes, and which are introduced here. These processes feature a network of nodes, each endowed with a state governing their random activation. Each activation triggers the fragmentation of the activated node state and the transmission of interaction signals to downstream nodes. In turn, the signals received by nodes are aggregated to their state. Our main contribution is a proof of the Poisson hypothesis for the replicameanfield version of any network in this class. The proof is obtained by establishing the propagation of asymptotic independence for state variables in the limit of an infinite number of replicas. Discretetime Galves–Löcherbach neural networks are used as a basic instance and illustration of our analysis.
38. Opinion Dynamics 6 We introduced a nonlinear and continuoustime opinion dynamics model with additive noise and statedependent interaction rates between agents. The model features interaction rates which are proportional to a negative power of the opinion distances. We establish a nonlocal partial differential equation for the distribution of opinion distances and use Mellin transforms to provide an explicit formula for the stationary solution of the latter, when it exists. Our approach leads to new qualitative and quantitative results on this type of dynamics. To the best of our knowledge these Mellin transform results are the first quantitative results on the equilibria of opinion dynamics with distancedependent interaction rates. The closedform expressions for this class of dynamics are obtained for the twoagent case. However, the results can be used in meanfield models featuring several agents whose interaction rates depend on the empirical average of their opinions. The technique also applies to linear dynamics, namely with a constant interaction rate, on an interaction graph.
39. Stochastic Modelling of FreeFloating CarSharing Systems 40 Carsharing systems (CSSs) have gained popularity during the last decade as a flexible, efficient and ecological alternative mode of transportation. But for the operator, managing such systems is far to be simple. Due to heterogeneity of demand and also randomness, the user may face a lack of resources: no car or no parking space available. And the operator has to design the system in order to improve it. The total number of cars impact the performance of the system. We address the dimensioning issue. For that, mathematical models are needed. In many cities, two systems coexist: stationbased and freefloating. The latter gives more flexibility to the user both to take or return the car. But he can reserve only the car for a short period, and not the parking space, as the car is parked on public space with no specific parking spaces. The car reservation is here to help the user. The aim of the paper is to study its influence on the system behavior. This study focuses on Communauto's Montreal freefloating carsharing system (FFCSS). Data analysis investigates the main features of the system based on user preferences. It allows proposing a mathematical modelling. Then we present two analytical approaches. First the meanfield method could be used for different variants, and w gives first insights on the optimal fleet size in a homogeneous framework. Second the general inhomogeneous model is described as a closed Jackson network with blockingrerouting policy. We prove that its state at stationarity is given by a productform distribution. It allows in future work to obtain an explicit largescale representation of the system which can be used both theoretically or numerically for optimization purposes.
40. Mean field analysis for bike and ebike sharing systems 17 Abstract : Electric bikes are deployed massively in preexisting bike sharing system in order to attract new users and replace cars on a larger scale (see [2]). But this causes interactions between the two populations of bikes. In this paper, we analyze a model of an homogeneous bike sharing system where two classes of bikes interact only through the finite capacity of stations. It models systems with both electric and normal bikes, these classes requiring different subscriptions. As far as we know (see [7]), it is the first stochastic largescale analysis for integrated ebike and bike sharing systems. The aim of the paper is to derive explicitly the limiting stationary distribution of the state of a station when the number of stations and the fleet size of each class increase at the same rate. Analysis for a spatially heterogeneous network is in preparation and discussed in Section 4.
41. Mean field analysis of stochastic networks with reservation 34 The problem of reservation in a large distributed system is analyzed via a new mathematical model. A typical application is a stationbased carsharing system which can be described as a closed stochastic network where the nodes are the stations and the customers are the cars. The user can reserve the car and the parking space. In the paper, we study the evolution of the system when the reservation of parking spaces and cars is effective for all users. The asymptotic behavior of the underlying stochastic network is given when the number $N$ of stations and the fleet size increase at the same rate. The analysis involves a Markov process on a state space with dimension of order ${N}^{2}$. It is quite remarkable that the state process describing the evolution of the stations, whose dimension is of order $N$, converges in distribution, although not Markov, to an nonhomogeneous Markov process. We prove this meanfield convergence. We also prove, using combinatorial arguments, that the meanfield limit has a unique equilibrium measure when the time between reserving and picking up the car is sufficiently small. This result extends the case where only the parking space can be reserved.
42. Mean field analysis of an incentive algorithm for a closed stochastic network 42 The paper deals with a loadbalancing algorithm for a closed stochastic network with two zones with different demands. The algorithm is motivated by an incentive algorithm for redistribution of cars in a largescale carsharing system. The service area is divided into two zones. When cars stay too much long in the lowdemand zone, users are encouraged to pick up them and return them in the highdemand zone. The zones are divided in cells called stations. The cars are the network customers. The meanfield limit solution of an ODE gives the large scale distribution of the station state in both clusters for this incentive policy in a discrete Markovian framework. An equilibrium point of this ODE is characterized via the invariant measure of a random walk in the quarterplane. The proportion of empty and saturated stations measures how the system is balanced. Numerical experiments illustrate the impact of the incentive policy. Our study shows that the incentive policy helps when the highdemand zone observes a lack of cars but a saturation must be prevented especially when the highdemand zone is small.
8 Bilateral contracts and grants with industry
8.1 Bilateral contracts with industry
Participants: Bartłomiej Błaszczyszyn, Ana Bušić, Antoine Brochard.
8.1.1 Contract with EDF
Collaborative research in the area of demand dispatch of flexible loads. PI: A. Bušić.
8.1.2 CIFRE with Huawei Technologies France
Contract with Huawei Technologies France started in 2018 and finished in 2021 for the coadvising by B. Błaszczyszyn of a PhD student Antoine Brochard. The PhD has been extended by Inria, who have hired the student until February 2022.
9 Partnerships and cooperations
9.1 ERC NEMO
Participants: François Baccelli, Bartłomiej Błaszczyszyn, Bharath Roy Choudhury, Simon Coste, Ke Feng, Sanket Kalamkar, Sayeh Khaniha, Ali Khezeli, Pierre Popineau.
NEMO, NEtwork MOtion cordis.europa.eu/project/id/788851, project.inria.fr/ercnemo is an ERC Advanced Grant (2019 – 2024, PI François Baccelli). It is an interdisciplinary proposal centered on network dynamics. The interdisciplinarity spans from communication engineering to mathematics, with an innovative interplay between the two. NEMO’s aim is to introduce dynamics in stochastic geometry. General mathematical tools combining stochastic geometry, random graph theory, and the theory of dynamical systems will be developed. NEMO will leverage interactions of Inria with Ecole Normale Supérieure on the mathematical side, and with Nokia Bell Labs and Orange on the engineering side. This year we hired postdoc Ke Feng.
9.2 International initiatives
9.2.1 GSSPPIFCAM
Participants: Bartłomiej Błaszczyszyn.

Title:
Geometric statistics of stationary point processes

Duration:
2018 > 2021

Coordinator:
B. Błaszczyszyn (Inria), Yogeshwaran D. (ISI Bangalore)

Description:
IndoFrench Centre for Applied Mathematics (IFCAM)

Scientific perimeter of the project

Decorrelation concept for general marked point
processes — essential element of limit theory.

Central Limit Theorems for general dependent
interacting particle systems.
 Examples: sequential adsorption, ballistic deposition, majority dynamics, epidemic models — allowing for nonPoisson locations of particles and their dependent initial states.
 Variance asymptotic required for CLTs: volumeorder under quasilocal perturbations and other asymptotic for various test statistics.
 Optimization of large particle systems: intensityoptimal and locally optimal, with applications in statistical physics, combinatorial optimization and communications.

9.2.2 RMF
Participants: François Baccelli.

Title:
Replica Mean Fields for Networks

Duration:
2021 >

Coordinator:
Senya Shlosman ("Senya" <shlosman@gmail.com>;)

Partners:
 Skoltech

Inria contact:
Francois Baccelli
 Summary:
9.3 International research visitors
9.3.1 Visits of international scientists
Inria International Chair
IIC MEYN Sean

Name of the chair:
Sean Meyn

Institution of origin:
University of Florida

Country:
Unites States

Dates:
From Tue Jan 01 2019 to Sun Dec 31 2023

Title:
Distributed Control and Smart Grid

Summary:
Prof. Meyn was at Inria Paris June  July 2021.
9.4 National initiatives
9.4.1 GdR GeoSto
Members of Dyogene participate in Research Group GeoSto (Groupement de recherche, GdR 3477) gdrgeostoch.math.cnrs.fr on Stochastic Geometry led by and David Coupier [Université de Valenciennes].
This is a collaboration framework for all French research teams working in the domain of spatial stochastic modeling, both on theory development and in applications.
9.4.2 GdR RO
Members of Dyogene participate in GdRRO (Recherche Opérationelle; GdR CNRS 3002), gdrro.lip6.fr/, working group COSMOS (Stochastic optimization and control, modeling and simulation), lead by A. Busic and E. Hyon (LIP 6); gdrro.lip6.fr/?q=node/78
9.4.3 ANR JCJC PARI
Probabilistic Approach for Renewable Energy Integration: Virtual Storage from Flexible Loads. The project started in January 2017. PI — A. Bušić. This project is motivated by current and projected needs of a power grid with significant renewable energy integration. Renewable energy sources such as wind and solar have a high degree of unpredictability and time variation, which makes balancing demand and supply challenging. There is an increased need for ancillary services to smooth the volatility of renewable power. In the absence of large, expensive batteries, we may have to increase our inventory of responsive fossilfuel generators, negating the environmental benefits of renewable energy. The proposed approach addresses this challenge by harnessing the inherent flexibility in demand of many types of loads. The objective of the project is to develop decentralized control for automated demand dispatch, that can be used by grid operators as ancillary service to regulate demandsupply balance at low cost. We call the resource obtained from these techniques virtual energy storage (VES). Our goal is to create the necessary ancillary services for the grid that are environmentally friendly, that have low cost and that do not impact the quality of service (QoS) for the consumers. Besides respecting the needs of the loads, the aim of the project is to design local control solutions that require minimal communications from the loads to the centralized entity. This is possible through a systems architecture that includes the following elements: i) local control at each load based on local measurements combined with a gridlevel signal; ii) frequency decomposition of the regulation signal based on QoS and physical constraints for each class of loads.
10 Dissemination
Participants: All Dyogene.
10.1 Promoting scientific activities
10.1.1 Scientific events: organisation
In 2021, A. Khezeli organized the weekly DYOGENE seminar. B. RoyChowdhury organized the NEMO reading group:
 January 22, 2021: "Uniqueness of stationary fixed point for certain exponential server nodes" by Michel Davydov, INRIA.
 March 26, 2021: "Local weak convergence of networks and graph spectra" by Simon Coste, INRIA.
 April 23, 2021: "Unimodular Hausdorff and Minkowski dimension for unimodular random discrete metric spaces, I " by Ali Khezeli, INRIA.
 May 7, 2021: "Unimodular Hausdorff and Minkowski dimension for unimodular random discrete metric spaces, II" by Ali Khezeli, INRIA.
 June 18, 2021: "Free point processes on unimodular second countable groups weakly factor onto Poisson point processes" by Roman Gambelin, INRIA.
 July 2, 2021: "Unimodular Hausdorff and Minkowski dimension for unimodular random discrete metric spaces, III" by Ali Khezeli, INRIA.
 September 30, 2021: "An introduction to the Aldous diffusion conjecture and algebraic measure trees" by Roman Gambelin, INRIA.
 October 14, 2021: " Shift coupling and Embeddings in Brownian motion" by Hermann Thorisson, University of Iceland.
 December 2, 2021: "The $k$parent spatial $\Lambda $Fleming Viot process: construction and limit as k tends to $+\infty $" by Apolline Louvet from MAP5, Université de Paris.
10.1.2 Invited talks
 F. Baccelli:
 Invited lecture at the Hausdorff Institute on High Dimensional Stochastic Geometry www.him.unibonn.de/programs/futureprograms/futuretrimesterprograms/interplayhighdimensionalgeometryprobability/workshophighdimensionalspatialrandomsystemsfebruary22262021
 Colloquium lecture at Georgia Tech, Department of Operations Research on replica meanfield limits for neural networks, March 16, 2021.
 Invited lecture at IPIT, Moscow on replica meanfield limits for neural networks, June 1, 2021.
 Invited lecture at the University of Melbourne, Department of Mathematics, on replica meanfield limits for neural networks, June 10, 2021.
 Keynote lecture at QEST, on Dynamics of Wireless Networks, Paris, August 25. 2021.
 M. Davydov: presentations on mean field limits at
 Forum des Jeunes Mathématiciennes et Mathématiciens (December 2021, Besançon).
 Stochastic Geometry Days 2021 (November 2021, Dunkerque, France).
 Congrès des Jeunes Chercheurs en Mathématiques Appliquées (October 2021, Palaiseau).
 Journées des Probabilités 2021 (June 2021, Guidel).
 PHD students seminar of Institut Mathématique de Toulouse (December 2021, Toulouse).
 Working group in probability, ergodic theory and dynamical systems (November 2021, Rouen).
 S. Khaniha:
 Stochastic Geometry Days, Dunkerque, November 15th19th, 2021. Title: Dynamical systems on the null recurrent Doeblin Graph.
 Forum des Jeunes mathématicien.ne.s (probabilités, statistique et applications), December 8th10th, 2021, Title: The Doeblin Graph of a null recurrent MC.
 A. Khezeli:
 22nd Applied Stochastic Processes Workshop, Alzahra University, Iran, February 2021. Title: Invariant Circle Packings.
 8th Conference of Frontiers in Mathematical Sciences, IPM, Iran, April 2021. Title: The GromovHausdorff metric and its generalizations.
 DYOGENE reading group, minicourse on unimodular dimensions (3 lectures) April 2021, Title: Unimodular Minkowski and Hausdorff dimensions.
 Combinatorics & Computing Seminars, IPM, Iran, June 2021, Title: Invariant Circle Packing
 DYOGENE seminar, October 2021, title: Stationary and PointStationary Circle Packings.
 M.Lelarge:
 invited talk at Toutelia 2021: statistical physics, probability and AI indico.math.cnrs.fr/event/6757/
 P. Popineau
 Presentation at IEEE Globecom 2021, Madrid 43, Dec. 2021.
 Séminaire du LINCS : On Velocitybased Association Policies for Multitier 5G Wireless Networks  20/10/2021.
10.2 Teaching  Supervision  Juries
10.2.1 Teaching
 Licence: B. Błaszczyszyn (Cours) Théorie de l'information et du codage 24 heqTD, L3, ENS Paris.
 Licence: A. Busic (Cours) and K. Scaman (Cours + TD) Structures et algorithmes aléatoires 60 heqTD, L3, ENS Paris.
 Licence: L. Massoulié (Cours) Social and Communication networks 60 heqTD, L3, l'X.
 Master: A. Busic (Cours) and M. Even (TD) Modèles et algorithmes de réseaux 60 heqTD, M1, ENS Paris.
 Master: A. Busic and L. Massoulié (Cours) Fondements de la modélisation des réseaux 37.5 heqTD, M2 MPRI.
 Master: M. Lelarge (Cours) Deep Learning Do it Yourself, M1, ENS Paris, X, XHEC. mlelarge.github.io/dataflowrweb/
 Master: L. Massoulié (Cours) Inference in large random graphs, M2 Université d'Orsay.
 Master: R. Gambelin, Chargé de TD "Méthodes quantitatives", L2 Eco & Gestion, l'université Paris Descartes;
 Préparation à l'agrégation: S. Coste (Cours+TD) Cours de statistiques, agrégation de sciences sociales, ENS ParisSaclay.
 MPSI: S. Coste (Interrogations) maths.
 Prépa: R. Gambelin, Colleur, ECG 1ère année au lycée Saint Louis de Gonzague (Franklin).
 Lectures on unimodularity by F. Baccelli at the Stochastic Geometry days, Dunkerque, Nov. 2021.
10.2.2 Supervision
PhD defended:
 Arnaud Cadas, Stochastic matching models and their applications to demandsupply balancing, Nov 2021, supervised by A. Bušić; www.theses.fr/s188918
 Hadrien Hendrikx, Accelerated methods for distributed optimization 31, Sep 2021, supervised by Laurent Massoulié; www.theses.fr/s209409
 Ludovic Stephan, Community detection by the block stochastic model Oct 2021, supervised by Laurent Massoulié;
PhD in progress:
 Claire Bizon Monroc, since Nov 2021, supervised by A. Bušić
 Matthieu Blanke since Sep 2021, supervised by M. Lelarge
 Antoine Brochard since Sep 2018, supervised by B. Błaszczyszyn
 Bharath Roy Choudhury since Sep 2019, supervised by F. Baccelli
 Eric Daoud since Sep 2021, supervised by M. Lelarge
 Michel Davydov, since 2020, supervised by F. Baccelli
 Bastien Dubail since Sep 2020, supervised by Laurent Massoulié
 Mathieu Even since Sep 2021, supervised by Laurent Massoulié
 Roman Gambelin since Sep 2020, supervised by B. Blaszczyszyn
 Cedric Gerbelot since Sep 2021, supervised by M. Lelarge
 Sayeh Khaniha, since Feb 2019, supervised by F. Baccelli
 Thomas Le Corre since Nov 2021, supervised by Bušić
 Maxime Leiber since Feb 2020, supervised by Laurent Massoulié
 Pierre Popineau, since Sep 2019, supervised by F. Baccelli
 David Robin, since Oct 2021, supervised by Laurent Massoulié
 Sébastien Samain since Nov 2016, supervised by A. Bušić
 Ilia Shilov since Sep 2019, supervised by A. Bušić
 Amaury Triboulin since Sep 2021, supervised by M. Lelarge
 Lucas Weber since Oct 2021, supervised by A. Bušić
10.2.3 Juries
 Marc Lelarge: jury CRCN INRIA Paris; jury prof. UFR Informatique  Université de Paris; member of the PhD committee for Louis Martin; Reviewer of the PhD thesis for Amine Echraib.
 B. Błaszczyszyn: jury M2 defense Sofia Robert.
 F. Baccelli: member of the thesis committee of D. Anada, INRIA Lyon, October 2021.
10.3 Popularization
10.3.1 Internal or external Inria responsibilities
Les réseaux de communications à l'INRIA. Rapport établi par F. Baccelli, I. Chrisment, J.M. Gorce et P. Mussi à la demande de B. Sportisse. Novembre 2021.
10.3.2 Articles and contents
Rapport de l'académie des sciences sur les réseaux de communications du futur, établi par S. Abiteboul, D. Andler, F. Baccelli, C. Bréchignac, G. Berry, S. Candel, M. Fink et E. Moulines, juillet 2021.
11 Scientific production
11.1 Publications of the year
International journals
 1 articleReplicameanfield limits of fragmentationinteractionaggregation processes.Journal of Applied ProbabilityJanuary 2022, 122
 2 articleUnimodular Hausdorff and Minkowski dimensions.Electronic Journal of Probability26noneJanuary 2021
 3 articleOn point processes defined by angular conditions on Delaunay neighbors in the Poisson–Voronoi Tessellation.Journal of Applied Probability584December 2021, 952965
 4 articleA computational framework for evaluating the role of mobility on the propagation of epidemics on point processes.Journal of Mathematical Biology8412January 2022, 4
 5 articleThe PairReplicaMeanField Limit for Intensitybased Neural Networks.SIAM Journal on Applied Dynamical Systems201January 2021, 165207
 6 articleOn the steady state of continuoustime stochastic opinion dynamics with powerlaw confidence.Journal of Applied Probability583September 2021, 746772
 7 articleOptimal stationary markings.Stochastic Processes and their Applications138August 2021, 153185
 8 articleModeling and Analysis of Vehicle Safety Message Broadcast in Cellular Networks.IEEE Transactions on Wireless Communications207July 2021, 40874099
 9 articleEmergence of extended states at zero in the spectrum of sparse random graphs.Annals of Probability494May 2021
 10 articleContinuum LineofSight Percolation on PoissonVoronoi Tessellations.Advances in Applied Probability532June 2021, 510–536
 11 articleBeam Management in 5G: A Stochastic Geometry Analysis.IEEE Transactions on Wireless CommunicationsSeptember 2021
 12 articleA product form for the general stochastic matching model.Journal of Applied Probability582June 2021, 449468
 13 articleHow Wireless Queues Benefit from Motion: An Analysis of the Continuum Between Zero and Infinite Mobility.IEEE Transactions on Wireless Communications2012December 2021, 81498162
 14 articleMulticlass Energy Packet Networks with finite capacity energy queues.Performance Evaluation152December 2021, 102228
 15 articlePrivacy Impact on Generalized Nash Equilibrium in PeertoPeer Electricity Market.Operations Research Letters49September 2021, 759766
 16 articleMaximum UtilityAware Capacity Partitioning in Cooperative Computing.IEEE Communications Letters2510October 2021, 33603364
International peerreviewed conferences
 17 inproceedingsMean field analysis for bike and ebike sharing systems.MAMA  ACM SIGMETRICS 2021 Workshop on MAthematical performance Modeling and Analysis492Beijing, ChinaJanuary 2022, 1214
 18 inproceedingsExpressive Power of Invariant and Equivariant Graph Neural Networks.ICLR 2021  International Conference on Learning RepresentationsVirtual, Unknown RegionMay 2021
 19 inproceedingsKullbackLeiblerQuadratic Optimal Control in a Stochastic Environment.Proceedings of The 60th IEEE conference on Decision and ControlCDC 2021  60th IEEE conference on Decision and ControlAustin (online), United States2021
 20 inproceedingsControl oriented modeling of TCLs.Proceedings of the 2021 American Control Conference (ACC)2021 American Control Conference (ACC)New Orleans, United States2021
 21 inproceedingsFast Stochastic Bregman Gradient Methods: Sharp Analysis and Variance Reduction.ICML 2021 38th International Conference on Machine Learning139Proceedings of the 38th International Conference on Machine Learningvirtual, United StatesJuly 2021, 28152825
 22 inproceedingsA Continuized View on Nesterov Acceleration for Stochastic Gradient Descent and Randomized Gossip.Advances in Neural Information Processing Systems 34NeurIPS 2021  35th Conference on Neural Information Processing SystemsSydney (virtual), AustraliaMorgan Kaufmann PublishersDecember 2021, 132
 23 inproceedingsConcentration of NonIsotropic Random Tensors with Applications to Learning and Empirical Risk Minimization.Proceedings of Machine Learning ReasearchConference on Learning Theory, 2021Boulder, United StatesPMLRAugust 2021
 24 inproceedingsImpossibility of Partial Recovery in the Graph Alignment Problem.COLT 2021  34th Annual Conference on Learning TheoryBoulder / Virtual, United StatesAugust 2021
 25 inproceedingsSharp threshold for alignment of graph databases with Gaussian weights.MSML21 (Mathematical and Scientific Machine Learning)Lausanne, SwitzerlandAugust 2021
 26 inproceedingsA FineGrained Analysis of Radar Detection in Vehicular Networks.IEEE GLOBECOM 2021  IEEE Global Communications ConferenceMadrid, SpainDecember 2021
 27 inproceedingsWhere to Deploy Reconfigurable Intelligent Surfaces in the Presence of Blockages?2021 IEEE 32nd Annual International Symposium on Personal, Indoor and Mobile Radio Communications (PIMRC)Helsinki, FranceIEEESeptember 2021, 14191424
 28 inproceedingsA Generalized Nash Equilibrium analysis of the interaction between a peertopeer financial market and the distribution grid.IEEE Smart Grid Comm 2021  IEEE International Conference on Communications, Control, and Computing Technologies for Smart GridsAachen, GermanyOctober 2021
Scientific book chapters
 29 inbookFundamental Design Principles for Reinforcement Learning Algorithms.325Handbook of Reinforcement Learning and ControlStudies in Systems, Decision and ControlSpringer International PublishingJune 2021, 75137
 30 inbookRandomised Geographic Caching and its Applications in Wireless Networks.Edge Caching for Mobile NetworksThe Institution of Engineering and Technology2021, 371404
Doctoral dissertations and habilitation theses
 31 thesisAccelerated Methods for Distributed Optimization.PSLSeptember 2021
Reports & preprints
 32 miscOn Point Processes Defined by Angular Conditions on Delaunay Neighbors in the PoissonVoronoi Tessellation.January 2021
 33 miscEnergy and Delay TradeOffs of EndtoEnd Vehicular Communications using a Hyperfractal Urban Modelling.January 2022
 34 miscMean field analysis of stochastic networks with reservation.January 2022
 35 miscA unified framework for coordination of thermostatically controlled loads.August 2021
 36 miscStochastic dynamic matching: A mixed graphtheory and linearalgebra approach.January 2022
 37 miscSparse matrices: convergence of the characteristic polynomial seen from infinity.June 2021
 38 miscA simpler spectral approach for clustering in directed networks.February 2021
 39 miscAccelerating Abelian Random Walks with Hyperbolic Dynamics.June 2021
 40 reportStochastic Modelling of FreeFloating CarSharing Systems.CIRRELTFebruary 2021
 41 miscAsynchrony and Acceleration in Gossip Algorithms.February 2021
 42 miscMean field analysis of an incentive algorithm for a closed stochastic network.January 2022
 43 miscOn Velocitybased Association Policies for Multitier 5G Wireless Networks.September 2021
 44 miscNonbacktracking spectra of weighted inhomogeneous random graphs.February 2021