Publications

2025

• From random matrices to systems of particles in interaction
  
  • Pesce Valentin
  , 2025. The goal of these expository notes is to give an introduction to random matrices for non-specialist of this topic focusing on the link between random matrices and systems of particles in interaction. We first recall some general results about the random matrix theory that create a link between random matrices and systems of particles through the knowledge of the law of the eigenvalues of certain random matrices models. We next focus on a continuous in time approach of random matrices called the Dyson Brownian motion. We detail some general methods to study the existence of system of particles in singular interaction and the existence of a mean field limit for these systems of particles. Finally, we present the main result of large deviations when studying the eigenvalues of random matrices. This method is based on the fact that the eigenvalues of certain models of random matrices can be viewed as log gases in dimension 1 or 2.
• A holographic global uniqueness in passive imaging
  
  • Novikov Roman
  Journal de l'École polytechnique — Mathématiques, École polytechnique, 2025, 12, pp.1069-1081. We consider a radiation solution $\psi$ for the Helmholtz equation in an exterior region in $\mathbb R^3$. We show that the restriction of $\psi$ to any ray $L$ in the exterior region is uniquely determined by its imaginary part Im $\psi$ on an interval of this ray. As a corollary, the restriction of $\psi$ to any plane $X$ in the exterior region is uniquely determined by Im $\psi$ on an open domain in this plane. These results have holographic prototypes in the recent work Novikov (2024, Proc. Steklov Inst. Math. 325, 218-223). In particular, these and known results imply a holographic type global uniqueness in passive imaging and for the Gelfand-Krein-Levitan inverse problem (from boundary values of the spectral measure in the whole space) in the monochromatic case. Some other surfaces for measurements instead of the planes $X$ are also considered. (10.5802/jep.306)
  DOI : 10.5802/jep.306
• Integrating Aggregated Electric Vehicle Flexibilities in Unit Commitment Models using Submodular Optimization
  
  • Arvis Hélène
  • Beaude Olivier
  • Gast Nicolas
  • Gaubert Stéphane
  • Gaujal Bruno
  , 2025. <div><p>The Unit Commitment (UC) problem consists in controlling a large fleet of heterogeneous electricity production units in order to minimize the total production cost while satisfying consumer demand. Electric Vehicles (EVs) are used as a source of flexibility and are often aggregated for problem tractability. We develop a new approach to integrate EV flexibilities in the UC problem and exploit the generalized polymatroid structure of aggregated flexibilities of a large population of users to develop an exact optimization algorithm, combining a cutting-plane approach and submodular optimization. We show in particular that the UC can be solved exactly in a time which scales linearly, up to a logarithmic factor, in the number of EV users when each production unit is subject to convex constraints. We illustrate our approach by solving a real instance of a long-term UC problem, combining open-source data of the European grid (European Resource Adequacy Assessment project) and data originating from a survey of user behavior of the French EV fleet.</p></div>
• Tensor rectifiable G-flat chains
  
  • Goldman Michael
  • Merlet Benoît
  Transactions of the American Mathematical Society, American Mathematical Society, 2025. A rigidity result for normal rectifiable $k$-chains in $\mathbb{R}^n$ with coefficients in an Abelian normed group is established. Given some decompositions $k=k_1+k_2$, $n=n_1+n_2$ and some rectifiable $k$-chain $A$ in $\mathbb{R}^n$, we consider the properties: (1) The tangent planes to $\mu_A$ split as $T_x\mu_A=L^1(x)\times L^2(x)$ for some $k_1$-plane $L^1(x)\subset\mathbb{R}^{n_1}$ and some $k_2$-plane $L^2(x)\subset\mathbb{R}^{n_2}$. (2) $A=A_{\vert\Sigma^1\times\Sigma^2}$ for some sets $\Sigma^1\subset\mathbb{R}^{n_1}$, $\Sigma^2\subset\mathbb{R}^{n_2}$ such that $\Sigma^1$ is $k_1$-rectifiable and $\Sigma^2$ is $k_2$-rectifiable (we say that $A$ is $(k_1,k_2)$-rectifiable). The main result is that for normal chains, (1) implies (2), the converse is immediate. In the proof we introduce the new groups of tensor flat chains (or $(k_1,k_2)$-chains) in $\mathbb{R}^{n_1}\times\mathbb{R}^{n_2}$ which generalize Fleming's $G$-flat chains. The other main tool is White's rectifiable slices theorem. We show that on the one hand any normal rectifiable chain satisfying~(1) identifies with a normal rectifiable $(k_1,k_2)$-chain and that on the other hand any normal rectifiable $(k_1,k_2)$-chain is $(k_1,k_2)$-rectifiable. (10.1090/tran/9392)
  DOI : 10.1090/tran/9392
• Dynamics of a kinetic model describing protein exchanges in a cell population
  
  • Magal Pierre
  • Raoul Gaël
  Journal of Mathematical Biology, Springer, 2025, 91 (6), pp.76. We consider a cell population structured by a positive real number, describing the number of P-glycoproteins carried by the cell. We are interested in the effect of those proteins on the growth of the population: those proteins are indeed involve in the resistance of cancer cells to chemotherapy drugs. To describe this dynamics, we introduce a kinetic model. We then introduce a rigorous hydrodynamic limit, showing that if the exchanges are frequent, then the dynamics of the model can be described by a system of two coupled differential equations. Finally, we also show that the kinetic model converges to a unique limit in large times. The main idea of this analysis is to use Wasserstein distance estimates to describe the effect of the kinetic operator, combined to more classical estimates on the macroscopic quantities. (10.1007/s00285-025-02295-w)
  DOI : 10.1007/s00285-025-02295-w
• Parameters estimation of a Threshold CKLS process from continuous and discrete observations
  
  • Mazzonetto Sara
  • Nieto Benoît
  Scandinavian Journal of Statistics, Wiley, 2025, 52 (4), pp.1670-1707. We consider a continuous time process that is self-exciting and ergodic, called threshold Chan–Karolyi–Longstaff–Sanders (CKLS) process. This process is a generalization of various models in econometrics, such as Vasicek model, Cox-Ingersoll-Ross, and Black-Scholes, allowing for the presence of several thresholds which determine changes in the dynamics. We study the asymptotic behavior of maximum-likelihood and quasi-maximum-likelihood estimators of the drift parameters in the case of continuous time and discrete time observations. We show that for high frequency observations and infinite horizon the estimators satisfy the same asymptotic normality property as in the case of continuous time observations. We also discuss diffusion coefficient estimation. Finally, we apply our estimators to simulated and real data to motivate considering (multiple) thresholds. (10.1111/sjos.70005)
  DOI : 10.1111/sjos.70005
• Inverse scattering for the multipoint potentials of Bethe-Peierls-Thomas-Fermi type
  
  • Kuo Pei-Cheng
  • Novikov Roman
  Inverse Problems, IOP Publishing, 2025, 41 (6), pp.065021. <div><p>We consider the Schrödinger equation with a multipoint potential of the Bethe-Peierls-Thomas-Fermi type. We show that such a potential in dimension d = 2 or d = 3 is uniquely determined by its scattering amplitude at a fixed positive energy. Moreover, we show that there is no non-zero potential of this type with zero scattering amplitude at a fixed positive energy and a fixed incident direction. Nevertheless, we also show that a multipoint potential of this type is not uniquely determined by its scattering amplitude at a positive energy E and a fixed incident direction. Our proofs also contribute to the theory of inverse source problem for the Helmholtz equation with multipoint source.</p></div> (10.1088/1361-6420/ade282)
  DOI : 10.1088/1361-6420/ade282
• PSWF-Radon approach to reconstruction from band-limited Hankel transform
  
  • Goncharov Fedor
  • Isaev Mikhail
  • Novikov Roman
  • Zaytsev Rodion
  Applied and Computational Harmonic Analysis, Elsevier, 2025. We give new formulas for reconstructions from band-limited Hankel transform of integer or half-integer order. Our formulas rely on the PSWF-Radon approach to super-resolution in multidimensional Fourier analysis. This approach consists of combining the theory of classical one-dimensional prolate spheroidal wave functions with the Radon transform theory. We also use the relation between Fourier and Hankel transforms and Cormack-type inversion of the Radon transform. Finally, we investigate numerically the capabilities of our approach to super-resolution for band-limited Hankel inversion in relation to varying levels of noise.
• Byzantine-Robust Gossip: Insights from a Dual Approach
  
  • Gaucher Renaud
  • Dieuleveut Aymeric
  • Hendrikx Hadrien
  , 2025. Distributed learning has many computational benefits but is vulnerable to attacks from a subset of devices transmitting incorrect information. This paper investigates Byzantine resilient algorithms in a decentralized setting, where devices communicate directly in a peer-to-peer manner within a communication network. We leverage the so-called dual approach for decentralized optimization and propose a Byzantine-robust algorithm. We provide convergence guarantees in the average consensus subcase, discuss the potential of the dual approach beyond this subcase, and re-interpret existing algorithms using the dual framework. Lastly, we experimentally show the soundness of our method.
• Taking Advantage of Multiple Scattering for Optical Reflection Tomography
  
  • Wasik Thomas
  • Barolle Victor
  • Aubry Alexandre
  • Garnier Josselin
  , 2025. <div><p>Optical Diffraction Tomography (ODT) is a powerful non-invasive imaging technique widely used in biological and medical applications. While significant progress has been made in transmission configuration, reflection ODT remains challenging due to the ill-posed nature of the inverse problem. We present a novel optimization algorithm for 3D refractive index (RI) reconstruction in reflection-mode microscopy. Our method takes advantage of the multiply-scattered waves that are reflected by uncontrolled background structures and that illuminate the foreground RI from behind. It tackles the ill-posed nature of the problem using weighted time loss, positivity constraints and Total Variation regularization. We have validated our method with data generated by detailed 2D and 3D simulations, demonstrating its performance under weak multiple scattering conditions and with simplified forward models used in the optimization routine for computational efficiency.</p><p>In addition, we highlight the need for multi-wavelength analysis and the use of regularization to ensure the reconstruction of the low spatial frequencies of the foreground RI.</p></div>
• Optimisation of space-time periodic eigenvalues
  
  • Bogosel Beniamin
  • Mazari Idriss
  • Nadin Grégoire
  Annali della Scuola Normale Superiore di Pisa, Classe di Scienze, Scuola Normale Superiore, 2025. <div><p>The goal of this paper is to provide a qualitative analysis of the optimisation of space-time periodic principal eigenvalues. Namely, considering a fixed time horizon T and the d-dimensional torus T d , let, for any m ∈ L ∞ ((0, T ) × T d ), λ(m) be the principal eigenvalue of the operator ∂t -∆ -m endowed with (time-space) periodic boundary conditions. The main question we set out to answer is the following: how to choose m so as to minimise λ(m)? This question stems from population dynamics. We prove that in several cases it is always beneficial to rearrange m with respect to time in a symmetric way, which is the first comparison result for the rearrangement in time of parabolic equations. Furthermore, we investigate the validity (or lack thereof) of Talenti inequalities for the rearrangement in time of parabolic equations. The numerical simulations which illustrate our results were obtained by developing a framework within which it is possible to optimise criteria with respect to functions having a prescribed rearrangement (or distribution function).</p></div> (10.2422/2036-2145.202501_012)
  DOI : 10.2422/2036-2145.202501_012
• Efficient treatment of the model error in the calibration of computer codes: the Complete Maximum a Posteriori method
  
  • Kahol Omar
  • Congedo Pietro Marco
  • Le Maitre Olivier
  • Denimal Goy Enora
  International Journal for Uncertainty Quantification, Begell House Publishers, 2025. <div><p>Computer models are widely used for the prediction of complex physical phenomena. Based on observations of these physical phenomena, it is possible to calibrate the model parameters. In most cases, such computer models are mis-specified, and the calibration process must be improved by including a model error term. The model error hyperparameters are, however, rarely learned jointly with the model parameters to reduce the dimensionality of the problem. Sequential and non-sequential approaches have been introduced to estimate the hyperparameters. The former, such as the Kennedy and O'Hagan (KOH) framework, estimates the model error hyperparameters before calibrating the model parameters. The latter, such as the Full Maximum a Posteriori (FMP), introduces a functional dependence between the model parameters and the model error hyperparameters. Despite being more reliable in some cases (bimodality e.g.), the FMP method still fails to estimate correctly the posterior distribution shape. This work proposes a new methodology for treating the model error term in computer code calibration. It builds upon the KOH and FMP framework. Called the Complete Maximum a Posteriori (CMP) method, it provides a closed-form expression for the marginalization integral over the model error hyperparameters, significantly reducing the dimensionality of the calibration problem. Such expression re- lies on a set of assumptions that are more general and less stringent than the ones usually employed. The CMP method is applied to four examples of increasing complexity, from elementary to real fluid dynamics problems, including or not bimodality. Compared to the true reference solution and unlike the KOH and FMP, the CMP method correctly captures the shape of the posterior distribution, including all modes and their weights. Moreover, it provides an accurate estimate of the distribution tails</p></div> (10.1615/Int.J.UncertaintyQuantification.2025056317)
  DOI : 10.1615/Int.J.UncertaintyQuantification.2025056317
• Long time behavior of a degenerate stochastic system modeling the response of a population face to environmental impacts
  
  • Collet Pierre
  • Ecotière Claire
  • Méléard Sylvie
  Electronic Communications in Probability, Institute of Mathematical Statistics (IMS), 2025, 30 (none). We study the asymptotics of a two-dimensional stochastic differential system with a degenerate diffusion matrix. This system describes the dynamics of a population where individuals contribute to the degradation of their environment through two different behaviors. We exploit the almost one-dimensional form of the dynamical system to compute explicitly the Freidlin-Wentzell action functional. That allows to give conditions under which the small noise regime of the invariant measure is concentrated around the equilibrium of the dynamical system having the smallest diffusion coefficient. (10.1214/24-ECP650)
  DOI : 10.1214/24-ECP650
• Adaptive Destruction Processes for Diffusion Samplers
  
  • Gritsaev Timofei
  • Morozov Nikita
  • Tamogashev Kirill
  • Tiapkin Daniil
  • Samsonov Sergey
  • Naumov Alexey
  • Vetrov Dmitry
  • Malkin Nikolay
  , 2025. This paper explores the challenges and benefits of a trainable destruction process in diffusion samplers -- diffusion-based generative models trained to sample an unnormalised density without access to data samples. Contrary to the majority of work that views diffusion samplers as approximations to an underlying continuous-time model, we view diffusion models as discrete-time policies trained to produce samples in very few generation steps. We propose to trade some of the elegance of the underlying theory for flexibility in the definition of the generative and destruction policies. In particular, we decouple the generation and destruction variances, enabling both transition kernels to be learned as unconstrained Gaussian densities. We show that, when the number of steps is limited, training both generation and destruction processes results in faster convergence and improved sampling quality on various benchmarks. Through a robust ablation study, we investigate the design choices necessary to facilitate stable training. Finally, we show the scalability of our approach through experiments on GAN latent space sampling for conditional image generation. (10.48550/arXiv.2506.01541)
  DOI : 10.48550/arXiv.2506.01541
• Crediting football players for creating dangerous actions in an unbiased way: the generation of threat (GoT) indices
  
  • Baouan Ali
  • Coustou Sebastien
  • Lacome Mathieu
  • Pulido Sergio
  • Rosenbaum Mathieu
  Journal of Quantitative Analysis in Sports, De Gruyter, 2025. We introduce an innovative methodology to identify football players at the origin of threatening actions in a team. In our framework, a threat is defined as entering the opposing team's danger area. We investigate the timing of threat events and ball touches of players, and capture their correlation using Hawkes processes. Our model-based approach allows us to evaluate a player's ability to create danger both directly and through interactions with teammates. We define a new index, called Generation of Threat (GoT), that measures in an unbiased way the contribution of a player to threat generation. For illustration, we present a detailed analysis of Chelsea's 2016-2017 season, with a standout performance from Eden Hazard. We are able to credit each player for his involvement in danger creation and determine the main circuits leading to threat. In the same spirit, we investigate the danger generation process of Stade Rennais in the 2021-2022 season. Furthermore, we establish a comprehensive ranking of Ligue 1 players based on their generated threat in the 2021-2022 season. Our analysis reveals surprising results, with players such as Jason Berthomier, Moses Simon and Frederic Guilbert among the top performers in the GoT rankings. We also present a ranking of Ligue 1 central defenders in terms of generation of threat and confirm the great performance of some center-back pairs, such as Nayef Aguerd and Warmed Omari.
• Partial regularity for optimal transport with $p$-cost away from fixed points
  
  • Goldman Michael
  • Koch Lukas
  Proceedings of the American Mathematical Society, American Mathematical Society, 2025. We consider maps $T$ solving the optimal transport problem with a cost $c(x-y)$ modeled on the $p$-cost. For Hölder continuous marginals, we prove a $C^{1,\alpha}$-partial regularity result for $T$ in the set $\{\lvert T(x)-x\rvert&gt;0\}$.
• Signature volatility models: pricing and hedging with Fourier
  
  • Abi Jaber Eduardo
  • Gérard Louis-Amand
  SIAM Journal on Financial Mathematics, Society for Industrial and Applied Mathematics, 2025, 16 (2). We consider a stochastic volatility model where the dynamics of the volatility are given by a possibly infinite linear combination of the elements of the time extended signature of a Brownian motion. First, we show that the model is remarkably universal, as it includes, but is not limited to, the celebrated Stein-Stein, Bergomi, and Heston models, together with some path-dependent variants. Second, we derive the joint characteristic functional of the log-price and integrated variance provided that some infinitedimensional extended tensor algebra valued Riccati equation admits a solution. This allows us to price and (quadratically) hedge certain European and path-dependent options using Fourier inversion techniques. We highlight the efficiency and accuracy of these Fourier techniques in a comprehensive numerical study. (10.1137/24M1636952)
  DOI : 10.1137/24M1636952
• Probing the speckle to estimate the effective speed of sound, a first step towards quantitative ultrasound imaging
  
  • Garnier Josselin
  • Giovangigli Laure
  • Goepfert Quentin
  • Millien Pierre
  Inverse Problems and Imaging, AIMS American Institute of Mathematical Sciences, 2025 (1930-8337). <div><p>In this paper, we present a mathematical model and analysis for a new experimental method [Bureau and al., arXiv:2409.13901, 2024] for effective sound velocity estimation in medical ultrasound imaging. We perform a detailed analysis of the point spread function of a medical ultrasound imaging system when there is a mismatch between the effective sound speed in the medium and the one used in the backpropagation imaging functional. Based on this analysis, an estimator for the speed of sound error is introduced. Using recent results on stochastic homogenization of the Helmholtz equation, we provide a representation formula for the field scattered by a random multi-scale medium (whose acoustic behavior is similar to a biological tissue) in the time-harmonic regime. We then prove that statistical moments of the imaging function can be accessed from data collected with only one realization of the medium. We show that it is possible to locally extract the point spread function from an image constituted only of speckle and build an estimator for the effective sound velocity in the micro-structured medium. Some numerical illustrations are presented at the end of the paper.</p></div> (10.3934/ipi.2026001)
  DOI : 10.3934/ipi.2026001
• Learning homogenized hyperelastic behavior for topology optimization of lattice structures
  
  • Ribeiro Nogueira Breno
  • Allaire Grégoire
  , 2025. (10.5281/zenodo.14900138)
  DOI : 10.5281/zenodo.14900138
• Objective assessment of cardiac function using patient-specific biophysical modeling based on cardiovascular MRI combined with catheterization
  
  • Gusseva Maria
  • Castellanos Daniel Alexander
  • Veeram Reddy Surendranath
  • Hussain Tarique
  • Chapelle Dominique
  • Chabiniok Radomír
  AJP - Heart and Circulatory Physiology, American Physiological Society, 2025, 329 (5), pp.H118-H1191. Synthesizing multi-modality data, such as cardiovascular magnetic resonance imaging (MRI) combined with catheterization, into a single framework is challenging. Different acquisition systems are subjected to different measurement errors. Coupling clinical data with biomechanical models can assist in clinical data processing (e.g., model-based filtering of measurement noise) and quantify myocardial mechanics via metrics not readily available in the data, such as myocardial contractility. In this work we use a biomechanical modeling with the aim 1) to quantitatively compare model- and data-derived signals, and 2) to explore the potential of model-derived myocardial contractility and distal resistance of the circulation (Rd) to robustly quantify cardiovascular physiology. We used 51 ventricular catheterization pressure and cine MRI volume datasets from patients with single-ventricle physiology and left and right ventricles of patients with repaired tetralogy of Fallot. Ventricular time-varying elastance (TVE) metrics and linear regression were used to quantify the relationship between the maximum value of TVE (Emax) and maximum time derivative of ventricular pressure (max(dP/dt)) in data- and model-derived pressure and volume signals at p&lt;0.05. Pearson’s correlations were used to compare model-derived contractility and data-derived Emax and max(dP/dt), and model-derived Rd and data-derived vascular resistance. All data and model-derived linear regressions were significant (p&lt;0.05). Model-derived max(dP/dt) vs. data-derived Emax produced higher R2 than data-derived max(dP/dt) vs. data-derived Emax. Correlations demonstrated significant relationships between most data- and model-derived metrics. This work revealed the clinical value of biomechanical modeling to assist in clinical data processing by providing high-quality pressure and volume signals, and to quantify cardiovascular pathophysiology. (10.1152/ajpheart.00232.2025)
  DOI : 10.1152/ajpheart.00232.2025
• Finite elements for Wasserstein $W_p$ gradient flows
  
  • Cancès Clément
  • Matthes Daniel
  • Nabet Flore
  • Rott Eva-Maria
  ESAIM: Mathematical Modelling and Numerical Analysis, Société de Mathématiques Appliquées et Industrielles (SMAI) / EDP, 2025, 59 (3), pp.1565-1600. Wasserstein $\bbW_p$ gradient flows for nonlinear integral functionals of the density yield degenerate parabolic equations involving diffusion operators of $q$-Laplacian type, with $q$ being $p$'s conjugate exponent. We propose a finite element scheme building on conformal $\mathbb{P}_1$ Lagrange elements with mass lumping and a backward Euler time discretization strategy. Our scheme preserves mass and positivity while energy decays in time. Building on the theory of gradient flows in metric spaces, we further prove convergence towards a weak solution of the PDE that satisfies the energy dissipation equality. The analytical results are illustrated by numerical simulations. (10.1051/m2an/2025035)
  DOI : 10.1051/m2an/2025035
• Volume growth of Funk geometry and the flags of polytopes
  
  • Faifman Dmitry
  • Vernicos Constantin
  • Walsh Cormac
  Geometry and Topology, Mathematical Sciences Publishers, 2025, 29 (7), pp.3773-3811. We consider the Holmes-Thompson volume of balls in the Funk geometry on the interior of a convex domain. We conjecture that for a fixed radius, this volume is minimized when the domain is a simplex and the ball is centered at the barycenter, or in the centrally-symmetric case, when the domain is a Hanner polytope. This interpolates between Mahler's conjecture and Kalai's flag conjecture. We verify this conjecture for unconditional domains. For polytopal Funk geometries, we study the asymptotics of the volume of balls of large radius, and compute the two highest-order terms. The highest depends only on the combinatorics, namely on the number of flags. The second highest depends also on the geometry, and thus serves as a geometric analogue of the centro-affine area for polytopes. We then show that for any polytope, the second highest coefficient is minimized by a unique choice of center point, extending the notion of Santaló point. Finally, we show that, in dimension two, this coefficient, with respect to the minimal center point, is uniquely maximized by affine images of the regular polygon. (10.2140/gt.2025.29.3773)
  DOI : 10.2140/gt.2025.29.3773
• A holographic uniqueness theorem for the two-dimensional Helmholtz equation
  
  • Nair Arjun
  • Novikov Roman
  The Journal of Geometric Analysis, Springer, 2025, 35 (4), pp.123. We consider a plane wave, a radiation solution, and the sum of these solutions (total solution) for the Helmholtz equation in an exterior region in $\mathbb R^2$. We consider a straight line in this region, such that the direction of propagation of the plane wave is not parallel to this line. We show that the radiation solution in the exterior region is uniquely determined by the intensity of the total solution on an interval of this line. In particular, this result solves one of the old mathematical questions of holography in its two-dimensional setting. Our proofs also contribute to the theory of the Karp expansion of radiation solutions in two dimensions. (10.1007/s12220-025-01949-x)
  DOI : 10.1007/s12220-025-01949-x
• Provable non-accelerations of the heavy-ball method
  
  • Goujaud Baptiste
  • Taylor Adrien
  • Dieuleveut Aymeric
  Mathematical Programming, Springer Verlag, 2025. In this work, we show that the heavy-ball ($\HB$) method provably does not reach an accelerated convergence rate on smooth strongly convex problems. More specifically, we show that for any condition number and any choice of algorithmic parameters, either the worst-case convergence rate of $\HB$ on the class of $L$-smooth and $μ$-strongly convex \textit{quadratic} functions is not accelerated (that is, slower than $1 - \mathcal{O}(κ)$), or there exists an $L$-smooth $μ$-strongly convex function and an initialization such that the method does not converge. To the best of our knowledge, this result closes a simple yet open question on one of the most used and iconic first-order optimization technique. Our approach builds on finding functions for which $\HB$ fails to converge and instead cycles over finitely many iterates. We analytically describe all parametrizations of $\HB$ that exhibit this cycling behavior on a particular cycle shape, whose choice is supported by a systematic and constructive approach to the study of cycling behaviors of first-order methods. We show the robustness of our results to perturbations of the cycle, and extend them to class of functions that also satisfy higher-order regularity conditions. (10.1007/s10107-025-02269-2)
  DOI : 10.1007/s10107-025-02269-2
• Busemann Functions in the Wasserstein Space: Existence, Closed-Forms, and Applications to Slicing
  
  • Bonet Clément
  • Cazelles Elsa
  • Drumetz Lucas
  • Courty Nicolas
  , 2025. The Busemann function has recently found much interest in a variety of geometric machine learning problems, as it naturally defines projections onto geodesic rays of Riemannian manifolds and generalizes the notion of hyperplanes. As several sources of data can be conveniently modeled as probability distributions, it is natural to study this function in the Wasserstein space, which carries a rich formal Riemannian structure induced by Optimal Transport metrics. In this work, we investigate the existence and computation of Busemann functions in Wasserstein space, which admits geodesic rays. We establish closed-form expressions in two important cases: one-dimensional distributions and Gaussian measures. These results enable explicit projection schemes for probability distributions on $\mathbb{R}$, which in turn allow us to define novel Sliced-Wasserstein distances over Gaussian mixtures and labeled datasets. We demonstrate the efficiency of those original schemes on synthetic datasets as well as transfer learning problems. (10.48550/arXiv.2510.04579)
  DOI : 10.48550/arXiv.2510.04579