Pietro Caputo
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Pietro Caputo has a total of 66 co-authors. Here are the top co-authors by number of shared publications:
- Fabio Martinelli
Università Roma Tre
16 shared publications - Fabio Lucio Toninelli
TU Wien, Vienna, Austria
15 shared publications - Matteo Quattropani
Dipartimento di Matematica e Fisica, Università di Roma Tre, Rome, Italy
8 shared publications - Alistair Sinclair
University of California, Berkeley, USA
5 shared publications Fabio Martinelli
5 shared publications- Daniel Parisi
Department of Mathematics, University of Roma Tre, Largo San Murialdo 1, 00146 Roma, Italy
4 shared publications
Latest Publications
- Open Access
Kac's Program and Relative Entropy Decay for Nonlinear Spin-Exchange Dynamics
We introduce and analyze a nonlinear exchange dynamics for Ising spin systems with arbitrary interactions. The evolution is governed by a quadratic Boltzmann-type equation that conserves the mean magnetization. Collisions are encoded through a spin-exchange kernel chosen so that the dynamics converge to the Ising model with the prescribed interaction and mean magnetization profile determined by the initial state. We prove a general convergence theorem, valid for any interaction and any transport kernel. Moreover, we show that, for sufficiently weak interactions, the system relaxes exponentially fast to equilibrium in relative entropy, with optimal decay rate independent of the initial condition. The proof relies on establishing a strong version of the Kac program from kinetic theory. In particular, we show that the associated Kac particle system satisfies a modified logarithmic Sobolev inequality with constants uniform in the number of particles. This is achieved by adapting the method of stochastic localization to the present conservative setting.
Open AccessKac's Program and Relative Entropy Decay for Nonlinear Spin-Exchange Dynamics
We introduce and analyze a nonlinear exchange dynamics for Ising spin systems with arbitrary interactions. The evolution is governed by a quadratic Boltzmann-type equation that conserves the mean magnetization. Collisions are encoded through a spin-exchange kernel chosen so that the dynamics converge to the Ising model with the prescribed interaction and mean magnetization profile determined by the initial state. We prove a general convergence theorem, valid for any interaction and any transport kernel. Moreover, we show that, for sufficiently weak interactions, the system relaxes exponentially fast to equilibrium in relative entropy, with optimal decay rate independent of the initial condition. The proof relies on establishing a strong version of the Kac program from kinetic theory. In particular, we show that the associated Kac particle system satisfies a modified logarithmic Sobolev inequality with constants uniform in the number of particles. This is achieved by adapting the method of stochastic localization to the present conservative setting.
- 11 July 2025
- Published by Springer Nature in Journal of Statistical Physics
- Vol. 192 (7) , 1-28
Relaxation Time and Topology in 1D O(N) Models
- 11 July 2025
- Published by Springer Nature in Journal of Statistical Physics
- Vol. 192 (7) , 1-28
Relaxation Time and Topology in 1D O(N) Models
- Open Access
A universal cutoff phenomenon for mean-field exchange models
We study a broad class of high-dimensional mean-field exchange models, encompassing both noisy and singular dynamics, along with their dual processes. This includes a generalized version of the averaging process as well as some non-reversible extensions of classical exchange dynamics, such as the flat Kac model. Within a unified framework, we analyze convergence to stationarity from worst-case initial data in Wasserstein distance. Our main result establishes a universal cutoff phenomenon at an explicit mixing time, with a precise window and limiting Gaussian profile. The mixing time and profile are characterized in terms of the logarithm of the size-biased redistribution random variable, thus admitting a natural entropic interpretation.
Open AccessA universal cutoff phenomenon for mean-field exchange models
We study a broad class of high-dimensional mean-field exchange models, encompassing both noisy and singular dynamics, along with their dual processes. This includes a generalized version of the averaging process as well as some non-reversible extensions of classical exchange dynamics, such as the flat Kac model. Within a unified framework, we analyze convergence to stationarity from worst-case initial data in Wasserstein distance. Our main result establishes a universal cutoff phenomenon at an explicit mixing time, with a precise window and limiting Gaussian profile. The mixing time and profile are characterized in terms of the logarithm of the size-biased redistribution random variable, thus admitting a natural entropic interpretation.
- Open Access
Factorizations of relative entropy using stochastic localization
We derive entropy factorization estimates for spin systems using the stochastic localization approach proposed by Eldan and Chen-Eldan, which, in this context, is equivalent to the renormalization group approach developed independently by Bauerschmidt, Bodineau, and Dagallier. The method provides approximate Shearer-type inequalities for the corresponding Gibbs measure at sufficiently high temperature, without restrictions on the degree of the underlying graph. For Ising systems, these are shown to hold up to the critical tree-uniqueness threshold, including polynomial bounds at the critical point, with optimal $O(\sqrt n)$ constants for the Curie-Weiss model at criticality. In turn, these estimates imply tight mixing time bounds for arbitrary block dynamics or Gibbs samplers, improving over existing results. Moreover, we establish new tensorization statements for the Shearer inequality asserting that if a system consists of weakly interacting but otherwise arbitrary components, each of which satisfies an approximate Shearer inequality, then the whole system also satisfies such an estimate.Open AccessFactorizations of relative entropy using stochastic localization
We derive entropy factorization estimates for spin systems using the stochastic localization approach proposed by Eldan and Chen-Eldan, which, in this context, is equivalent to the renormalization group approach developed independently by Bauerschmidt, Bodineau, and Dagallier. The method provides approximate Shearer-type inequalities for the corresponding Gibbs measure at sufficiently high temperature, without restrictions on the degree of the underlying graph. For Ising systems, these are shown to hold up to the critical tree-uniqueness threshold, including polynomial bounds at the critical point, with optimal $O(\sqrt n)$ constants for the Curie-Weiss model at criticality. In turn, these estimates imply tight mixing time bounds for arbitrary block dynamics or Gibbs samplers, improving over existing results. Moreover, we establish new tensorization statements for the Shearer inequality asserting that if a system consists of weakly interacting but otherwise arbitrary components, each of which satisfies an approximate Shearer inequality, then the whole system also satisfies such an estimate. - 22 January 2025
- Published by American Mathematical Society (AMS) in Transactions of the American Mathematical Society
- Vol. 378 (5) , 3551-3571
Entropy and curvature: Beyond the Peres-Tetali conjecture
We study Markov chains with non-negative sectional curvature on finite metric spaces. Neither reversibility, nor the restriction to a particular combinatorial distance is imposed. In this level of generality, we prove that a 1-step contraction in the Wasserstein distance implies a 1-step contraction in relative entropy, by the same amount. Our result substantially strengthens a recent breakthrough of the second author, and has the advantage of being applicable to arbitrary scales. This leads to a time-varying refinement of the standard Modified Log-Sobolev Inequality (MLSI), which allows us to leverage the well-acknowledged fact that curvature improves at large scales. We illustrate this principle with several applications, including birth and death chains, colored exclusion processes, permutation walks, Gibbs samplers for high-temperature spin systems, and attractive zero-range dynamics. In particular, we prove an MLSI with constant equal to the minimal rate increment for the mean-field zero-range process, thereby answering a long-standing question.
- 22 January 2025
- Published by American Mathematical Society (AMS) in Transactions of the American Mathematical Society
- Vol. 378 (5) , 3551-3571
Entropy and curvature: Beyond the Peres-Tetali conjecture
We study Markov chains with non-negative sectional curvature on finite metric spaces. Neither reversibility, nor the restriction to a particular combinatorial distance is imposed. In this level of generality, we prove that a 1-step contraction in the Wasserstein distance implies a 1-step contraction in relative entropy, by the same amount. Our result substantially strengthens a recent breakthrough of the second author, and has the advantage of being applicable to arbitrary scales. This leads to a time-varying refinement of the standard Modified Log-Sobolev Inequality (MLSI), which allows us to leverage the well-acknowledged fact that curvature improves at large scales. We illustrate this principle with several applications, including birth and death chains, colored exclusion processes, permutation walks, Gibbs samplers for high-temperature spin systems, and attractive zero-range dynamics. In particular, we prove an MLSI with constant equal to the minimal rate increment for the mean-field zero-range process, thereby answering a long-standing question.
- Open Access
- 1 January 2025
- Published by Institute of Mathematical Statistics in Electronic Journal of Probability
- Vol. 30 (none) , 1-29
Entropy contractions in Markov chains: Half-step, full-step and continuous-time
This paper considers the speed of convergence (mixing) of a finite Markov kernel P with respect to the Kullback-Leibler divergence (entropy). Given a Markov kernel one defines either a discrete-time Markov chain (with the n-step transition kernel given by the matrix power Pn) or a continuous-time Markov process (with the time-t transition kernel given by et(P−Id)). The contraction of entropy for n=1 or t=0+ are characterized by the famous functional inequalities, the strong data processing inequality (SDPI) and the modified log-Sobolev inequality (MLSI), respectively. When P=KK∗ is written as the product of a kernel and its adjoint, one could also consider the “half-step” contraction, which is the SDPI for K, while the “full-step” contraction refers to the SDPI for P. The work [18] claimed that these contraction coefficients (half-step, full-step, and continuous-time) are generally comparable, that is their ratio is bounded from above and below by two absolute constants. We disprove this and related conjectures by working out a number of different counterexamples. In particular, we construct (a) a continuous-time Markov process that contracts arbitrarily faster than its discrete-time counterpart; and (b) a kernel P such that Pm+1 contracts arbitrarily better than Pm. Hence, our main conclusion is that the four standard inequalities comparing five common notions of entropy and variance contraction are generally not improvable (even in the subclass of factorizable Markov chains). In the process of analyzing the counterexamples, we survey and sharpen the tools for bounding the contraction coefficients and characterize properties of extremizers of the respective functional inequalities. For example, we show that while the MLSI extremizers always have full support (unless the MLSI constant equals twice the spectral gap), the SDPI extremizers can have partial support. As our examples range from Bernoulli-Laplace model, random walks on graphs, to birth-death chains, the paper is also intended as a survey on computing MLSI, SDPI and other constants for these types of commonly occurring Markov chains.
- 1 January 2025
- Published by Institute of Mathematical Statistics in Electronic Journal of Probability
- Vol. 30 (none) , 1-29
Open AccessEntropy contractions in Markov chains: Half-step, full-step and continuous-time
This paper considers the speed of convergence (mixing) of a finite Markov kernel P with respect to the Kullback-Leibler divergence (entropy). Given a Markov kernel one defines either a discrete-time Markov chain (with the n-step transition kernel given by the matrix power Pn) or a continuous-time Markov process (with the time-t transition kernel given by et(P−Id)). The contraction of entropy for n=1 or t=0+ are characterized by the famous functional inequalities, the strong data processing inequality (SDPI) and the modified log-Sobolev inequality (MLSI), respectively. When P=KK∗ is written as the product of a kernel and its adjoint, one could also consider the “half-step” contraction, which is the SDPI for K, while the “full-step” contraction refers to the SDPI for P. The work [18] claimed that these contraction coefficients (half-step, full-step, and continuous-time) are generally comparable, that is their ratio is bounded from above and below by two absolute constants. We disprove this and related conjectures by working out a number of different counterexamples. In particular, we construct (a) a continuous-time Markov process that contracts arbitrarily faster than its discrete-time counterpart; and (b) a kernel P such that Pm+1 contracts arbitrarily better than Pm. Hence, our main conclusion is that the four standard inequalities comparing five common notions of entropy and variance contraction are generally not improvable (even in the subclass of factorizable Markov chains). In the process of analyzing the counterexamples, we survey and sharpen the tools for bounding the contraction coefficients and characterize properties of extremizers of the respective functional inequalities. For example, we show that while the MLSI extremizers always have full support (unless the MLSI constant equals twice the spectral gap), the SDPI extremizers can have partial support. As our examples range from Bernoulli-Laplace model, random walks on graphs, to birth-death chains, the paper is also intended as a survey on computing MLSI, SDPI and other constants for these types of commonly occurring Markov chains.
- Open Access
- 12 October 2024
- Published by Springer Nature in Communications in Mathematical Physics
- Vol. 405 (11) , 1-48
Nonlinear Dynamics for the Ising Model
We introduce and analyze a natural class of nonlinear dynamics for spin systems such as the Ising model. This class of dynamics is based on the framework of mass action kinetics, which models the evolution of systems of entities under pairwise interactions, and captures a number of important nonlinear models from various fields, including chemical reaction networks, Boltzmann’s model of an ideal gas, recombination in population genetics and genetic algorithms. In the context of spin systems, it is a natural generalization of linear dynamics based on Markov chains, such as Glauber dynamics and block dynamics, which are by now well understood. However, the inherent nonlinearity makes the dynamics much harder to analyze, and rigorous quantitative results so far are limited to processes which converge to essentially trivial stationary distributions that are product measures. In this paper we provide the first quantitative convergence analysis for natural nonlinear dynamics in a combinatorial setting where the stationary distribution contains non-trivial correlations, namely spin systems at high temperatures. We prove that nonlinear versions of both the Glauber dynamics and the block dynamics converge to the Gibbs distribution of the Ising model (with given external fields) in times $$O(n\log n)$$ and $$O(\log n)$$ respectively, where n is the size of the underlying graph (number of spins). Given the lack of general analytical methods for such nonlinear systems, our analysis is unconventional, and combines tools such as information percolation (due in the linear setting to Lubetzky and Sly), a novel coupling of the Ising model with Erdős-Rényi random graphs, and non-traditional branching processes augmented by a “fragmentation” process. Our results extend immediately to any spin system with a finite number of spins and bounded interactions.- 12 October 2024
- Published by Springer Nature in Communications in Mathematical Physics
- Vol. 405 (11) , 1-48
Open AccessNonlinear Dynamics for the Ising Model
We introduce and analyze a natural class of nonlinear dynamics for spin systems such as the Ising model. This class of dynamics is based on the framework of mass action kinetics, which models the evolution of systems of entities under pairwise interactions, and captures a number of important nonlinear models from various fields, including chemical reaction networks, Boltzmann’s model of an ideal gas, recombination in population genetics and genetic algorithms. In the context of spin systems, it is a natural generalization of linear dynamics based on Markov chains, such as Glauber dynamics and block dynamics, which are by now well understood. However, the inherent nonlinearity makes the dynamics much harder to analyze, and rigorous quantitative results so far are limited to processes which converge to essentially trivial stationary distributions that are product measures. In this paper we provide the first quantitative convergence analysis for natural nonlinear dynamics in a combinatorial setting where the stationary distribution contains non-trivial correlations, namely spin systems at high temperatures. We prove that nonlinear versions of both the Glauber dynamics and the block dynamics converge to the Gibbs distribution of the Ising model (with given external fields) in times $$O(n\log n)$$ and $$O(\log n)$$ respectively, where n is the size of the underlying graph (number of spins). Given the lack of general analytical methods for such nonlinear systems, our analysis is unconventional, and combines tools such as information percolation (due in the linear setting to Lubetzky and Sly), a novel coupling of the Ising model with Erdős-Rényi random graphs, and non-traditional branching processes augmented by a “fragmentation” process. Our results extend immediately to any spin system with a finite number of spins and bounded interactions. - Open Access
Entropy Contractions in Markov Chains: Half-Step, Full-Step and Continuous-Time
This paper considers the speed of convergence (mixing) of a finite Markov kernel $P$ with respect to the Kullback-Leibler divergence (entropy). Given a Markov kernel one defines either a discrete-time Markov chain (with the $n$-step transition kernel given by the matrix power $P^n$) or a continuous-time Markov process (with the time-$t$ transition kernel given by $e^{t(P-\mathrm{Id})}$). The contraction of entropy for $n=1$ or $t=0+$ are characterized by the famous functional inequalities, the strong data processing inequality (SDPI) and the modified log-Sobolev inequality (MLSI), respectively. When $P=KK^*$ is written as the product of a kernel and its adjoint, one could also consider the ``half-step'' contraction, which is the SDPI for $K$, while the ``full-step'' contraction refers to the SDPI for $P$. The work [DMLM03] claimed that these contraction coefficients (half-step, full-step, and continuous-time) are generally within a constant factor of each other. We disprove this and related conjectures by working out a number of different counterexamples. In particular, we construct (a) a continuous-time Markov process that contracts arbitrarily faster than its discrete-time counterpart; and (b) a kernel $P$ such that $P^{m+1}$ contracts arbitrarily better than $P^m$. Hence, our main conclusion is that the four standard inequalities comparing five common notions of entropy and variance contraction are generally not improvable. In the process of analyzing the counterexamples, we survey and sharpen the tools for bounding the contraction coefficients and characterize properties of extremizers of the respective functional inequalities. As our examples range from Bernoulli-Laplace model, random walks on graphs, to birth-death chains, the paper is also intended as a tutorial on computing MLSI, SDPI and other constants for these types of commonly occurring Markov chains.Open AccessEntropy Contractions in Markov Chains: Half-Step, Full-Step and Continuous-Time
This paper considers the speed of convergence (mixing) of a finite Markov kernel $P$ with respect to the Kullback-Leibler divergence (entropy). Given a Markov kernel one defines either a discrete-time Markov chain (with the $n$-step transition kernel given by the matrix power $P^n$) or a continuous-time Markov process (with the time-$t$ transition kernel given by $e^{t(P-\mathrm{Id})}$). The contraction of entropy for $n=1$ or $t=0+$ are characterized by the famous functional inequalities, the strong data processing inequality (SDPI) and the modified log-Sobolev inequality (MLSI), respectively. When $P=KK^*$ is written as the product of a kernel and its adjoint, one could also consider the ``half-step'' contraction, which is the SDPI for $K$, while the ``full-step'' contraction refers to the SDPI for $P$. The work [DMLM03] claimed that these contraction coefficients (half-step, full-step, and continuous-time) are generally within a constant factor of each other. We disprove this and related conjectures by working out a number of different counterexamples. In particular, we construct (a) a continuous-time Markov process that contracts arbitrarily faster than its discrete-time counterpart; and (b) a kernel $P$ such that $P^{m+1}$ contracts arbitrarily better than $P^m$. Hence, our main conclusion is that the four standard inequalities comparing five common notions of entropy and variance contraction are generally not improvable. In the process of analyzing the counterexamples, we survey and sharpen the tools for bounding the contraction coefficients and characterize properties of extremizers of the respective functional inequalities. As our examples range from Bernoulli-Laplace model, random walks on graphs, to birth-death chains, the paper is also intended as a tutorial on computing MLSI, SDPI and other constants for these types of commonly occurring Markov chains. - Open Access
Repeated Block Averages: entropic time and mixing profiles
We consider randomized dynamics over the $n$-simplex, where at each step a random set, or block, of coordinates is evenly averaged. When all blocks have size 2, this reduces to the repeated averages studied in [CDSZ22], a version of the averaging process on a graph [AL12]. We study the convergence to equilibrium of this process as a function of the distribution of the block size, and provide sharp conditions for the emergence of the cutoff phenomenon. Moreover, we characterize the size of the cutoff window and provide an explicit Gaussian cutoff profile. To complete the analysis, we study in detail the simplified case where the block size is not random. We show that the absence of a cutoff is equivalent to having blocks of size $n^{Ω(1)}$, in which case we provide a convergence in distribution for the total variation distance at any given time, showing that, on the proper time scale, it remains constantly 1 up to an exponentially distributed random time, after which it decays following a Poissonian profile.Open AccessRepeated Block Averages: entropic time and mixing profiles
We consider randomized dynamics over the $n$-simplex, where at each step a random set, or block, of coordinates is evenly averaged. When all blocks have size 2, this reduces to the repeated averages studied in [CDSZ22], a version of the averaging process on a graph [AL12]. We study the convergence to equilibrium of this process as a function of the distribution of the block size, and provide sharp conditions for the emergence of the cutoff phenomenon. Moreover, we characterize the size of the cutoff window and provide an explicit Gaussian cutoff profile. To complete the analysis, we study in detail the simplified case where the block size is not random. We show that the absence of a cutoff is equivalent to having blocks of size $n^{Ω(1)}$, in which case we provide a convergence in distribution for the total variation distance at any given time, showing that, on the proper time scale, it remains constantly 1 up to an exponentially distributed random time, after which it decays following a Poissonian profile. - Open AccessRelaxation time and topology in 1D $O(N)$ modelsWe discuss the relaxation time (inverse spectral gap) of the one dimensional $O(N)$ model, for all $N$ and with two types of boundary conditions. We see how its low temperature asymptotic behavior is affected by the topology. The combination of the space dimension, which here is always 1, the boundary condition (free or periodic), and the spin state $S^{N-1}$, determines the existence or absence of non-trivial homotopy classes in some discrete version. Such non-trivial topology reflects in bottlenecks of the dynamics, creating metastable states that the system exits at exponential times; while when only one homotopy class exists the relaxation time depends polynomially on the temperature. We prove in the one dimensional case that, indeed, the relaxation time is a proxy to the model's topological properties via the exponential/polynomial dependence on the temperature.Open AccessRelaxation time and topology in 1D $O(N)$ modelsWe discuss the relaxation time (inverse spectral gap) of the one dimensional $O(N)$ model, for all $N$ and with two types of boundary conditions. We see how its low temperature asymptotic behavior is affected by the topology. The combination of the space dimension, which here is always 1, the boundary condition (free or periodic), and the spin state $S^{N-1}$, determines the existence or absence of non-trivial homotopy classes in some discrete version. Such non-trivial topology reflects in bottlenecks of the dynamics, creating metastable states that the system exits at exponential times; while when only one homotopy class exists the relaxation time depends polynomially on the temperature. We prove in the one dimensional case that, indeed, the relaxation time is a proxy to the model's topological properties via the exponential/polynomial dependence on the temperature.